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ANALYTICAL DESIGN OF CONTROLLERS IN SYSTEMS 
WITH RANDOM ATTRIBUTES 
1. STATEMENT OF THE PROBLEM, METHOD OF SOLVING 


N. N. Krasovskii and E. A. Lidskii 


(Sverdlovsk) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 9, 
pp. 1145-1150, September, 1961 

Original article submitted March 18, 1961 





The systems are considered which are undergoing random variations during the control process 
and also are submitted to random interference. The method is studied of determining the con- 
trolling rule such that it ensures stability and also minimizes a given integral criterion of 
quality. The formulation of the problem as well as the approach to its solution are given in 
Part I, the latter being based on the method of the Lyapunov function [1, 2] brought up to date 
to comply with the principles of dynamic programming [3, 4). 


1. Preliminary Notes 

1. Consider a system under control (see the diagram). Here z(t) is the regulating vector quantity, Z(t) its 
prescribed value, x(t) = z(t)—zg(t) the discrepancy vector, g(t) the given input to the controlled device, (¢ + &) 
the output of the regulator, y(t) the interference, (t) the factor describing the random variations of the con- 
trolled object A. 





The system consisting of the controlled object A and the regulator B should produce the prescribed motion 


Z(t) when the given signal g(t) is acting on it. In the absence of the discrepancy [x(t) = 0], this is effected by the 
action ¢ of the regulator. 


The construction of the action ¢ such that it would bring about the motion z9(t) prescribed in advance is out- 
side the scope of the present paper. We shall just consider the quantity ¢ as known. If some coordinates x; of the 
vector x(t) do not vanish, then the member B will, in 


addition, generate a stabilizing signal € regulating 
the transient process, The analytic design of the quan- 





















































ft) wt) tity € [5] will also be our aim in this work. The rule 
1} of control £ (x, 1) shall be determined via minimizing 
e g(t) » a(t) 2 Ezy) i 20) (2), 2,2) a given integral criterion of quality. 
' ‘ aft) =2 y(t) +2(t) We shall assume that the equations of the transient 
ate lave a1) process in the system can be written in the form 
y(t) 
dz 
c } z(t) (9) FM LT In, 7 (¢), El+ % (1,1) 
_i 











E — E [Ziy---) Dns n) (1.2) 


in terms of the coordinates x; of the discrepancy vector x(t). 


The functions ¢; are assumed known; they are continuous, The fact that the main part of the system, the 
member A,can be subjected to random variations during the regulation process is a special characteristic of this 
system. This is taken into account by introducing a random variable 1(t) as one of the arguments of the functions 
¢j- The statistical properties of the random component (t) and of the interference y(t) are described in section 2. 


It follows directly from the definition of the quantities x(t) and € (x, 9) that for the prescribed motion z(t) 
the equalities x = 0 and € = 0 are fulfilled. If the discrepancy x(t) does not vanish at the initial time instant t= to, 
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then a transient process takes place in the system, the process being described by the stochastic equations (1.1); the 
member B produces the signal € following the measuring of the values of n(t) and x(t) which have occurred, the 
member B working as an ideal regulator [6] in accordance with the equation (1,2). This process is a probabilistic 
one [7]. 


2. Properties of Random Functions n(t) and y(t) 

‘Fhe behavior of the object at the time instant tis determined by the values which are — being taken 
by the random function (t). An object with electrical resistances can serve as an example, the resistances being 
functions of a randomly varying temperature n(t). It shall be assumed that 9(t) is a Markovian random process [7], 


The symbol P [L / Q) denotes the conditional probability of the event L when Q has taken place, the symbol 
M {2 / Q} denotes the conditional expectation of the random quantity 1 when Q has taken place. o(e) are under- 
stood to be expressions of higher order of smallness than the small quantity ¢, and O{e) a quantity of the same order 
of smaliness as ¢. The random variations (t) are described by means of the functions q(a) and q(a, 8) in the 


following way: 


Pin(sy) =a(¢<cr<—t+ At)/y() =a) = 1 — g(a) At + 0 (Ad), (2.1) 





Pin(e) < By (te) falt< vw <t+ Ad)/y() =a) = (2.2) 
= q (a, B) At + o (Ad). 


It can be seen from the definition of the functions q(a) and q(a, 8) that q(a, —@)= 0, q(a, ~)= q(a) are 
true. When n(t) can take only one of the k values n = {ay,...., at , then in order to describe fully the process, 
it is enough to know the transition matrix || p;j/] , where 


P [ny (¢ + At) = aj/n (t) = og] = pysAt + 0 (Ad), (2.3) 
k m 
q (ai) om Py Pij, (a, B) = = Pi for Gm <S B < Om+1« 
’ j=1 
A 
If the probability density exists for the function q(a, 6) thatisif g(a, B) = \ p (a, v) dv, then 


es, Bs 
P [Bi < n (t+ At)h< Be, n(t 4+-At)4a/y (t)=a]= At p (a, v) dv+o (At). (2.4) 


BP shall also be assumed that it is possible to measure the occurring values of 1 (t) confining ourselves, however, 
to the case when the transmission of the signal n (t) in the channel takes place without either distortions or delays* 
until 45 teaches the member B(see the figure). 


The interference y at the input will be described by means of random pulse functions in the following way: 

‘the quantities x(t), (t) and € (t) have taken place, then during a small time interval t = r<t+ At in the 
equation (1.1) one has to take into account having in the interference y,; a pulse component with the probability 
p(At) =a. At + o( At); the latter produces a stepwise change Ayx; * vu of the discrepancy where »; is a random 
quantity and uy; a known function. The average value of the random quantity vj is taken as zero; its variance 
M {v= oj = 0 and the correlation coefficients kj; (where M{v4vj} = kjjo joj) are assumed as known. We shall 
[Seen ae oe ening ence Of tho fnteriesence described above A + and oj — 0 in such a way that 
of = const [9]. In the third part of this work when considering linear systems, the assumption of the limiting case 
X - wand oj - 0 shall be dropped. 


* Should there be distortions in the quantity 1 (t) during its transmission to the block B, this will not involve any essen- 

tial complications provided that the distortions at different time instants are not correlated. In fact, if in such a case 

the distortion »-- is described by the distribution function F (B /a} = P(yn’ < 8/7 (/) =a}, one can obtgin from 

the law of large numbers (7) that if the control rule is taken as £ = €* [x, n* ], the process will proceed as if there 
co 


were no distortion but with a control rule ¢ (z, y) = Af [E* (z, n°)/n] = ) t* (2, v) dF (v, 1). 
The latter integral should be regarded as Stieltjes integral [8]. —” 
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3. One can verify® that with random functions n and y as described above, a Markov random process [7] 
takes place in the system (1,1)-(1.2), the process occurring in the phase space {x;,..., Xp» 1}. This signifies 
that the initial conditions Xo, ..., Xno» No» to determine the random functions x;(t) = x4(Xp, Mo, to, t) (1 = 1,.0..M), 
N(t) = 1(NMq, to, t) (t = ty), for t > te, such that the probability distribution of the quantities { x;(t)}, n(t) is not 
affected by the knowledge of the actual behavior of the random process for T < ty. 


3. Formulation of the Problem 


1, Let a function w[x;,..., Xp»1, &] be given, positive definite in the x; and such that a criterion of 
quality of the transient process is defined with its acid 





Ik [Z5, Nol — \M {o/z, = Z19,---» Za = Ings N= No for t=t= 0} dt, (3.1) 
0 


where 


@ = [2 (t),-.-, tn(t), 0 (4), & (2 (4), n(¢))I- 


The problem can be formulated in the following way. It is required to find the function & = (x, wees Xo 1) 
such that if x(t) is the solution of the system (1.1)-(1.2), the following conditions must be satisfied: a) the given 
motion x = 0 is stable in probability [12]; b) the deviations { xj) } which occur in the system will generate a process 
asymptotically stable in the probability [12] with respect to x = 0; c) the integral (3.1) attains for a selected control 
rule & = € (x, n) its minimum value (for all admissible initial conditions x», 19) thatisitisless than for any other 
choice of — (from a certain well defined family of functions {€ (x, n)}). 


Note. The problem is only formulated when the intensity of the interference y approaches zero with the 
discrepancy x vanishing. When this condition does not take place, the process can persist in the neighborhood of 
the given motion x = 0 with an accuracy of an avoidable error A> 0 (Zzj(!)< A*) . In this case the conditions 
a) and b) should be valid when ¢ > A and the average value M{w} under the integral sign in (3.1) must now be 
computed for realizations x(P)(t) such that xi (t)>A* when &SrT<t. 


2. The problem under consideration belongs to the class of optimum control problems, The optimum selection 
of parameters in linear systems based on the Lyapunov functions method was investigated by N. G, Chetaev [2, 13]. 
The problem of least action time was formulated, and the first results given by A. A. Fel'dbaum [14]. The same 
problems were investigated by A. Ya. Lerner in a number of papers (see, for instance,[15]), Some very deep results 
in the mathematical theory of optimum control are due to L, S. Pontryagin, V. G, Boltyanskii, R, V. Gamkrelidze, 
and E. F. Mishchenko [16] whose investigations became the basis for a considerable number of studies (see the studies 
of L. I, Rozonoér [17]), Fundamental theoretical inquiries with regard to the existence of optimal solutions were 
initiated by A. F, Filippov [18], Many authors investigated the sampled-data systems. Interesting results in this field 
were obtained by Ya, Z, Tsypkin [19]. Also the publications of R, Kulikowski [20], R. Bellman [21], R. E. Kalman 
and J, E, Bertram [22], J. La Salle [23] should be noted. The statistical aspects of optimum control were examined 
by L, S, Pontryagin and E. F, Mishchenko [24] on the basis of the maximum principle, and on the basis of the maxi- 
mum rule of dynamic programming*® * by R, Bellman and his colleagues [4], 


* An exact definition of a solution of stochastic equations (1.1) and (1.2) shall not be given here, only the following 
simple and clear interpretation will be used. One assumes that to each realization n(PXt) and y Xt) of the random 
functions n(t) and y(t) there corresponds the realization of the solution x(t) satisfying the equations (1,1) and (1.2), 
with n = n Pt), within the intervals t, <t < t,+1, which are free from pulses entering the realization of y Pt); 
at the instant t = t, of the pulse occurring, the realization of x(P); (t) suffers a jumplike change A,2'7)= vj, (x, &, 9) - 
By giving a full and exact development to the points 1-3 of the section 2 it is possible by using the terminology of 
the theory of random processes [7-11] to verify the existence of a random process x(t) which is a solution of the equa- 
tions (1.1)-(1.2) in an exact meaning of the word. 

** In the latter paper the problems have certain characteristics which make them different from the usual problems 
of optimum filtering, prediction, etc., explored in [25-27]; they also differ from the statistical problems considered 
for instance in the papers [28-29]. 
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The present work deals mainly with the studies of the problem formulated in para. 1 of this section, and it is 
based on the method of the Lyapunov functions [30-31] and on the use of methods of dynamic programming. In 
Part 2 we shall describe how to find approximations to optimum functions and how to solve the problem, This 
procedure shall be based on a continuous deformation of the system via the introduction of a parameter [30]. In 
Part 3 a solution is given to the problem of minimum squared error in the case of linear systems. Studies of the 
optimal integral criteria of quality in linear systems were initiated by A. A. Krasovskii and A. A. Fel'dbaum (see 
[32] and the literature there). The problem investigated by A. M. Letov in its deterministic form is now generalized 
by the approach adopted in our work. 


4. General Approach to the Method of Solution 

We shall assume that a function v(x, . .., Xp» 1) has been found to be such that it will satisfy the following 
conditions” : a) the function v(xj, ..., Xp, ) is positive definite for all —o<x; <o (i= 1,...,M), that 
is there exists a continuous function w(x;,..., Xn) such that v(x;,..., Xp. 1) = W(X...» Xp) > O for all 9 
and x # 0, and also lim w(x;,..., Xp) = @ if 2 2? —» oo ;b) the function v(x, 7) admits a higher limit,” * that 
is v(Xy,..+.%p, 7) = W(xy,..es X,), where the continuous function W satisfies the equality W(0,... , 0)= 0; 
c) the generalized derivative [12, 31] (dM {(v)} / dt)-. when (1.1) and (1,2) are valid, satisfies for E = _°(x, n) 
the conditions | dM {v}/dt | > vk (k = const) and also 





(oe Yee +o [z, 1 E°] = 0, (4.1) 


with the left-hand side attaining its minimum when € = €°, that is 





(“),. + Iz, n, E°)= min [(), + (x, n, &) | , (4.2) 


Remark. In the terminology of the theory of random processes the derivative dM{ V} /dtis determinedby the 
infinitely small differential operator [10] of the Markov process x(t), n(t). This can be illustrated in the following 
way. Suppose the values x(t) = x, n(t)= 1 were realized, From the point (x, 9) there originates a pencil of trajec- 
tories {x (r), n(t)} (r =), generating the random quantity 


fo (e (¢+ Ad), n (t+ Ad) /2z() = 2,9 (0 = OI. 


In order to compute the derivative dM {v}/dt we evaluate the mathematical expectation M (v (z(t + Ad), 
m (t+ At) / z(t) = z, 4 () =), from which we compute v(x, 9), divide the result by At > 0 and proceed to 
the limit with At — 0. 


The function v°(x,  ) satisfying the conditions a) - c) shall be called an optimum Lyapunov function, From 
the results of the theory of dynamic programming [3-4] and from the theorems on stability in probability [12], it 
follows that the function v(x, 7) together with the corresponding function € %x, 1) from the formulas (4.1) and (4.2) 
provide a solution of our problem,*** that is when & = E(x, 7) all the requirements a) - c) of the section 3 are 
fulfilled. Thus the functions v° and £ ° satisfying the conditions a) - c) solve our problem, The actual construction 
of the functions v° and € ° can be found in Part 2 and Part 3. 


The authors wish to express their gratitude to A, M. Letov for discussing the present work with them. 


* Any initial deviations x are assumed permissible. If deviations in some domain G are only allowed, one has to 
modify the subsequent considerations accordingly. 

** A quadratic form v = Db4;x;x; satisfies, for example, the conditions a) and b) if its coefficients are bounded and 
if they satisfy the Sylvester's criterion of constant sign, that is by (n)>e.),;(m) [2 >e,---, lo; (WIT>e, © * 
= const = 0, uniformly in 7. 

*** In (12) stochastic systems were considered of a more specialized type, nevertheless the theorems on stability in 
probability remain also valid in the case now being considered by us. For the unavoidable error A > 0 (see the 
Note in section 3) the requirements a) -c) should be modified for the domain z{> 4’. 
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The paper analyzes the passage of random signals through a time discriminator and an integrating 
amplifier. Based on the analysis performed in this paper, the properties of a pulse system are detez - 
mined and their relationship to the properties of asynchronous detector are demonstrated. For a 
pulse system we find the equivalent circuits with constant parameters in which the processes are 
continuous. Replacing the pulse system with an equivalent circuit simplifies the analysis and 
leads to simple expressions for the characteristics of the original system. These expressions take 
into account the special features of operation of the pulse system with a thoroughness that is ade- 
quate for practical purposes. 

In part I of this paper we derive the recursion relationship in matrix form for determining 
the coordinate lattice functions which characterize the random processes in the pulse system. 


Introduction 

A pulse system consisting of a time discriminator and an operational amplifier plays an important part in 
automatic radar units and radionavigation systems, In particular, such a system is an important element in an auto- 
matic range finding unit for which certain theoretical problems were analyzed in [1-5] et al. The use of frequency 
methods in analyzing a linear pulse control system with time selectors was studied in [6], 


The engineering principles of designing time discriminators and certain characteristics of time discriminators 
have been dealt with in a number of papers. Individual sets of data on these problems are cited in the books [7-10] 
etal.” Our literature does not yet contain papers in which a) the scattered data on time discriminators is system - 
atized and generalized, and b) a complete theory for the operation of these devices in various automatic systems is 


presented. 


In order to clarify the properties of the pulse system under study when random signals pass through it we make 
use of the paper by V. S. Pugachev on the canonical expansion of random functions [13], The use of other methods 
for analyzing random processes in a pulse system with time selectors gives rise to difficulties, since these processes 
are nonstationary by their internal nature. However, these nonstationary processes can be described by a system of 
stationary lattice functions. These provide the possibility of finding equivalent stationary continuous systems that 
approximately simulate the dynamics of a pulse system with time selectors. 


Based on this we developed a method for analyzing the dynamics of pulse systems with time selectors which 
is presented in [14, 15]. The application of this method made it possible to discover many important details which 
are characteristic of the operation of pulse systems with time selectors. The conversion of signals in pulse systems 
with a time discriminator, as we shall demonstrate below, has much in common with conversionin a synchronous 
detector. 





* Cf. also [11] and [12], 
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Problems of synchronous detection, which is also called selective detection, were first studied in [16], Syn- 
chronous detection usually is applied in those cases where higher noise stability requirements are imposed on a radio- 
engineering system. For example, it is widely used in experimental radiospectroscopy units [17, 18]. Synchronous 
detection in these units assures a high sensitivity of the electronic devices. In radioastronomy the combination of a 
highly stable synchronous detector and a narrow-band output unit assures isolation of very weak signals whose power 
is only thousandths of the noise power [19, 20]. At the present time many forms of synchronous detection are known 


which are applied in various automatic systems. In this regard the theory of synchronous detection requires further 
development and refinement. 


We should, however, pay attention to the fact that although the conversion of random signals in time dis- 
criminators and an integrating amplifier has much in common with conversion in a synchronous detector, there are 
substantial differences between these types of conversion; we shall demonstrate these differences in this paper. 


The transient responses in a pulse system consisting of a time discriminator and an integrating amplifier are 
characterized by large time constants; i.e., the specified system has a sufficiently large inertia. Therefore only the 
low-frequency components of the input signal can have an effect on the operation of such a system. In solving 
applied problems and performing engineering computations we usually originate from this fact and replace the pulse 
system with an equivalent continuous system which passes the low-frequency components of the input signal.” In 
reality we encounter the opposite picture. Due to the presence of a time discriminator in the system the high -fre- 
quency components of the input signal produce the principal effect on the system; this contradicts the statement 
made above. In fact, there is no contradiction here. This can be explained physically on the basis of the fact that 
the high-frequency components are converted into low-frequency components in the time discriminator; these latter 
components cannot be filtered by the narrow-band filter. Here frequency conversion in the system is achieved by 


means of linear elements. This is the principal difference between the specified system and other systems with 
identical transient responses. 


Using a time discriminator and an integrating amplifier, we can achieve the isolation and analysis of the 
information which is carried in radar signals, A penetration into the essence of the physical phenomena which occur 
in the pulse system dealt with in this paper will make it possible to design highly effective units for isolating and 
analyzing radar information and transmitting it to other elements in an automatic system, If a time discriminator 
is included in an automatic system as a detecting element, then its properties have an appreciable effect on the 
operation of the automatic system as a whole. In part III of this paper we shall turn our attention a) to the design 
of a continuous system which is equivalent to the pulse system, and b) to the derivation of approximate expressions 
which are simple in structure but consider the specifics of the pulse system in a sufficiently complete manner. 


1. The Mathematical Description of the Dynamics of the Pulse System 

In order to determine the basic properties of the pulse system discussed above we shall limit ourselves to the 
solution of the linear problem. In accordance with this we shall study a system consisting of a time discriminator 
and an inertial section in which the processes can be described by the equations 





(T, (2) D + 1} U (2) eax Ky (LE, Yo (t) + EMI, (1.1) 
[t. D + 1] Ug (t) = K,U, (0). 


Here D = d/dt is the differentiation operator, t is the “running” time, U, is the output signal of the time discriminator, 
U, is the output signal of the operational amplifier which is an inertial section with the time constant rg. (If we in- 
crease the time constant T2 of the inertial section, then its properties will approach the properties of an integrating 
section.) Kg and Ky are the corresponding gains. 


The second equation in (1,1) describes the processes in the vacuum-tube irtegrator whose operator transfer 
function is written as 


RK, 
Ky i (D) = TRH ERD TT 





gititie 
Cf. for example, [21]. 
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where R is the resistance of the integrator input circuit, C is the capacitance of the capacitor in the feedback cir- 
cuit (Fig. 1), From this equation it follows that 


Tt. =(K,+1)CR. 


More detailed data on a vacuum-tube integrator and operational amplifiers is given in [22] 
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Fig. 1. 1) From starting pulse generator; 2) selector pulse generator; 3) 
time discriminator; 4) vacuum -tube integrator. 




















The circuit for the pulse system is shown in Fig. 1. Here Al and A2 are coincidence amplifiers (time selector), 
DD is a difference detector. The coincidence amplifiers Al, A2 and the difference detector DD form the time dis- 
criminator circuit. Each coincidence amplifier passes only that portion of the input signal which coincides in time 
with the corresponding selector pulse. The selector pulse generator forms two pulse trains that are time shifted by 
the amount op» with respect to one another, The difference signal obtained at the output of the time discriminator 


is applied to the integrator CR based on the tube T. The operation of the time discriminator is described in greater 
detail in part I [5]. 


The useful signal E,yo(t) and the noise E,(t) are applied to the input of the system. We shall denote the in- 
put signal by 


Ein (t) = & Yo (t) + EF nil). (1,2) 


The first equation in (1.1) describes the conversion of the input signal E;, (t) in the time discriminator. Ac- 
cording to this equation the input signal is first subjected to time conversion described by the pulse function v(t), 
The graph of this function is shown in Fig. 2b. After time conversion the signal is subjected to additional conversion 
in a fourpole with a variable parameter 1, that varies stepwise. The second equation (1.1) describes the subsequent 


conversion of the signal in the inertial section with the time constant rT2 after the signal has passed through the time 
discriminator. 


The pulse function y(t) is written as 


Volt) = D>) [H(t — me —kT) — H(t —2 — me —AT)). (1.3) 


k=0 


Here a is the duration of the pulses which form the useful component, T is the repetition period for the reference 
pulses, H(t) is a unit switching function which satisfies the conditions 


eG 5 for t< 0, 
sdeosess T for t>0. (1.4) 
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The lattice function 1, (k= 0, 1, 2,3,..., #) describes the time position of the kth pulse of the useful 
component relative to the reference pulse with the number k, 


The pulse function y(t) which describes the conversion of the input signal in the time discriminator can be 
represented by the following formulas: 


¥@ =n () — 72 (0). (1.5) 


vi () = >) LH (t — 0, — kT) — H(t — ay —%& — &T)), 
k=o (1.6) 


12 (t) = 2; [H (t — a — % — kT) — H (t — 209 — 0 —k7)). 
k=0 


Here oo is the duration of the tracking (selector) pulse, The lattice function %,(k= 0,1, 2,..., ©) deter- 


mines the time displacements of the tracking pulses relative to the reference pulses, The tracking error in the case 
under study is determined from the equation 


Aa [n} = tn +5 — On— (1.7) 


We note the fact that in our papers [5] the pulse functions wo(t), y(t) and y2(t) are written in a different 
form. The difference is very slight and consists of the fact that the time origin in this present paper is chosen in 
such a way that it coincides with the leading edge of the reference pulse with the number 0. The form of notation 


for (1.3) and (1.6) proves to be more convenient in studying random processes in the dynamic system under study 
here. 
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Figure 2 shows graphs of the pulse functions Ey g(t), y(t) and Egy oft) y(t). The quantity 9, determines the 
position of the leading edge of the tracking pulse I with the number n relative to the reference pulse with the same 
number. Analogously, the quantity 1, determines the position of the useful-component pulse with the number n 
relative to the corresponding reference pulse. 


Figure 3 shows the graph of one realization of the random function E,(t) which describes the noise at the in- 
put of the time discriminator. The same figure shows the graph of the function E,(t)y(t) which expresses the trans- 


formation of the specified realization of the random process Ep(t) by the system of time selectors (coincidence am-~- 
plifiers). 
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Figure 4 shows the block diagram of the pulse system. The time discriminator is represented by 2 serially con- 
nected fourpoies which are inscribed in the dotted rectangle. The conversion of the input signals in the time selector 


and the difference stage is represented by a fourpole with a variable gain Koy (t). The subsequent gain of the signal 
in the circuits with varying parameters is represented by a fourpole with the operator transfer function 


1 
Q,(D, t) = muy D+i : (1.8) 


Since a continuous signal consisting of the useful component and the noise appears at the input, it follows that 


the parameters of the time discriminator will vary in accordance with the application of the selector pulses to the 
system; i.e., 


T (Q) = to Go (t) + tgs (2), (1.9) 
PM=-nd+rwr, 2) =A) —-e id. (1.10) 
In expanded form the pulse function is written as 
@ (t) = >) 1A (t — Oy — kT) — H (t — 205 — Oe— &7)I. (1.11) 
k=0 
Figure 5 shows graphs of the pulse functions y(t), g(t) and ¢;,(t). 
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We shall assume that the noise component E,(t) is described by a stationary random function for which the 
integral canonical form is written as° 


* Cf. [13], pp. 315-323. 


1030 


























+00 
E ,(t) = met+ \ Y, (a) edo, 


pn (1.12) 
ke (t—t)= | Sa) emterda, 





K.(@,@') = M lY.(@) ¥-@’)] = & @) 6 @ —@’). 
Here k,(t—t’) and Se(w) respectively denote the correlation function and the spectral density of the noise E,(t). The 
constant Me is not passed by the time discriminator; the output signal which it produces is equal to zero. Therefore, 
it will not be a simplification if we assume that the component mg in the first of Eq. (1.12) will be equal to zero;i.e., 
= M [En(t)] = 0. (1,13) 


We shall divide both sides of the first equation in (1,1) by 1,(t) and the second equation by rg. As a result 
we will obtain: 


ID + a (9) Us (Q) = “2 [Beye () + En), 


(1,14) 
a2,U; (t) + (D+ Qo) Uz (t) = 0. 


Here 


1 
4a; (t) = = aw 9. OS + A ee (1.15) 
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The right side of the first equation in (1.14) includes the constant Ty together with the variable r,(t). This 
can be explained on the basis of the fact that the right hand side of this equation is non-zero only over time inter- 
vals during which pulses which can be described by the function ¢go(t) exist. During the intervals between pulses 
Y(t) when the constant 1,(t) acquires the value p(t) the right side of the first equation in (1,14) goes to zero. 


In order to abridge the notation and achieve compactness in our presentation we shall write the equations in 
matrix form 


[DI +a()) UW) = cay (t) + E ()). (1.16) 
Here 
,(), 0 U, (t 
a (t) -|""" we |. U (t) -| ph * . (1.17) 
E .; 
f(t) = —? . £0 =| "|. (1.18) 
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The quadratic matrix a(t) can also be written as 


a (1) = a%po (t) + a’, (0), (1.19) 
where 
4 1 
o a) “1 0 1 -| a, 0 |. 
a G22, |. r Qo, Ase (1.20) 


The coefficients a}, and al; are defined in terms of the time constants T 49 and Ty, using the formulas 


eh 


a 
Tio * 


| 
, oe oe (1.21 
ay Ti ) 


1 
ll 
2. Equations for the Coordinate Functions which Characterize the Random Process in 
the Time Discriminator and the Operational Amplifier® 

We introduce the coordinate functions for the time discriminator and the operational amplifier: 4,(t) and 4(t), 
respectively. Based on (1,14) the equations for #,(t) and 4(t) are written as 








ID + ay (1, () == 7 eo, 


(2.1) 
az,D, (t) + ID + ase] D2 (t) = 0. 
The initial conditions for Eq, (2.1) are written as 
®; (0) = M2 (0) = 0 for t= 0. (2.2) 
We shall write Eq. (2.1) in matrix form: 
IDI +410) == ¥@,9 10 (2.3) 
The initial conditions are 
® (0) =0 for ¢=0. (2.4) 
Here 
® (t) = | ‘ i (2.5) 
iat 
Y (a, 2) -| 0 |. (2.6) 
We shall find the solution for Eq. (2.3) over the time interval 
nT <t<(n+1)T. (2.7) 


The time interval (2.7) is subdivided into three regions: 


Region J: n7 <i <q nT +9,, 
Region JJ: n7? +0, <t <Q nT +O, + 200, (2.8) 
Region IJIl: n7 +0, + 200 <t < (n+ 1) 7. 


* The general method for determining the coordinate functions for pulse systems and its application to individual 
particular cases was presented in detail in [15]. 
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We shall find the expression for the matrix function 4(t) in each region, 


3. Determining the Coordinate Function #(t) in Region I 
Equation (2.3) in region I is written as 





IDI +a] O() =0, nT <tc nT +%p. (3.1) 


We shall write the initial conditions for the instant t = nT: 
® [n] = ® (n7). (3.2) 
Solving Eq. (3.1) and taking into account the initial conditions (3.2), we obtain 


@ (t) =e" 70 [nn], 











(3.3) 
: nT <t < nT + On. 
4. Determining the Coordinate Function (t) in Region II 
We shall write Eq, (3.2) for the region II: 
[DI + a] ® (t) = fe Y (w, t) LH (t —0, — nT) — 2H (t—co— 0, —n7)}, (4.1) 
nT +0, <t < nT +0, + 200. 
The matrix column Y(w, t) is written as 
ele! (4.2) 
real 
The initial conditions for the instant t = 9, + nT are found from Eq, (3.3): 
® (n7 + 0,) = *"*"@ [n}. (4.3) 
The solution of Eq. (4.1) for the initial conditions (4.3) is written as follows: 
D (t) = enn [n] + W (0), (4.4) 
nT + On < t Cal + On + 2a0. 
The function 
Ww 
W (t) -| w, ~ (4.5) 
is the response of the pulse system to the input signal Y(w, t), This function satisfies the equation 
| [DI + a°] W (Z) = Ko (@. ) 1H (¢ —0, — nT) — 2H (t — co —%,—2n7)), 
| Tio (4.6) 


| nT +0, <t <nT +0, + 200 
' for the initial conditions 
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W (nT + 0,) = 0. (4.7) 


In order to solve Eq, (4.6) for the initial conditions (4.7) we shall write the function Y(w, t) in a different form: 


eiul 
Y (o, 4) -| 0 |- ef2'Vo, (4.8) 
where 
eet 0 
eat — | 0, eiut |. (4.9) 
Yo -| ' (4.10) 
0 . 


We shall rewrite Eq, (4.6) in new variables: 


IDI +a) W(t) = at e''Yo [H (¢ —0, — n7) — 2H (t — a, — 6, — n7)), (4.11) 


nT +0, << t < nT + On + 2a. 


The solution of Eq, (4.11) for the initial conditions (4.7) is written as 


K,c12T n A2%, 





W (t) a SAFE oF izle doa ae peerer Po] 
(4,12) 
es 2 eifit { elite 8, T) — een Y, 
or 
W (t) = efO Men Z (1), (4.13) 
where 
a K. {2(i—9,—nT ) —a"(i—8,—nT) 
Z(t) t co a’) {le ang sae ] ag al 
See a) ee ed oO T))) Y,. (4.14) 
5. Determining the Coordinate Function 6(t) in the Region III 
We shall write Eq, (2.3) for the region III: 
IDI + a’) ® (i) = 0, nT +, + 200 < t < (mn + 1) 2 (5.1) 


This equation coincides in form with Eq, (3.1) for region I, but it differs from it in its initial conditions which 
for Eq, (5.1) are written as 


wo yore Z,; (nT + 0, + 2a0) (5.2) 
ah 2a%a, —-a'd. igo(nT+8,,) 9!) 
@ (nT + 0, + 200) = ¢ e” "® [n] + e Zy (nT +m + 2ao)h” 


The quantities Z,(nT + 9, + 2a) and Z,(nT + 9, + 2c) depend on the parameter i w in accordance with 
(4.14) and (4,9), Thus, we can introduce the substitutions: 


Ri (i@) = Z, (nT + On + 200), Re (iw) = Ze (nT +O, + 20) (5.3) 
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and 


Ri (io) 
R (io) = Rs (iw) |" (5.4) 


We shall write the expression for the initial conditions (5.2) in terms of the new substitutions: 


OD (nT + Oq + 200) = cen [n} + c"T +n R (iw). (5.5) 
We shall find a solution of Eq. (5.1) for the initial conditions (5.5) 


@ (t) o—_ et 2ae OT) 20" 28 nt [n] + ciinT +4, Vga —2a—8 nT) R (iw), 


(5.6) 
nT +0, + 200 < t < (n + 1) T. 


6. The Recursion Relationship for Determining the Coordinate Lattice Function which 
Characterizes the Random Processes in a Time Discriminator and an Integrating Am- 
plifier 


Assuming that in Eq. (5.6) we have t = (n + 1)T and writing ®[n + 1)= (nT + T), 








e—C§aT — c—O'(T—20,— 8p) p—20° 2,9 —0'8,, . (6.1) 


we obtain the recursion relationship 


®D [n + 1) = "" @ [nn] + eT Hn eT en (i'w), (6,2) 


from which it is possible to find the coordinate lattice function #[n]. 


The recursion relationship (6.2) can be simplified on the basis of the following concepts. For a time dis- 
criminator and an operational amplifier the following conditions are always satisfied: 


|2ajao|<li, fauT|<4, (6.3) 
where the vertical lines denote the absolute values of the following products: a) the elements of matrix a® mul- 


tiplied by the duration 2d, of the tracking pulses, and b) the elements of the matrix a! multiplied by the repetition 
period T of the reference pulses. 


When inequalities (6.3) are satisfied the following relationships are valid: 
= I — 2a%a,, 
en ~ I — aD, (6.4) 
eu tee) — J — a! (T — 0, — 2a,). 
Therefore 
CaT =~ aT = 2a°a, + a! (T — 2a), (6.5) 


i.e., when inequalities (6.3) are satisfied the matrix C,(n= 0, 1, 2,3,..., ) is practically independent of the 
number n, Taking this into account, we rewrite the recursion relationship (6.2) in the form 


® [n + 1] =e" @® [nn] + eR (iw). (6.6) 
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The resulting recursion relationship for the coordinate lattice function ¢[n] will be used in the subsequent 


parts of this paper to determine the correlation functions and spectral densities corresponding to the output signals 
of the pulse system. 
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OPTIMUM FILTER DISCRIMINATION OF TELEGRAPH SIGNALS 


R,. L. Stratonovich 


(Moscow) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No, 9, 
pp. 1163-1174, September, 1961 

Original article submitted March 27, 1961 


Equations determining the a posteriori probabilities of telegraph pulse signal values for the case 
of “white noise” interference are derived. The article also describes the block diagrams of 
optimum receivers that are able to perform optimum nonlinear filter discrimination and to 
provide filtered-out pulse signals at the output. 


The theory of Markov conditional process [1, 2] can be successfully used in many practical cases for the 
construction and calculation of optimum filter-discrimination systems. Some of the applications were described in 
(3}. The present article is concerned with the application of this theory to another important particular case, namely, 
to the case of the filter discrimination of telegraph signals. 


A telegraph signal should be understood as a two-position pulse signal, i.e., a signal which assumes the values 
s(t)= 4a. Various assumptions can be made with respect to the statistical properties of the pulse durations and the 
moments of transition from one state into another. The durations of pulses and intervals can be dependent or in- 
dependent, they can have different distributions, and, in particular, they can be determined (nonrandom durations). 
Without making our considerations very general, we shall analyze only the case of independent durations of pulses 
and intervals, 


We shall note, however, that different generalizations of the method are possible in the case of dependent 
durations of pulses and intervals that are, however, correlated in a Markov chain, 


Assume that a mixture of the signal and noise r(t) = s(t) + n(t), which is to be subjected to optimum processing, 
arrives at the receiver's input. We shall consider a disturbance n(t) in the form of white Gaussian noise, which has 
the autocorrelation function nn, = N5(r). The methods for generalizing these results for the case of correlated 
noise are given in the Appendix. 


1, Equations for the A Posteriori Probabilities 
We shall begin with the derivation of some general equations for determining the a posteriori probabilities, 
which we shall need later. 





We shall denote by ..., Ty, Tz... . the sequential instants of time when the telegraph signal s(t) passes 
from the state —a into the state +a. These instants of time will be called “points of upward transition.” . . . 04, 
O2,..., we shall denote the sequential instants of time of the reverse transition, i.e., the “points of downward 
transition,” If we choose a certain conditional initial instant of time ty = 0 and if we count the transition points 
from this instant, we have: 


0<t% <0, <...< Te <Oe <..-. 


It can be considered that this process begins with the instant of time t = ty (if necessary, the change of limits 
t)-> —e can be performed after solving the problem). 


Let Wor(Ty, G4. ~~ +» TMe 9x4) be the common a priori distribution density for a sufficiently large number 
of sequential transition points. If the input process r(t) is observed during the time interval 0 < t < T, then, ac- 
cording to the reverse probability formula, the a posteriori distribution density for the indicated transition points 
will be 
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W 0 (t;, G1, «ee Tu? ou) = 
T 


= Kw,, (1,91, ..., Ty» Oy) OXP {— sx Ir (t) — s(t)P dt}, (1) 
0 
where K is a constant which depends on T and r(t). 
In deriving Eq. (1), we used the well-known expression for the functional of the white noise probability: 


T 


W In (t)) = exp {—sy\ n? (t) at} : (2) 


where n(t) = r(t)—s(t) was substituted (see, for instance, [4]. 
We shall now calculate the a posteriori probability that, at the instant of time t = T, the signal will be in the 
upper position: 
w, (T) = P{s(T) = a|r(),0<t< 7), (3) 
and also the probability of the opposite event: 


w, (7) = P{s(T) = —a|r(y), OK t <7). (4) 


If a pulse exists at the instant of time t = T, this pulse can have different ordinal numbers. Let P;, be the 
probability that the kth pulse will arrive at the instant of time t= T. Then, it is obvious that 


M 
w, (T) = pA Py. (6) 


k=1 


Here, M is a certain sufficiently large (nonrandom) number, the exact value of which is not important for 
theoretical considerations, since it does not figure in the final results. The probability that the number of pulses 
will exceed the M value is assumed to be negligibly small. 


The a posteriori probability P;, can be found from the distribution density (1) by integration with respect to 
all possible positions of the transition points: 


Py = \ saa ( Wps (Tis Or, - » +» Tay Om) At, . . ., doy. (6) 
t%}%<T, o,>T 


Therefore, according to (1) and (5), we obtain: 


M T 
Ww, (7);= K > | “er { Wpr(T, _ -sOm) exp {| za (r—sjtat| x 
k=) ty<T %>T 0 
x dt, — ., d0y. 


However, for o},> T, the integral \ (r—s)*dt isindependent of Tr41,%%}1,-+-Om , 


9 


\. ‘ {woe (Ti, . . +» Om) UAP whe (r — s)? at| dress . « Oy = 


T 
= Wor (Ti, . « «» Tx) EXP {- ow \ (r — s)? at\ 
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and 


M 
wi (T) = KS) \ +? \ Wr (ei «= 0) exp {— shy (r = aa x 


k=itg<T.  o%>T (7) 


x dt, 2 «pi lee 


Similarly, we find the a posteriori probability that an interval between pulses will exist at the instant of time 
t=T: 


M T 
w, (7) = K >} \ es \ wort «sw exp[— shy { (r— Path x 
k=] y—4< T t\%>T 0 


(8) 
x dt, ee e, dt. 


We shall now use more detailed a priori data on the statistical distribution of the transition points. If the 
number pairs (ry, 04), i= 1, 2,..., form a Markov chain, then, 


Wor (Tr, On, « - «5 Tey Oe) = Wor (Ty, Fy) Wor (Te, Op | Tis Oi) - - - 
9 
- «+ «Wor (Tx, Or | Tk—1» Ox—1)s ” 


where Wpr (Ty, 04) is the initial distribution, and wWpr(T4, 04 | 74 - 4» 04 -y) are the transition probabilities. If we 
assume that the durations of intervals and pulses are independent, the sequence of the transition points T,, 04, T2, 
02, .. . Will form a simple Markov chain, and the following expressions will hold: 


Wor (Ty, + ++ Tey Oe) = Wor (11) Wye (O, | T%) . . . Wr (Tr | Oxy) Wye (Ox | Te), 
(10 
Wpr (T,, *“* #9 Tr) = Wor (t,) Wpr (0; | T;) “ee Wpr (Tx | Ox), ) 
where Wpr (1) is the initial distribution, and wp; (0; | 3) and Wpr(T1»94 - 4) are the transition probabilities. 
In the equations 
Wpr (O4| %) = Pr (04 — Ta), Wor (| O11) = P2 (% — F-4) a1) 


they are expressed in terms of the probability distribution density p;(T,,) for the pulse duration Tpu and a similar 
density p, (T;,,) for the intervals. In the case of steady -state signals, de fenctions Pa(Tpy) and p2(Tj;,) are in- 
dependent of t; in the general case, such a dependence can be taken into account by considering the functions 


Wor (64 | Ta) = Pr (O14 — Ta, Ti), Wpr (| Or1) = Po (Ts — O41, O41) 


Limiting our consideration to the case of a simple Markov chain, we shall substitute (10) in (7) and (8). If 
we denote 


M 
Ve (T%) = > 4 a \ Wpr (T) Wor (0, | T)) . « Wpr (Te |Ox—s) X 
k=l O¢_4<ty 


™ 
xX exp \—s¥7 \ (x — s)? at} dx, — —- (12) 
0 
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M 
Vy (Ox) = >) \ sb \ Wpr (41) Wpr (G; | T2) - . - Wpr (Gea | Te-1) X 
Kol tye 4S Op} 
*~—1 


x oxp|— sy \ (— 9) at} dt... dtr, 


0 


the results of the substitution will be given by 


T 
w, (7) = K \ { Vg (Tx) Wpr (Ox | Te) exp | — on \ (r — a)? at dvd. 
Tk 


te<T , oy>T 
(13) 
. T 
w,(T) = K \ { Vy (Ox—1) Wr (Te | Ox—1) EXP [— oN { (r + a)* at| x 
Se-1<T , %e>'T oe-1 


4 do,_,d tr. 


Equations (12) determine the v; 2 functions in terms of the arguments rt = rT; and o = 0}. 3. If we substitute 
r and o for the arguments of these functions, we readily see that they do not depend on k, since the summation 
was performed with respect to k, The first of the equations (12) can be written in the following form: 


- 


t M 
a(t) = (do 5) 0. . (wor (11)... wpr (6 | tes) Wpr (4 | 0) x 
0 k=1 tp <6 


x exp|— oN w\tr— s)* at} ar, ae ae 
Hence, by using the second equation, we have: 
v(t) = (>, (6) Wpr (t | 0) exp {— oy \(r + a)? at| do. (14a) 
Similarly, we obtain: 


v, (0) = \ 2 (Tt) Wor (0 | T) exp |—j aw \( —a)* at| ax. (14b) 


cow" a 


Thus, the summation in (12) is eliminated. 
The latter equations serve for determining the functions v; and vz, which can now be used for calculating w; 
and wa. By introducing the density probabilities (11) and also the integral distribution functions 


Py (2) =(p, (at, Pa (z) =\ ra (at, as) 
the basic equations (13) and (14) can be written thus: 


v, (t) = (>, (¢ — t) exp Lt (cr — a)? ar' v, (t) dt, 
0 t oan 


t t 
ve (t) =\ps (¢ — t) exp —sy\r + a)? dt’\ he (t) dt, 
0 


t 
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a) 


ib) 


5) 


6) 


T T 
w, (T) = K\ P, (T — 1) exp - oy \ (r — a)? a’ Vv, (t) dt, 


a 


(17) 


r r 
w, (7) = Kk \ P,(T — thexp — sw \ (r 4 ay a’| v, (t) dt. 


If we explicitly consider the initial instant of time tg, the following initial conditions must be added to Eqs, (16): 


V, (0) > Wp (0), V2 (t) + Wp, (t) for 9, T—> bo. 


However, under steady-state conditions, the role of the initial conditions is not essential, and they can be left 
out of the consideration, assuming that the lower integration limit is tg = —« in (16) and (17). 


2. Block Diagram of the Optimal Receiver 








The form of Eqs. (16) and (17) determines the block diagram of the receiver that is to perform the optimum 
nonlinear transformation of the received signal r(t) for the purpose of separating the useful signal. The v,,2(t) and 
W1,2 (t) functions can be realized in the circuit as variable electrical signals. Instead of the v,(t) function, it is 
convenient to consider the functions 


i 


v4.2 (t) exp {ay \ (7 + a?) ar’ 


or the functions 


{ 
Vi2 (t) = 2 (t) exp in \(r + a)? a’| ° (18) 
9 
If we introduce the notation 
: 
E() =24\r(t) ae, (19) 


Eqs. (16) and (17) will assume the following form: 


V; (ft) = e& (P: (¢ — t) e~ POW, (1) dt, 


0 


V; (t) = V2 (t — 2) V, (2) dt, 


(20) 
w, (t) = K’e®\ P, (t — 2) eB, (1) dr. 


9 
t 


ws (t) = K’ (7, (t — t)V, (x) dt. 


n 


Figure 1 shows a block diagram which corresponds to the adduced equations, The device which produces a 
E(t) or e*(®) signal comes first. Then there follows a system of units with feedback, which produces the V,(t) and 
V2(t) signals. This system consists of linear units with the transfer functions p;(t—1r) and p,(t—r) and of units which 
perform multiplication by the functions e&(t) and e~©(*), In correspondence with Eqs. (20), the V(t) and V2(t) signals 
are then transmitted to units with a similar structure, but without feedback, which correspond to the transfer functions 
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P,(t—r) and P,(t-r). If the signal value is estimated with respect to the maximum probability, the difference be- 
tween the obtained functions must be fed to a threshold circuit (relay) P, The output signal so(t) = a sgn [wy(t)—we(t)] 
actually represents the filtered-out useful signal. In this, the K' factor, which can always be eliminated by using 
the normalization condition w;(t) + we(t) = 1, is of no importance. Instead of V,(t) and V2(t), we can also obtain 
functions that are proportional to these signals by choosing suitable amplification factors. 











Fig. 1 


The totality of the p; and p2 units constitutes a closed-loop system with positive feedback. Self-oscillations, 
which are modulated by the e*® signal, will be excited in this system. Such conditions are normal for the operation 
of this system. In order to avoid an unlimited rise of the oscillations which occur in linear systems, it is advisable 
to provide a nonlinear element, the gain regulator GR, which limits the amplitude. It is controlled by a certain 
signal | v, | (or | vz | ), which is averaged over several periods. 


It is understood that various receiver modifications can be provided for different actual cases in dependence 
on the technological conditions. 


3, Filter Discrimination of Specially Shaped Telegraph Signals 
Let us consider a telegraph signal with a somewhat more particular shape, for which the above results can be 
modified without affecting the basic principles. 





Assume that the signal transitions from one state into another can occur only at the points t,, = (m—1)Tg+t, 
where m = 1, 2, ... is an integer, and Ty is a fixed quantity (period) that is assigned in advance. The precise posi- 
tions of the transition points are, however, not known, since the value of t, (0 < ty = Tg ), which coincides with the 
coordinate of the first of the indicated points, is not known. 


The transition probabilities entering (10) are now different from zero only at the points ty,. In this case, al- 
though the above-presented theory can be strictly applied, it is, however, more convenient to use it in a somewhat 
different form. 


Each of the periods is characterized by a constant value of the signal sp = s (t)= 4 a, tr) -y. We shall assume 
that the quantities s), , . . . represent a Markov chain with the a priori transition probabilities w,,, (sj | sj - ). Con- 
sequently the multidimensional a priori distribution of signal values during the first m periods will have the fol- 


lowing form: 


Wor (8. - - «> #m4i) = Wpr (8;) Wr (82 | 8;) - «Wye (S41 | fm)- (21) 


We shall assume that the value of t; is independent of s;, %, . . . and that it has the a priori distribution den- 
sity Wor(ts)- 


We shall find the a posteriori distribution Wor (tm» $4, - » » » S44) Of the random quantities t,, = t, + mT9-Ty 
$4, » + + » Sm 41 that corresponds to the observation in interval 0 < t < T, while the end of the observation time and 
tm pertain to the same period (m—1)T» < T < mTp». 


If, as before, the noise is assumed to be white noise, then, from (21) and by analogy with (1), we have: 
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on 


si- 


\- 
at 


1) 


-Tp 
and 





Wie (Ly Sty ~~ +) Smt alt = AWpy (tm — mTo + To) Wor (s;) Wr (Sa | 5). . . 
T 
. «Wor (Sm4i| Sm) Exp {- aw \ (r — s)? at} (22) 


If we are interested in the signal value during the mth period or during the neighbouring periods, the sum- 


mation in (22) must be performed with respect to the remaining signal values s, ®,... . It is convenient to in- 
troduce the function 


D (tm, Sm) = >) Wor (tm — MT'o + To) Wor (8;) Wr (Se | 8) - - - 
Bigrers &m—1 
. + Wr (Sm | Sm—1) exp |— ay ( (r — arat. (23) 


By means of this function, the probabilities 


Wys (tiny Sm) = b> Wye (bane Sr + + +» Sandals 
My m—mt, Smt 
Woe (tins Sm4 ir = > > Ws (tm, Sing + + +9 Sm4i)s (24) 
Eigeces Say 
can be written as 
: T 
Wys (tm, Sm) = K er Wr (Sm | Sm—1) exp =k ( (r — Sm)* dt} x 
m—1 'm—1 (25) 
x V (tm—1, Sm—1) for tm >T, 
T 
{ 2 atl 
oe (bs tts) = K Stop (Sms 8m) OP {— say Cr — sad? dtl (tn 8 
7 ” (26) 


for tm <7. 


In order to construct an optimum discriminating system, it is necessary to know the a posteriori probabilities 
that the signal will have a positive or negative value at the instant of time T. Let us denote these probabilities by 
w(T, a) and wT, —a). 


In view of the fact that the points t,, and T pertain to the same period (m~1)T, < t < mTp, their sequence 
can have different orders. In correspondence with two different possibilities (t,, > T and t,, < T), the above proba- 
bilities are composed of two components: 


mT, 


T 
w (T, 8) = \ Lwpe (tm, Sm)rIaq etm + | UWdpe (ty Smt) le gle. on 
(m—1)T, 


By substituting (25) and (26) in this expression and by combining both terms, we find: 


T T 
w(T, 8) =K >) | woe (s| 8’) exp {= alr} (r — a at! » (x, s’) dt. (28) 
s T—T, : 
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Thus, the probabilities in which we are interested are expressed by the function (23). Let us introduce the 
equation which serves for calculating the latter function. If we increase the number m by 1 in (23), we obtain: 


D (tmii» Sm4i) = >» Wpr (tm+1 — MT 0) Wpr (8;) Wr (S2| 8) . - - 


pe oho ‘ ‘mH (29) 
- + Wpr (Sm4i | Sm) OXP | oN \ (r — 8)? dt —sy \ (7 — Sma) dt . 
0 bn 


By comparing (29) and (23), we arrive at the following relationship: 


(may 


V (tm+i» Sm+i) = eXP {| * \ = smti)*ath x 
‘maT 
. (30) 
~ >Wpr (Sm+1| Sm) ¥(tm41 — To, Sm), 


5m 


which is the usual relationship in the theory of conditional Markov processes. This equation is similar to Eq, (14), 
and it serves for the step-by-step determination of the v(t, s) function. 


Equations (28) and (30) indicate what optimum transformations of the input signal must be performed in order 


to obtain the a posteriori probabilities w(T, a) and w(T, —a) at the system's output, which are to be used for deter- 
mining the filtered-out value s(t) of the signal at the instant of time t = T. 


It should be borne in mind that the a posteriori information will be more accurate and that the filter dis- 
crimination error will be smaller if we know the behavior of the signal r(t) for sometime after the moment t= t 
for which the useful signal value is determined (t® < T). _In this case, the filtered-out signal will have a certain 
time lag T-t° with respect to the received signal, which is compensated by an improvement in the quality of fil- 
ter discrimination. In the case under consideration, it is convenient to choose Ty as the time lag. Let us derive the 
optimum equations for filter discrimination for this lag. We shall find the a posteriori signal probabilities at the 
instant of time which precedes (by one period) the end of observation. If we denote these probabilities by 


wie (1 — To) = P{s(T — Tr) = ta|r(t),0<4+<7), (31) 
we find by analogy with (27): 
mT,—T, T 
W1,2 (T — To) we \ Wps (tm—1, Sm-1) dtm, + \ Wps (tm; Sm) T dtm = 

T—T, mT,—T. 

“re ( (32) 

= \ Ws (tm, Sm—1)T dtm + \ Wps (tm, Sm)T dtm, 
4 mT,--T, 


where Sp -1= 44, Sp = 24. 
From (22), by summation with respect to all s with the exception of sp -1 and s,,, we readily find: 


T 
Wps (tm, Sm) = Kv (tm, $m) 2 exp |— ay \ (7 — Sm-+1)* dt ie (Sm+i|$m) for tm < 7’, 


#m-+1 tm 
; : (33) 
Wy (tims Sm—1)'? = Kv (tm — To, Sm—1) 7 exp | — oN \ (7 = s,,)* aut X Wpr (8m|$m—1) for tm>T. 
*on tmn—T 
By substituting the latter expressions in (32) and by combining both terms, we find: 

T T 
wiz (T’ — To) = K ( v(t, +a) Sexp\— aw \ (r — s)* dt, wpe (8| +.) dx (9 

T—T, ad 1 
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The obtained equations (30) and (34) can be used for solving the problem of optimum filter discrimination. 
We shall apply them to the important particular case where the telegraph messages during the neighbouring periods 
are a priori independent and equally probable. This means that 


Wor (+a}a)=+, Wr (ta|—a)=, Wye (+0) =>. (35) 


For this, we obtain from (30) and (34): 


t 
2 (t)=vi(t, +0) = Sexp|— \ (r ¥ ay? dt") (oy (tT) + v, (t — To)}, (36) 


t"T, 


T T 
Wi (T — To) = 5K \ Vive (T) | exp|- ow \ (r — a)* at} + 


Tents 


T 
+ exp {— aw \ (t -+ a)? at] dt. 67) 


As before, we shall introduce the functions V;,» (t) and E(t) by means of (18) and (19), Then, Eqs. (36) and 
(37) will assume the following form: 


V(t) = seta B79 [Vy (t— Ty) + Valt —Ty)I, 


i 7- . 
Vo(t)= 5 Vi (t— To) + Velt—T oh, @8) 
T T 
w,.2(1 —T,) = Keb \ V-EOW, dt + K’ \ V,.2(t)dt. 
T—T, T—T, 


These relationships represent variants of Eqs. (20) in application to our case. They can be used for constructing 
the optimal receiver whose block diagram is shown in Fig. 2. Instead of linear units with the transfer functions p;(t—r) 
and p2(t—r), the diagram now contains lag units LU with the lag time Ty, Instead of units with the transfer func- 
tions py(t—7) and p2(t—T), we now have units with the following transfer function: 


ae” We G, for O0<t—rt<T,. 
ve o={5 for t—1<Oandt—t>T7,. 


(39) 
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Moreover, in comparison with the circuit of Fig. 1, the circuit of Fig. 2 is able to perform additional sum - 
mation and subtraction operations, The circuit contains four units for multiplication by exp (+ E); however, by 
modifying the circuit in the manner shown in Fig. 3, it is possible to reduce 
the number of such units to three for the same number of other operations. 
Se s,(t-) It is understood that other variants of the circuit can be obtained, 





As before, the signal which is proportional to the difference w;(t—T »)— 
—Wz (t~T9) between the probabilities is fed to the relay which supplies the 
filtered -out signal s9(t-T ) that has been delayed by one period. 





Rit- 
In conclusion, we shall consider the case where the telegraph signal has 
Fig. 4 a high-frequency filling, i.e., where it can be written in the form 
S(t) = s(t) cos (w,t + 9%), (40) 


where s(t) is a signal of the type considered above. In a slightly more complicated form, the theory presented above 
can also be applied to this case. We shall indicate the direction to be followed in solving this problem in the case 
where the frequency w», and the phase 9p, are sufficiently stable. 


The signal r(t), which coincides, or almost coincides, with s(t) + n(t), where n(t) is a certain new noise, is 
first separated from the received signal R(t) = S(t) + N(t) by means of phase detector PD (Fig. 4). This signal is 
supplied to the optimal filter OF, which is of the same type that has been described above. The result of filter 
discrimination s(t — T,), which is close to s(t—Tg), is obtained at the output. The indicated filtered-out signal 
can be used for producing the reference voltage necessary for the phase detector's operation. This voltage is prod- 
uced by multiplying R(t—T ,) by s9(t—T») and by transmitting the product through a selective filter, SF (for instance, 
an oscillating circuit), which is tuned to the frequency wp. 


Thus, the optimal filtering system that was calculated on the basis of the above theory represents a portion of 
a more complex system. The same situation can also be found in other cases. 


APPENDIX 
Generalization for the Case of Correlated Pulses and Correlated Noise 
If the correlations existing between the durations of intervals and pulses must be taken into account, it is 
necessary to consider the sequence 14, 04, 02, .. . a8 a complex Markov chain and to use, for instance, Eq, (9) 
instead of Eq. (10), The theory presented here can be extended to this case, but the corresponding expressions be- 
come more cumbersome. Thus, the function v = v(r, 0) will now depend on two arguments, and it will satisfy 


the equation 
9(t,0)= \{ 





wy, (t,o | t’,0") exp {~ P| (r + a)? dt— ui (r— a)tat} x 


° ? 


x v(t’, 0’) dr’de’, 


which can be considered as combining both of the equations (14). 
The probabilities w;,. (T) are given by the following equations: 


SK 
wi (T) = K ( ( exp {sy \ (r— a)rat} v (Ty, ,)dt,d3,, 
ty<T. 1%,>T T 


i 
wi =K | | wp ttetteas ce) exp f—gy | ot arrarl >» 
oy i<T. tH>T “om 
Xv (Ty 4s Gp_,) dt, ds _, dt,. 
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While the expressions 


const w,,, (t|o) = Wy, (c| t) exp |— sy t 


x \ (r+ ayrat} 
const w,, (s|t) = w,, (t19) exp {—; nT \! (r— a)rat} 


2: 


before played the role ofa posteriori (un-normalized) transition probabilities, the following expressions will now rep- 
resent such transition probabilities: 


const w,,, (t,0 |t’, 6’) = Wy, (t, 6| t’,0") exp \— sy \ (r + a)*dt — 


= Bana 


We shall touch upon the problem of extending this theory to the case of correlated noise. Generally speaking, 
the correlation inherent in noise results in the fact that the transition points do not form a Markov chain a posteriori. 
In this case, one must consider complex a posteriori chains, It is of interest to analyze the particular case where the 
noise correlation time considerably exceeds the sum Tcor «< Tpy + Tin Of the average pulse and interval durations 
(however this time can be comparable to Tpu OF Tip). In this case, the two-dimensional a priori Markov chain (9) 
remains the same a posteriori, and the a posteriori transition probability is given by 


const w,,, (t, s| t’,o') = Wy, (T, 6|t’, 0’) exp {— 


-w\ 


3’ 
a 
i 
~ 2N 
t 


Here, the kernel A(t, t’) is the inverse of the noise autocorrelation function 


aN\\! r(t)--alA(e’, [+ (t’)+-a]—dtat” 


[r(t) + a] A(t, ¢’) [r (t’) — a] dtdt’ — 


\ ir (t) — a) A (t,t) Ir (t’) — atat'| 


- 


ace 8 ate 


\4 (t, t’) mn (t') n (t") dt’ = 8 (t — 0”). 
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A STABILITY CRITERION BASED ON THE METHOD 
OF TWO HODOGRAPHS 


W. Oppelt 


(Darmstadt, Germany) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 9, 
pp. 1175-1178, September, 1961 

Original article submitted March 7, 1961 


A stability criterion for automatic control systems is considered, which is based on the separation 
of the hodograph (locus) of an open system into two separate hodographs. The criterion may be 
applied to continuous linear, nonlinear, and pulse control systems. 


To test automatic control systems for stability, methods are used based on the analysis of the configurations of 
hodographs (of complex frequency characteristics) of closed or open systems. For open systems one examines the 
configuration of the hodographs of Fy = Xoy1/Xjn» Where Xpyz and Xjn denote the output and input magnitudes of the 
Open system upon harmonic excitation of the input. Following Nyquist [1], one examines the position of the hodo- 
graph of F, relative to the “critical point” P, = 1 (Fig. 1). In what follows we limit ourselves to the case of smooth 
hodographs, which occurs most frequently in practice. In this case the closed system is stable, if the hodograph of 
F, of the open system intersects the positive real axis in the interval between 0 and 1, To insure a sufficient amount 
of stability, it is required in most cases that the intersection point be between 0.1 to 0.4. If the hodograph Fy» passes 
through the point 1, then there arises an undamped vibration with frequency w, correspondingto the given point of the 
hodograph. In this case the system will be unstable at the boundary of the stability region. 


The system is unstable if the hodograph intersects the real axis to the right of the point 1. 


If the hodograph Fp, does not go as it should, then the elements of the control system should be changed so as to 
obtain a better positioning for the hodograph. The influence of the individual elements on the character of the curve 
is difficult to track down in general. This is easier to accomplish if one joins all constant elements of the control 
circuit into a link with frequency characteristic Fs and all the variable elements into a link with characteristic Fp. 


Only a moderate number of elements of the system should be included in the link with characteristic Fp, since 
a simple hodograph should be formed for Fp for which it is easy to establish the connection between its configuration 
and the parameters of the variable elements. 


The frequency characteristic of the open system Fo, on which the stability criterion of Nyquist is based, is 
determined by the expression Fy = -FpFs. Therefore the stability of the closed system may be determined on the 
basis of the consideration of the hodographs of Fp and Fs. The hodograph of Fp may be obtained as the quotient 
Xout /X (Fig. 2). 


For the input magnitude x;, we then find xjp/x = —1/Fs. The hodographs of Fp and —1/Fgs may be taken as 
measures for the input and output magnitudes respectively of the open system if one considers points on the two 
curves which correspond to the same frequency. 


The passage of the Nyquist curve through the positive real axis corresponds to the intersection of the hodographs 
Fr and —1/Fs at points with identical frequency, with a straight line issuing from the origin (dotted line in Fig. 2). 
In this case if the value of F, is less than the value —1/Fs, then the open system is stable. To obtain a sufficient 
amount of stability the ratio of the magnitudes | Fp | /| -1/Fs | should, as before, lie between 0.1 and 0.4. 


The criterion of Leonhard [2] (for linear systems) was given in this form. The limits of its applicability were 
investigated by Cremeter and Kolberg [3]. It is particularly significant that this method of two hodographs Fp and 
—1/Fg may serve as basis for the investigation of nonlinear and pulse control systems. 
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Fig. 1. Open linear control system and the hodograph of 
its frequency characteristic Fo. 
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Fig. 2. Separation of a linear control system into mildly 
variable links (Fp), rigidly prescribed links (Fs) and 
elements changing sign. The corresponding separation of 
the hodograph of Fy into two curves: Fp and —1/Fs. 


Nonlinear Control Systems 








We concentrate all linear elements into the link Fs, which corresponds to the hodograph —1/Fs. The link Fp 
is nonlinear; the hodograph of the equivalent impedance of the nonlinear element serves as hodograph of Fp in this 
case. The equivalent impedance is obtained from the first harmonic in the Fourier series expansion of the nonlinear 
characteristic of the variable output magnitude xo, when the input magnitude x varies harmonically. Thus this 
method is approximate. 


The equivalent impedance in the general case depends on the amplitude x, and the frequency w; for simple 
nonlinearities it depends only on X». Fig. 3 shows the case of ideally nonlinear link with friction, The intersection 
points of both hodographs correspond to damped vibrations in the open system. The frequency w and amplitude xX» 
of these undamped vibrations may be reckoned with reference to the marks on the hodographs. The stability of the 
vibration may be determined if one considers the motion of the system near the point of intersection of the hodo- 
graphs, The point A corresponds to stable damped vibrations. Undamped vibrations may not be sustained at the 
point B, since a slight change in amplitude leads either to a stable vibration at a higher frequency (if the amplitude 
is increased), or to damping of the vibration (if it is decreased). 


If the equivalent impedance associated with a complicated nonlinearity depends not only on the amplitude x», 
but also on the angular frequency w, then our considerations are not changed. 


The method of investigation of the stability of nonlinear control systems with the aid of equivalent impedance 
has been developed in many papers; above all in those of L. S, Gol'dfarb [4] and Kochenburger [5]. 


Pulse Control Systems 
The method of two hodographs may also be applied in the construction of stability criteria for pulse control 
systems, if the pulse element is a fixed element. In Fig. 4 such a system is represented, The linear elements are also 





1049 








grouped into a link which corresponds to the hodograph —1/Fs. The characteristics of the union of a pulse element 
with a delay element may be represented with the help of an equivalent impedance [6], Upon harmonic excitation 
at the input of these elements there appears a step curve at the output, from the Fourier expansion of which one 
separates out the first harmonic. Then the hodograph of Fp (Fig. 4) is constructed. It underlies all pulse control 
systems, but the curve ~—1/Fs changes from case to case. 
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Fig. 3. Separation of a nonlinear system into linear links (Fs) and 
nonlinear links (Fp). 
































Fig. 4. Separation of a linear pulse system control system into a 
pulse element, a fixed element, and linear links (Fs). 


Since a control system containing a pulse element and a fixed element is linear, Fp depends not only on the 
amplitude, but also on the ratio T/T, where T = 2n/w is the period of vibration, and T, is the interval between 
pulses, The intersection point of both hodographs, as before, again corresponds to undamped vibrations, The fre- 
quency we of the undamped vibrations may be reckoned with reference to the marks on the curve —1/Fp. By the 
marks of T/T, on the hodograph Fp one may determine this ratio, and thus determine the critical interval be- 
tween pulses, In the same manner one determines conditions on the boundary of the stable region, since in linear 
systems the presence of undamped vibrations indicates the boundary of vibratory instability. Depending on the mutual 
position of the curves, a decrease or an increase in the interval between pulses will lead to a stable state. Figure 4 
depicts the rare case when an instability arises for small T,. 


For moderate magnitudes of T/T, the method becomes unreliable since then the first harmonic in the Fourier 
series yields a very rough approximation. However for ratios greater than T/T, = 3, one nevertheless obtains ad- 
missible relations (see [7]). The method may be formulated directly for hodographs of Fo; i.e., for the Nyquist 
stability criterion; in this, Fp is constructed according to Fp and Fs. However in this case the method loses clarity. 
It was given in this form in the papers of Brown and Murphy [8]. 
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ON THE STATISTICAL AUTONOMY OF DYNAMIC PROCESSES 
IN OPTIMALIZATION OBJECTS CONTAINING CONTROL 
SYSTEMS 


R. I. Stakhovskii 


(Moscow) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 9, 
pp. 1179-1186, September, 1961 

Original article submitted February 22, 1961 


It is demonstrated that for a broad class of n-dimensional optimalization objects which con- 
tain control systems and are controlled by the adjustment of controllers there exist conditions 
for which the systems that perform optimalization with respect to each individual input are 
statistically autonomous, These conditions can be found by automatic search for the ex- 
tremum of the average value of the performance criterion if we study the magnitude of the 
correlation shifts as independent variables. 


Introduction 

Many optimalization objects are technological objects in which the most important parameters Xj (i = 1,...,n) 
(temperature, pressure, etc.) are stabilized by means of control systems, This permits the average values of the 
stabilized parameters to be maintained equal to the specified values irrespective of the appearance of random noise, 
Usually the specified values Xjg are determined on the basis of engineering considerations in such a way that a cer- 
tain performance criterion (the efficiency, the productivity, the cost, etc.) reaches an extremal value. However, 
since with the passage of time the optimalization-object characteristics vary in unforseen fashion, it follows that 
the specified mode of operation for the object ceases to be optimal. The problem arises of performing a continuous 
automatic selection of Xj9 in such a way that the performance criterion Q has an optimal value [1]. 





The realization of optimalization requires the measurement of the quantities which characterize the interrela- 
tionship between Xj9 and Q. Usually in order to determine the characteristics of this interrelationship (i.e., the quan- 
tities dQ/dXj9) the optimalization object is subjected to small trial perturbations at the inputs Xj9, and the response 
of the object output (i.e., the variation of Q[1]) is studied. Such a method is possible but not alwaysconvenient be- 
cause of engineering considerations and because of the effects produced by the random noise which gives rise to 
random fluctuations of Xj and Q and thus masks the effects of the trial perturbations. 


Therefore it is of importance to seek a method for determining the interrelationship between the random 
processes without requiring the application of any trial input. Here we make use of the statistical characteristics 
of the processes expressed in the form of mutual correlation functions [2] in order to determine the characteristics 
of the object which is subjected to random noise. This method permits us to investigate the interrelationships with- 
out disrupting the normal operation of the object. A knowledge of the object characteristics makes it possible to 
control them using the natural fluctuations of the system [3] as the trial movements. 


However, in the case of multiple-loop systems the result of the correlation will contain more than just useful 
information on the interrelationship between the processes. Besides the useful component, it will contain an un- 
known component which depends on the effect of extraneous random factors if the latter are correlated with the 
processes with respect to which the mutual correlation function is determined. For example, it would be possible 
to study the fluctuations of the stabilized parameters Xj (i.e., the dynamic errors x;(t) of the controllers) by treating 
them as random trial steps which cause random fluctuations q(t) of the performance criterion [4]. Then the mutual 
correlation characteristics for the quantities x;(t) and the fluctuations of the performance criterion 4(t)—Rxjq(T) 
could be used to determine a quantity analogous to 8Q/0Xj9 which can easily be used as the control signal during 
the search process [1]. However, in the case of multiple-loop systems x;(t) are correlated with each other. More- 
over, the causes for the appearance of q(t) in the general case are not just the quantities x;(t) but the noise signals 
proper which are correlated with x;(t). 
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In this regard it is useful to study the conditions for which the quantity Ry ig (T ) would be proportional to 
8Q, /2Xjo (Qo is the average value of the performance criterion); i.e., we wish to study the conditions for which 
this quantity would change sign for the transition of each individual function Qo(Xj9) through the extremum, Under 
these conditions other input signals are treatedas aparticular case of noise. This permits us to “decouple” the op- 
timalization system so that we can treat individual inputs. Here the indicated conditions are the autonomy con- 
ditions (in the statistical sense) for systems which perform optimalization with respect to the individual inputs. 


Thus we are treating a particular case of the optimaliza- 
Ait) p(t) p(t) tion problem; the special feature of this problem is the pres- 
| passeee ive Teg | TR SGT Hae | ence of noise in an object which is controlled by means of 


Optimalization object | Qf) setting the controllers. 
4,(t) n(t) 
% }- Figure 1 shows an optimalization object which is con- 
ae a MnP) trolled by n control systems, This means that the object con- 
z(t) z(Y tains n technological parameters, each of which is stabilized 
Ang git) near a specified value Xj9. Each control system consists of 
a pick-up device for a stabilized parameter, which produces a 
Fig. 1 quantity representing the stabilized parameter X;(t) at its out- 
put. After comparing X;(t) with the standard value the system 
produces a quantity x;(t) which is amplified byablock with 
the transfer function Wj(p) and is applied to an integrator (an actuator) with the transfer function kj/p. The integrator 
output is a certain point in the object at which the control signal is applied. The section of the control loop which 
closes the system consists of that sector of the object being optimalized whose boundaries are defined by the points 
at which the control input is applied and the stabilized parameter is measured. The transfer function must be such 
that the control system is stable. All n control systems are mutually coupled; i.e., 
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The optimalization object is subjected to m random noise signals Pi(t). The noise signals are assumed to be 
stationary random processes. We introduce a certain performance criterion Q which in general case is also a random 
process. The magnitude of the criterion can be measured directly in the object, or it can be determined continuously 

by a high speed computer on the bases of measurements of the instantaneous 
Ae Sith values of definite output coordinates of the objects. In particular, these 
oat e-r- may be the stabilized parameters X;(t). If a computer is used, then it is 
ad ) Object | : a(t included in the concept of the optimalization object. As far as the aver- 
: caer i sity age value of the performance criterion Qy (its mathematical expectation) 
is concerned, we make the assumption that it is a nonlinear function of 
Fig. 2 the adjustments Xjo, . . . , Xpo Of the controllers which has one extremum, 
, The value of Qy is extracted from Q(t) by means of an averaging filter Wr. 
The fluctuations q(t) = Q(t)—Q, of the performance criterion are obtained 
by summing Q, and the inverted value of Q(t). 


Due to the effect of the random noise signals P(t) which are applied to the object, the object is in a state for 
which continuous random fluctuations occur about a certain average mode which is defined by the adjustments Xj, 
of the controllers. Here P;(t) can be treated as input perturbations, and X;(t) (or the dynamic errors xj(t) and Q(t)) 
can be treated as output random processes. Thus we can write 














X4 (t) —_ Xin + XY (t), 
Q() = @+q), (1) 
P; (t) = Po + 7p; (d), 


where Xj is the value for the setting of the ith controller, Q, is the mathematical expectation of the performance 
criterion, Pjp is the mathematical expectation of the jth random noise signal, x;(t) is the dynamic error of the ith 
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controller which is a centered random process, q(t) are the fluctuations of the performance criterion relative to the 
median level Qo, and p;(t) are the fluctuations of the jth random noise signal relative to Pjg. 


The problem consists of finding the extremum point for the function Qy(X4o, . . . » Xo) using only the values 
xj(t) and q(t). 


We shall assume that Pj(t), x4(t), (t) are small, and that Xj9 (before search is begun) are independent of time. 
These propositions make it possible to study an essentially nonlinear optimalization object as an object which is 
"linear in the small” if the set Xj9 does not correspond to the extremum point of the function Q)(Xqp, . . . . Xno)- 

In that case, in accordance with the theory of linear transformations of random functions [5], it is possible to perform 
a separate study of the passage of mathematical expectations and centered components through the system. Since 
according to the definition accepted for the time being mathematical expectations are independent of time, it follows 
that the nonlinearity of the object "in the large" (for mathematical expectations) is of no significance in deter - 
mining the mutual correlation functions for the centered components. “ 


One of the purposes of the present analysis is the realization of a method for determining a signal proportional 
to the partial derivative 0.Q, / 2X49. We shall study the problem of determining the interrelationships between x;(t) 
and q(t) as the problem of determining certain components of the correlation matrix for the vector random function 
X Q(t) whose components are xj;(t) and q(t). Here X Q(t) is the result of transforming the vector random function 'P (t) 
whose components are pj(t) (Fig. 2). This problem was solved on the assumption that the object is linear and is 
characterized by the matrix of "transfer" operators Aj; which consists of m(n + 1) terms [5]. A "transfer" operator 
Ajj is defined as an operator which represents the effect of the function p;(t), which is treated as an input perturba - 
tion, on the ith component of the function Xo(t) (i= 1,...,n+1). 


The difference between this teatment and the analysis performed in [5] resides in the fact that in addition to 
the matrix of “transfer” operators Ajj the object has a matrix of the operators Aj7 (i, 7 = 1,...,n+1). The 
operators Aj, are defined as operators which represent the effect of the input i on the output 2, The operators Aj; 
can be subdivided into two groups, First, there exist n “useful” operators Ajq which depend on the unknown quan- 
tities 9Q, /8Xj9. Second, there exist operators which reflect the interrelationship between the dynamic errors of 
the control systems. Moreover, if the noise signals pj(t) are correlated, then there also exist the operators Aj, (j, r= 
=1,...,m) (6}). 


In fact, if there are two correlated processes, for example p;(t) and p,(t), then in each of them it is possible to 
isolate the part which is not correlated with the other process. Here the correlated part is obtained using a certain 
equivalent operator which connects the measurement points p;(t) and p,(t). This means that without violating the 
general nature of our analysis we can study the case of statistically independent noise signals Aj, if we take into 
account the presence of the operators Ajr. 


All of the concepts cited above do not make it possible for us to use the formulas derived in [5]. 


On the Existence of the Autonomy Conditions 
We shall study the case when the only input is the noise pj(t) and we shall be concerned with the properties of 
the mutual correlation function Rxjq(T). Here (making use of the linearization method [5]) we shall assume that 
the system is linear for small values of the fluctuations p;(t), x;(t), q(t), The assumption that these quantities are 
small is justified by the very existence of control systems which reduce the 
fluctuations of the stabilized parameters to a minimum. 





Ad 
| We shall transform the network in Fig. 2 (treating only the jth noise 
signal alone) with respect to the ith “useful” operator (Fig. 3). It should be 
noted that the fluctuations q(t) consist of two components: one of them is 








Ay hv caused directly by the dynamic error x;(t) which is transformed by the oper- 
0 ator Ajq, and the other is caused by the noise p;(t) transformed by the gener- 
lized operator Ajq. 
A gi) * iq 
Lyyp)(t) q 











Thus we can write 


Fig. 3 
3 a5 (t) = (Argus (0) + 25 (0), 2) 
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z;(t) is that component of the oscillations of the performance criterion which is caused by the noise p;(t) which 
al transformed by the operator Ajq; x4(jy(t) is that component of the dynamic error of the ith controller which is produced 
by the jth noise signal; qj(t) is that component of the performance criterion fluctuations which is caused by the jth 
noise signal, 


The transformation shown in Fig. 3 can always be performed for any form of the linear operators Aj, Ajj, Ajz - 
The operator Ajq, which is a particular case of the operator Aq; (J = q), can be written as 


C 
Aiq = see > Ala 


where Aig is the operator for a “unit” section which has the same dynamic properties as the sector i—q and is linear 
for small values of the dynamic errors x;(t), In other words, Ajq has a static gain equal to unity. Then expression 
(2) can be rewritten as 


ie 


qj (t) = * {Aig Lj) (t)}+ 2; (2). (3) 


We shall multiply both sides of expression (3) by xj(j)(t + T) and shall find the mathematical expectation of 
both parts 


M (x45) (¢ + 1%) 95 (t)) = Raja; (2) = = tM [205 (t+ t) (Aigrus (O}) + 


+ M [xyy (t + 2) 201 = 20 Te Ria, (1) + Riga; (9s 


where a ie (x) is the useful component of the function Raggy (7) and is proportional to 0 Qo / AX ip; 
) (t) is the interfering component of the function R,. ia”? which must be determined by setting Ajq = 0. 


Before we proceed to a further analysis we shall determine the properties of the generalized operators A ji 
We shall juxtapose each operator with a pulse weighting function Bij(- It can be assumed that they characterize 
the segments of the control system between the point at which the arbitrary noise is 
applied and the point at which the stabilized parameter is measured, If the control 
—[ooi0 LS system is astatic and stable, then fij(0) = 0 and jj) = 0, where fj; (t) is the response 


























(t) of the above-mentioned segment to a unit step perturbation. If the control system is 
a static and its open-loop system has a high gain, then hjj (@) » 0. Then we can write 
2,(t} co 
(t) | ~ 
& hy (00) = ( gy (tat = 0. (4) 
Fig. 4 : 


This means that the graphs of the pulse weighting functions for all of the above-mentioned segments of the 
object intersect the time axis at least once. 


We shall now study the mutual correlation function between the results of transforming a certain random 
process y(t) using two sections (Fig. 4), one of which has a weighting function g(t) with the property (4). In that 
case we can write 


Z_(t) = \ y(t —9) gp (0) a8. (5) 


or 8 


Multiplying both sides of expression (5) by z(t + +) and finding the mathematical expectations, we obtain 


M [xy (t + 1) 2 (0)) = Rea,(t) = M La (¢ + ) Vy (t — 9) ge @) a0). 
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Changing the sequence of the operations of integration and mathematical expectation, we find 


620) dOM [z, (t + t) y(t —9)] = | Ray (x — 8) gq (0) dd. 


0 


Riz, (t) ms 


ou 3Q 


We shall now take into account the fact that g,(t) has the property (4) and shall find the magnitude of the 
integral of the function Rz 42,7): 


( R,,,,(t) dt = \ de\ Rey (x — 0) ga (6) a0. (6) 


We shall transform expression (6); for this purpose we alter the sequence of integration 


y R,,2,(t) dt = ( 6s (0) dé ( Rew (x — 0) dr. (7) 
—oo } bie 


The inner integral in the right side of expression (7) is equal to a constant, since it represents an area bounded 
by the graph for the mutual correlation function Rz y(t) shifted along the r axis by an amount equal to the variable 
parameter 9, The magnitude of the area is equal to a finite quantity, since Rz,y(—) = Rz,y(@) = 0 and Rz,y(r) is 
a bounded function. In this regard the inner integral in the right side of expression (7) can be treated as a constant 
multiplier c, Thus, 


| Reset) dt = cl g.(0) dd = 0. (8) 


Assuming now that y(t) = pj(t), z(t) = z;(t) and z(t) = x;¢j)(t), we find that Rsqy’T) has the property (8). This 
means that there exists at least one finite value 1 jo for which Rx, ay Tj9)= 0. The resulting derivation is valid for 


any noise P;(t). Making the transition to the case where m noise signals are involved, we note that both x;(t) and 
a(t) can be written in the following form in accordance with the superposition principle: 





x; (t) = D} zu (t), (9) 
j=1 
a(t) = 3 9;(t) => la { Ajgtigyy (t)) + 2; (t)| . (10) 
j=1 j=1 io 


Thus in determining the function Ral) {for this purpose it is necessary to multiply the right sides of ex- 


pressions (9) and (10) and then average the results] we find that 2m* components appear; each of these will represent 
the result of correlation between any pair xij) and q,(t). It is easy to demonstrate that of the 2m* components of 
the function qT) only 2m will be nonzero. In fact, the results of the transformations of statistically independent 
processes are statistically independent. Therefore we can consider just the components (2Q)/3 XR xiqft) 
and Rigi (r). Thus, - 


9Q, Pp . . o Ri K 
Rzq(t) = pa # Req; (t) + pa Ria = 7 Rea (t) + Rayq (7), (11) 
j= (0 j=1 io 


*79j 


Rey () = y Rija;(t)» Rega (*) =D Rea, (2)- 


J=1 
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We can demonstrate that the mutual correlation function Rx q(T) has the property (8). In fact, 


oo co m © 


\ Rigg (t) dt = \ D Fea, (ar => \ Rig, (t)de = 0. 


—oo —oo j= j =1 —oo 


This property of the function Ra (r) makes it possible for us to state that for each function Rasq(T) there 


exists at least one value Tj9 such that 


6 


0 pl 12 
Rs,q (Tio) = x Req (Tio): (2) 


This allows us to use Rx q(T io) as a measure of the magnitude of the sought for partial derivative. 


Determining the Conditions for Statistical Autonomy by the Method of Automatic 








Search 


We shall demonstrate that the quantities Tj», for which automatic systems that achieve optimalization with 
respect to separate inputs are autonomous in the statistical sense, can be found by means of automatic search for the 
minimum magnitude of the steady-state value of the performance criterion Qpogr s+, treated as a function of the cor- 
relation shifts r;. Under these conditions the entire automatic optimalization system is of the form shown in Fig. 5. 
The entire system consists of n separate systems which per~- 
form optimalization with respect to each separate input of 
the system shown in Fig, 1; here the input is defined as the 
iY) Object setting of the corresponding controller. Each individual 
system consists of a block that achieves controlled delay, a 
multiplier, a low-pass filter, and an integrator [7]. The 
average value Q, is measured at the output of the low-pass 
averaging filter. 


Ko 





If we close each separate automatic optimalization 
system (for certain arbitrary values of r;) through its ac- 
tuator which consists of an integrator with a time constant 
appreciably greater than the time constant of the low-pass 
averaging filter, then the stable state of each individual 
system (for the condition that negative feedback is used) 
is obtained if the input voltage of each integrator is equal 
Fig. 5 to zero on the average. This will correspond to the case 
where the left side of expression (11) is equal to zero, 
Therefore a condition for the stable state of each individual 
automatic optimalization system will be 





R t 
( i ) =— nn 2 (i.e Ie. 10 (13) 
io /Xio=Xio st.st, Ria 


The system of equations (13) will define a steady -state point in the space of the independent parameters X{o. 

This steady-state point will correspond to the specific value Qoggr st; i.e., it will correspond to the average value of 
the performance criterion in a steady-state mode. It is natural that when 1, varies the quantity Qpgr sr, will also 
vary, since the steady -state point in the space of the parameters Xjg varies. Thus it is possible to study the quantity 
Qost.st, 28 a function of the correlation shifts r; in the separate automatic optimalization system. If we assume for 
simplicity that the function Qogt st.(71, - - - » Tn) has one extremum point, then this extremum will coincide with 
the unknown extremum point for the function Qg(X1p, . . . . » Xo); here the values rj will coincide with rio. In 
fact, for rj = Tyo, Ry T to) = 0; therefore, (8Q,/0 Xioxso = xjgst.st > 0. In order to find the extremal point we use 


a single-channel automatic optimalizer AO. The case where the function Qogr s¢.(T1, +++» Tm) has several ex- 
tremums can be reduced to the case of a single extremum by an appropriate choice of the initial correlation shift. 
The problem of choosing the initial shifts is not studied in this paper. 
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SUMMARY 
1, For a wide class ofn-dimensional objects which are optimalized by changing the settings of controllers and 
which are subjected to random noise signals, there exist values of the correlation shifts in the mutual correlation 
functions for the dynamic errors of the controllers and the fluctuations of the performance criterion for which the 
systems performing automatic optimalization with respect to separate inputs are statistically autonomous, Thus an 
m-dimensional optimalization system reduces to m unidimensional systems. 


2. These values of the correlation shifts can be found by means of automatic search for the extremum of the 
average value of the performance criterion if we treat the magnitudes of the correlation shifts as independent vari- 
ables. 


3. When the search method cited above is used (this method does not require the application of trial inputs to 
the objects) the extremum for the function Qo(Xio, . . . , Xng) is reached automatically. 
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ESTIMATING THE EFFECT OF LEVEL QUANTIZATION 
ON PROCESSES IN DIGITAL AUTOMATIC SYSTEMS 
WHEN A RANDOM INPUT SIGNAL IS USED 


V. A. Volkonskii 


(Moscow) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 9, 
pp. 1187-1193, September, 1961 

Original article submitted February 18, 1961 


The paper determines the probability characteristics for the level quantization noise in digital 
automatic systems when a Gaussian random process is applied to the system input and the quan- 
tization step is small compared to its mean square deviation, 


In digital automatic systems (DAS) quantities are subjected to quantization with respect to level and with res- 
pect to time. Level quantization consists of rounding off the quantity (i.e., it consists of substituting the closest 
value which is a multiple of the “quantization step” A for it). In [1] we provide an estimate of the “quantization 
noise” — the maximum deviation of the output signal of a DAS from the output signal of the same system under 
conditions where no level quantization is performed (the limiting system for 4-> 0), Although this estimate cannot 
be improved for certain "worst" input functions, it is nevertheless true that in the majority of cases it proves to be 
much too high. Therefore it is logical to examine how the probability characteristics of the quantization noise vary 
in the case where a random process is applied to the input and the quantization step A tends toward zero, 


We shall study the DAS systems shown in Figs. 1 and 2 where K and L denote linear elements, and AD denotes 
a unit which performsa level quantization(conversion of an analog quantity to a digital quantity), Here we shall 
study only the case where time quantization is not performed or when its effect can be neglected. 
























































(t 
a(t) ‘ap + Fi Gt" 
x(t) y(t) 
AD ~~ K _ 4 | AD owe 
a 
Fig. 1 Fig. 2 


We shall use x(t) and y(t) to denote the input and output signals. We shall assume that x(t) is a Gaussian 
stationary random process with zero mathematical expectation, a disperion 0%, and a normalized correlation function 
t,- We shall find the approximate form of the correlation function for the process: 


ds) = ya) — yo (O, 
where yo (t) is the output signal from the limiting system (the system without quantization), 
In addition we shall introduce the following substitutions * : 


* x"(t) and y'(t) denote the first derivatives of the processes x(t) and yo(t); 19" denotes the second derivative of the 
correlation function r, at the point t = 0. 
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6, ° , ” ‘ ’ 5, 
o% = Dy (t), B= =, 62 == Dz (t) = — rio, 0? = Dy, (t), k= = l=. 


The quantities fx and fy are the mean square frequencies of the processes x(t) and yo(t). If gy (w) is the 
spectral density of the process x(t), then 


_V/ 6 2%) 
nat ta = do. 


—oo « 








Here X, shall denote the result of applying level quantization to the signal x; i.e., it is the quantity closest to x 
which is a multiple of the quantization step: 
~ ° 1 
Za = Es + a A, 
where [x] denotes the integer part of the number x and y , x denotes the difference 


YVat=2z —Z,. 


(1) 
1. The Case of an Open-Loop System 





Assume that for large 8 the distribution of the quantity *4* (!) 
segments [~1/2, 1/2]. Therefore for B-> , 





is close to a uniform distribution over the 


Fol 


1 
ar P¥az (t)— 
and 


z?dz= 


TF Ld 
—_— 
3 = 


ar Dvaz(t) 1. 


We shall use the substitution F(t) =‘M [paz (0), pax (t)]. In the Appendix we prove [cf. Eq. (13a)] 
that 


R(t) = (2) Xt exp [— (eng) 


m? + n2 


= ]su [(2xB)*r,mn]. (2) 


We shall now study the normalized correlation function T(t) = ar Ri). Assume that | r, | < 1 for t = 0. 
Then from formula (2) we find that for 8-> « and t # 0, we have 





F(t) ~ 5 exp [— (28) (1 —7,)]" 


(3) 
this quantity is already negligibly small far 8>>1/ —2— 


— However, the function (sm) (the correlation 
wad 

function of the process 2, (t) = “eS x az (sr))) coverges for 6 -» ©, as is indicated by formula (2), to 
the function 


foo} 
—= >) ax exp (rt*k?), 
x i” 0 


*a-~b denotes that (a/b) +1. 








where ro" = f*, is the mean square frequency in the spectrum of the process x(t) (the terms m =n with in formula 
(2) drop out), 


Here the spectral density of the process ¢, (A) will converge to the function 


6 V 2x 
9 (A) = VEY te oxp(— ne) 


Assume the quantization noise y ,x(t) is applied to the input of a linear element that converts the signal x(t) 
to the signal 


Kz (t)= \ k(s) x(t —s)ds, 


—oo 


and assume that the transfer function for this system is k* (A). Then the spectral density of the process Ky Ax (t) will 
be close to the function 


sax | (A) (sa) - 


For 8B + « 





® (sar) > 90) = % Ae 


ty 


Therefore the spectral density of the process Ky ,x(t) can be assumed equal to 





Bea cet ie a) (#) 
where 
Cc 2Vz Gr 0.8. (5) 


The error of this approximation can be estimated, for example, as follows. If a linear system does not pass 
harmonics with frequencies exceeding f ymax (i.e., if K (A) = 0 for | A] >fmmax)» then the relative error does not 
exceed the value 





(4 ) 
?\ B / V2 on 2 0. , | 
1 ay Spe (Sis ) (4% 3) ~ =(—) ” 
Thus, for example, if 
fx max< Bfz, (7) 


then the use of formula (4) leads to an error that does not exceed 1%, 
Correspondingly, the correlation function p; of the process Ky,x(t) can be assumed equal to* 


* Moreover, it can be proved that the normalized process (V/A) K4=(¢) converges to a certain stationary process 
for A - 0 (or B -» «); thus it is possible to find the probability distribution for it. We do not cite this proof here. 
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k(s) k(t — s)ds. 


jg 


Note also that even for 68 = 1 the quantization noise » , x(t) can be assumed to be uncorrelated with the input 
process x(t). In fact, in the Appendix we prove [cf. (14a)] that the mutual correlation function for the processes 


 ax(t) and x(t) is R(t) = Q (8) o%r,, where @Q (3) —=> c e—2(=8k)* and that the quantity Q(B) is 
already negligibly small for 6 = 1. ' 


2. The Case of a Closed-Loop System 
Assume that for passage through a linear element the function x(t) is subjected to the operator L with the transfer 
function 1*(w), Then the equation for determining y(t) is written as 





ys (t) = L [2 (t) — paz (t) — Ya (t) + paya (O), 
hence we find that 


ya (t) = K [x (t) — faz (t) + Paya (I), (8) 
L* (@) 

1 -+/* (@) 

it is an operator which must be applied to the input signal in order to obtain the output signal for the limiting system 


where K = (E + L)L is a linear operator with the transfer function k* (@) = (E is a unit operator); i.e., 


yo (t) = Ka (t). (9) 
Therefore 


5a (t) = yo (t) — yo (t) = — Kaz (t) +- Kwaya (0. (10) 


In the right hand side of equation (10) the operation Ky q is applied to the processes x(t) and y,(t). If both 
of these processes were to be Gaussian, then it would be possible to clarify the properties of the process 5,(t) by 
applying the conclusions of section 1, Although the process y ,(t) is not Gaussian, it is nevertheless close to the 
Gaussian process y¢(t) for small A. It can be proved that for natural assumptions 


Kpaya (t) = Kayo (4) (1 + aa (2), (11) 


where a a(t) is a stationary process for fixed A, and a ,(t) tends toward zero with respect to probability when A=,’ 
Based on relationship (11) we can replace Kipyy,(t) by Awayo (t). in (10). 


Now making use of the results of section 1, we find that the correlation function for the process can be assumed 
equal to 


c= \ k (s) k (t — 8) ds. 
Sy 2. 
Thus for a sufficiently small step A the error 5,(t) can be assumed equal to the result of passing the difference be- 


tween the quantization noise in the output signal of the limiting system and the input signal through an inertial 
system : 


53 (2) ~ Kayo (t) — Karz (0). (10a) 


* For example, this can be proved for the case where a linear system has a bounded frequency pass band; this con- 
dition is practically always satisfied. Under these conditions the quantity a ,(t) has a dispersion of the order of 
Vs. 
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Note that if the normalized mutual correlation function pr for the process x(t) and y(t) does not exceed a certain 
number p < 1 in absolute magnitude, then the normalized mutual correlation for their quantization noises is al- 


ready a negligibly small quantity for $ > V ~s . In fact, by analogy with relationship (3) it is easy to find 
that for 8 -> « the mutual correlation function ry == + Max (0) payo (2) satisfies the relationship 
= 6 . ¥ 
bt ~ =z exp {— (2m8)* (1 — px)} (—00 <t < 0) 
2 
i—p° 
can be assumed uncorrelated, and the correlation function for the process 5,(t) can be assumed equivalent to the 


sum of the correlation functions for the processes Kipayy (¢) and Arp,az (t). Applying the conclusions of sec- 
tion 1 to these processes, we find that the correlation function for the process 6 A (t) can be assumed equal to 





and is a negligibly small quantity for 8 > In that case the processes Kupayo (t) and Awaz (t) 


pic +=) Jk ke — s)ds. 





If we do not assume that 3 > then from relationship (10a) it is possible to find the estimate of 


= —p’ 
the dispersion for the process according to the formula 








Dd~4 (t) <(V VY [Kpayo (t)] + "= D [Kaz (t)|)*= 





k? (s) ds. 


=ca*( 4 vz) \ ( 


SUMMARY 
1. If the input of a digital automatic system is subjected to a Gaussian stationary process x(t) and the width 
fxmax Of the frequency band passed by the linear element does not exceed the mean frequency f of the input signal 
multiplied by 8, then the quantization noise y Ax(t) can be treated as white noise; i.e., it can be treated as a process 
with a constant spectral density 0.3(A°/a x) which is not correlated with the signal x(t). 


Here the correlation process for the function 6, (¢) = Awaz (t) can be assumed equal to 


Pt = 0,3 \ k(s)k(t—s)ds. 


e 
« —20 


2. In the case of a closed-loop system and identical conditions it can be assumed that the dispersion of the 
process 5 (t) = ya (t) — yo (t) is 


Dia (t) < 0,34" ( vr vz): \ k*(s) ds. 





If under these conditions the normalized mutual correlation function for the input signal x(t) and the output 
signal y(t) = Kx(t) of the limiting system does not exceed p <1 and the condition 8 > V = is valid, 
then the correlation function for the process 6 , (t) can be assumed equal to 


ro) 
1 


0,848 ( 7 + x) | k(s)k(t—s)ds. 


S 
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APPENDIX 
The Correlation Function for the Level Quantization Noise in the Case of a Gaussian 


Process, and the Mutual Correlation Function for the Process Plus the Noise 
Assume (€, 0) is a pair of Gaussian random functions with zero mathematical expectation and a second- 





moment matrix 


a, ro 
ro*, oi} 


and R = M (yt, bani (ef. (1)}. We shall introduce the 














We shall find the quantities # = M [E, », 1/ 
substitution B = o/ A, 

Assume p(x, y) are the distribution densities for the vectors (E, 1). Then the distribution densities for the 
vectors .~,&%, 1 and (&,%,%) will respectively be equal to 


= A 
P(z,y) = yp et id,yt kd), leh lylSz’ 
jk 


‘i A 
P (z, y) = Dip ytkd), lulz: 
k 


Assume now that (A, w) is the characteristic function of the vector (€, 1); i.e., 


2 
@ (h, @) = M [exp (i AE + on) = exp {- F OF + 2rhw + w} ; (12) 
Then 
a 2n 
p(z,y) = > Qmnn &XP \-i 9 (mz + ny} ’ 
where 
oie sam i (2x 2x \ 
NE ct 
and 
ry o* ky 
Bee, wW =>) \ age Pre 4, 
k —co 
where 
a, (A) sax? (* 2 ke) 
From this we obtain 
a4 4 
2 2 
= = 2 (—1)mt"t! on On 
R= \ \ zyp(z, y) dzdy= (=) > BO (WR (= m, n) : (13) 
4 A m.n+0 
ee” PE 
A 
ee A < (—tk @ 
ae - pm US A 2x 
it \ \ zyp (*, y) dydz-=— 2x » k OX sf (0, wy k) ° (14) 
-o A k+0 
2 
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Making use of the expression (12) for the function g(A,w) and performing simplifications, we obtain® 





= rod __4ym+n 
R an (2) > ay exp {- (2B)* =t “ sh (2B) *rmn]), (13a) 
™m n=! 
R = Q (6) o*, one 


where 


Q (B) mt on >> (—1)* etn Bay | 


k=1 


In particular, if € = x(0), n = x(t), where x(t) is a stationary Gaussian process, it follows that r = r, will be the 
normalized correlation function for the process x(t). Under these conditions formulas (13a) and (14a) yield the fol- 
lowing expressions for R = R(t) (the unnormalized correlation function for the process » ,x(t)), and R = R(t) (the 
mutual correlation function for the processes x(t) and » Ax (t) ). 


LITERATURE CITED 
1, Ya. Z, Tsypkin, “Estimating the effect of level quantization on the processes in digital automatic systems,” 
Avtomatika i Telemekhanika, 21, 3 (1960). 
2.  8B,R. Levin, The Theory of Random Processes and Its Application in Radio-Engineering [in Russian] (Soviet 
Radio Press, 1960), 





All abbreviations of periodicals in the above bibliography are letter-by-letter transliter- 
ations of the abbreviations as given in the original Russian journal. Some or all of this peri- 
odical literature may well be available in English translation. A complete list of the cover-to- 
cover English translations appears at the back of this issue. 





* In contrast to [2] the expression for the quantity r is derived without using Fourier transforms [cf. [2], formula 
(7.124)}, 
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SYNTHESIS METHOD FOR FINITE AUTOMATA 


V. G. Lazarev and E. I. Piil' 


(Moscow) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 9, 
pp. 1194-1201, September, 1961 

Original article submitted February 14, 1961 


A synthesis method is presented for discrete circuits with feedback, containing the operators 
AND, OR, NOR, d (transition operator), and feedback elements (relays, triggers, delay lines, 
etc.). In this connection selection criteria for various types of feedback are examined. 


In automation, remote control, communications, computer technology, multisequence circuits (finite automata) 
constructed from contactless elements find a wide application. To synthesize these circuits logical devices are wide- 
ly used [1-5]. 


It is possible to represent multisequence circuits by oriented (n, k) poles with s feedbacks, where n is the num- 
ber of inputs, k the number of outputs. In the case where combinations of input signals uniquely determine com- 
binations of output signals, there is no feedback. Then it is possible to reduce the construction of the circuit to the 
synthesis of an (n, k) pole with operators AND, OR, NOR, if the input and output signals are of the same type (poten- 
tials or impulses). If, however, the input signals are potentials and the output impulses then, besides operators AND, 
OR, NOR, the transition operator d is used [6]. 


In the case when combinations of output signals are not uniquely determined by combinations of input signals 
it is necessary to introduce additional signals at the input end to satisfy the given conditions, These additional sig- 
nals can be obtained with the help of feedback loops. Then, we introduce the necessary number s of these feedbacks 
and build the (n + s, k + 8) pole, where the outputs n+ 1,..., +8, are the inputs to the corresponding feedback 
elements (FE) and the additional inputs n + 1, ..., n +s are the outputs of FE.° 


The number of necessary feedback signals is determined by both the given operating conditions of the circuit 
and the inherent properties of the feedback schemes. 


We shall examine those FE which give rise to potential signals at the output. Then, depending on the nature 
of the inputs, the FE can be divided into two basic groups: 


1) FE reacting only to combinations of input signals, 
2) FE reacting to signals corresponding to some transition from one combination of inputs to another. 
FE in the first group are called FE of type A (FE-A); FE in the second group FE of type B(FE-B). 


Considering that FE introduce a delay r (rt < Tyjp,. where T,)j7 is the minimum duration of the combinations 
of input signals), we shall examine such FE-A the output signals of which repeat the input signals with delay 1, and 
such FE-B the output signals of which appear with delay r after the appearance of the switching-in signal at the in- 
put and vanish in time r after the appearance at the input of the switching -out signal. Electromechanical unpolarized 
relays and delay lines are examples of FE-A; thyratrons and triggers are examples of FE-B. 


Different types of FE require different approaches for switching them at proper moments into the circuit opera~ 
tion for the purpose of transforming unrealizable conditions to realizable ones. 


FE-A have been studied by many authors [1, 4, 7-11]. Estimates of the minimal necessary number of FE have 
been obtained [8, 9, 11] and conditions for realizing the switching tables have been formulated [10], 


* Input and output signals of FE correspond to one and the same ordinal number, therefore the additional output sig- 
nals have the numbersn+1,...,n+sandnotk+1,...,k+s. 


1066 





ons 


\= 


fa~ 


/e 





All these problems have not been sufficiently studied for applications to FE-B, although the latter are widely 
used in practice. 


The present paper formulates the conditions for switching FE-B into the circuit operation, gives an estimate 


for the minimum number of them for realizing the conditions, and proposes ways to synthesize finite automata using 
FE-B, 
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We assume that the operating conditions of the finite automata are given in the form of switching tables [7, 10]. 
Then to transform unrealizable switching tables to realizable ones it is necessary to switch in s FE. 


If we use FE reacting to a combination of input signals (FE-A), they can be switched-in independently of previ- 
ous switching -ins in sequence with such combinations of input signals which did not occur in the previous sequence. 


However, if we use FE reacting to transitions from one combination of input signals to another (FE~-B), they can 
be switched-in independently of previous switching-in signals, produced by the transitions, not occurring earlier. 


Since for a given number n + s of input signals the number of different transitions can be greater than the num- 
ber of combinations, then the number of sequences where it is possible to switch in FE-B is greater than the number 
of sequences where it is possible to switch in FE-A. 


For instance, for n + s input signals the number of combinations can be a maximum of 2"*5, however, the num- 
ber of transitions may be (2"*8—1)2"+8 if coincidence of transitions is allowed, i.e., if simultaneous changes of two 


ot more input signals are permitted; and, the number of transitions is (n + s)2"** when coincidence of transitions is 
not allowed. 


Therefore, in a number of cases the use of FE-B makes it possible to reduce the number of FE necessary for 
transforming unrealizable conditions to realizable ones. 


For instance, let the switching table* be given as shown in Fig. la. 


Then, if we use FE-A to transform the unrealizable switching table (Fig. 1a) to a realizable one it is necessary 
to use three FE-A (Fig. 1b); however, if we use FE-B then only two FE-B are necessary (Fig. 1c). 


However, we have not succeeded in changing the estimate previously obtained [11] for the necessary number 
of FE. In fact, if there exists m,,,, sequences with identical states of the input signals to every one of which for 
given conditions there must correspond different combinations of output signals (coincident sequences), then to dis- 
tinguish between these sequences it is necessary to introduce not less than 2°N states of the secondary input signals. 
Here sy is the lower estimate of the necessary number of secondary signals (or the number of FE), 


We have also not succeeded in decreasing the upper estimate, because in the case when m,= Mg =. . . = Mays 
when using FE-B we get for the upper estimate [11]: 


28u > 2mmax —1. 
* In the switching table we use the following notation convention: x; are the input signals, 2; are the output signals, 


= denotes the presence of a signal at the fundamental input, — denotes the presence of a signal at the secondary in- 
put or at the output. 
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Thus, without lowering the estimate of the number of FE, the use of FE-B makes it possible to reduce the num- 
ber of cases in which the upper estimate is realized. 


In those cases when the number of necessary feedback signals cannot be decreased, it is possible to reduce the 
number of states of these signals we need to use, i.e., a smaller number of their switching-ins and switching -outs. 


If we use FE-B we change to method of recording the signals entering FE. Since the output signals of the cir- 
cuit which are the inputs of the FE are produced by the transition of one combination of input signals to another, then 
it is necessary to describe these transitions, They are described by means of potential-impulse forms, introduced by 
A. D, Talantsev [6], which besides operators AND, OR, NOR, also switch in the transition operator d. Here we dis- 
tinguish between the transition from a unit state of the input signal to a zero state and vice versa: dx(t)= 1 is the 
transition from a unit state to a zero state, dXx(t)= 1 is the transition from a zero state to a unit state, dx(t) = 0 and 
dX(t) = 0 are not transitions.* 


Thus, for instance, the feedback signals needed for realizing the conditions given in the switching table (Fig. 2) 
is written in the following way (we denote the output signals entering FE by y,): 


Ys = Te Ta I4%,dz7, V Ty Ly 4 Ls dt», 
Ya = ry Ly % xz, dz, \/ ©, Ts T, aX, da *. 


Ys = Ly Lz Tq Tz Az2. 


The output signals of the fundamental outputs, having a potential character, are written in the usual form 
(3, 10): 


2y = My Bq Ly Tq Ly \f Ly Le Ly Ly Ly \/ Ly Ly Hy Ly Ts \/ Wy Le Ly Ly Up \f Ly 1p LyX Lg, 


Ze = Ty Lo X_%q Fy \/ Xy Le Ly Lg Ly \/ Uy Te Ts, Ly Ty \/ Ly Xp Ly Te Tz. 


The method of constructing circuits to obtain potential signals is well known [2, 12-14], It would be desirable 
to express the signals for switching in FE-B in the usual form to construct the entire circuit, in the first place to make 
use of already known methods, and in the second place to 
obtain a unified circuit for all k + s outputs of the (n +s, 

k + 8) pole. Such unification of circuits permits the reduc- 
tion to the smallest number of operators AND, OR, and NOR, 


Signals to switch in FE-B can be described by Boolean 
functions if we can find such a function F that y = dF, i.e., 
if there exists such a potential function F which after dif- 
ferentiation gives the necessary impulse signal y. Here the 
circuit takes the form shown in Fig. 3. Here every signal 
except y is a potential. 


“ 














The process of finding the function F is called integra- 
tion of the potential-impulse form y [6, 15]. 


Function F can be obtained immediately from the 
operating conditions of the circuit given in the tables, In 
fact, for instance, if in the transition from the fourth to the 
fifth sequence (Fig. 2) it is necessary to obtain the signal 
for switching in FE-B, then the function F; in the fifth sequence must necessarily have the value zero, and in the 
preceding fourth sequence must necessarily have the value one (if we use the transition from the unit to the zero 
state; analogously it is possible to use the inverse transition and to get the function F). 
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* For conciseness the parameter t will be omitted in what follows. 
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Thus, for our example function F; can be written as 
Fy =% Te Lg Xe Zs V Zy Lq Tz XM Zs: 


The state of the input signals at which the function F; must necessarily have unit and zero values is defined 
(Fy = 0 at XyXpXgXqXs and XyX~X~X4Xs). 


Thus we can obtain function F written in the same form as the fundamental output signal z, and therefore can 
construct the circuit by known methods from elements AND, OR, NOR [2, 12-14], The signals described by function 
F enter at the inputs of the operators d at whose outputs are obtained 
the signals necessary for switching in FE-B (Fig. 3). 


As an illustration we examine the construction order for a cir- 


cuit using FE-B, The operating conditions for the circuit are given in 
the switching table (Fig. 4a), 


1. In accordance with the above rules we switch in the secondary 


input signals x, and x, necessary for realizing the given conditions 
FE, (Fig. 4b). 


FE 2. Having assigned weights to the input signals: x,-1, @-2, 
: X3—4, x4—8, we obtain the necessary and conditional selection of num- 
bers for the output signals zy and 2, [10]: 





Fig. 3 


Z, = 3 (2, 6, 7, 10, 14, 15), 
% = 11, 13 (2, 6, 7, 10, 14, 15), 
Here the conditional numbers are the unused numbers. 


3. We consider two cases when for FE-B we use: a) triggers with separate inputs, and b) triggers with de- 
numerabie inputs. 


We obtain the necessary and conditional numbers for case (a). The signals reacting on the trigger must be 


separated into two groups: signals switching in the trigger (y'), and signals switching out the trigger (y"). We can 
determine the corresponding functions F* and F*. 
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From the table in Fig. 4c we find the necessary number of functions Fj: 
F}=3. 
Besides the unused ones, the numbers 4, 5, 12 are also conditional numbers for Fj because the appearance of 


signal ys in any sequence up to the moment the signal y; appears does not change the signal at the output of this 
trigger. 
Thus, 
Fj = 3 (2, 4, 5, 6, 7, 10, 12, 14, 15). 


Analogously we can show the equivalent states for functions Ff and Fj (since the signal x, is not switched out, Fj 
is not defined). 


FZ = 13 (0, 2, 6, 7, 8, 9, 10, 11, 14, 15), 
F, = 5 (2, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15). 
Besides the numbers 8, 9, 11, 12, 13 as conditional numbers for function F4 we can use the numbers 0, 1, 3. 
But, when considering the necessary number 1 it is also necessary to use the number 3 because otherwise the signal 


y4 will appear at the transition from state 1 to state 3, i.e., the trigger will switch in prematurely. Thus the num- 
bers 0, 1, 3 are interconnected. 





Fig. 5 


Since the use of such connected numbers complicates the determination of the corresponding function F, here 
and in what follows these states will not be considered. 


For case (b), when we use the trigger with a denumerable input, the switching-in and switching -out signals 
enter at one and the same input of the trigger, and therefore for every trigger we can find only one function F (Fig.4d). 


In this case every signal reacting at the trigger changes the signal at its output, and therefore, besides the un- 
used ones we need to take only the connected numbers as conditional numbers. Since these numbers are not to be 
used, then 
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4d), 


Be 





Fs = 3, 13 (2, 6, 7, 10, 14, 15) 
F,= 5 (2, 6, 7, 10, 14, 15), 


4. In each of the two cases we simplify the determination of the functions z and F with the use of the con- 
ditional numbers [10] and then construct the circuit from the elements AND, OR, NOR, and d, Figure 5a shows the 
circuit realizing the conditions given in the table of Fig. 4a for case (a), and Fig. 5b for case (b).” 
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Fig. 7 


It is interesting to note that for realizing the conditions given in the same switching table (Fig. 4a), when 


using FE-A it is necessary to introduce three feedback signals (Fig. 6a), and the entire circuit obtained becomes con- 
siderably complex (Fig. 6b). 


Thus the use of FE-B leads to the possibility of a simple circuit in a number of cases. 


In conclusion, however, we note that such a circuit representation (Fig. 3) is not always possible. It is known 
(6, 15] that for some potential-impulse forms we cannot find a unique function F whose differential equals y, Here 
sometimes y can be represented by 


y= VdF) G=1,...,m). 


In this case the circuit will have the form shown in Fig. 7. 


* In the diagrams we represent the operator OR (realizing Boolean addition) by the symbol + in a rectangle, and the 
operator AND (realizing Boolean multiplication by the symbol « in a triangle, 
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THE EFFECT OF PULSE NOISE ON TELEMETERING UNITS 
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(Moscow) 
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pp. 1202-1209, September, 1961 

Original article submitted January 23, 1961 


The paper performs a theoretical investigation of the nature of the effect produced by pulse noise 
with a logarithmically -normal amplitude distribution law on linear devices of two types and on a 
system consisting of a limiter and a filter. Simple computation formulas are derived for deter- 
mining the amplitude and frequency characteristics of the noise. 


Recently an ever increasing amount of attention has been devoted to the investigation of the interfering effect 
of pulse noise, in addition to the studies which are being performed on fluctuating noise. However, heretofore the 
problems involved in the noise stability of communication and remote control systems which are subjected to actual 
pulse noise have been inadequately studied. In the majority of published papers the analysis has been performed for 
idealized pulse noise signals which do not conform to actual ones. Such an approach can be explained by the absence 
of adequate data on statistical characteristics of pulse noise in practical channels. 


In addition, the investigation of noise performed by a number of authors [1,2] has made it possible to conclude 
that many types of industrial noise can be included in the class of pulse noise signals with a logarithmically -normal 
amplitude distribution law. This type of noise exists, in particular, in low-voltage power nets (0.4/6 kv) which are 
used for a number of industrial objects [1, 3]. It should be noted that the noise stability for transmitting information 
over such channels becomes especially important in view of the high noise level and the large attenuation. 


In computing the noise stability of remote control and communication systems when pulse noise is present it 
is necessary to determine the statistical amplitude, time and frequency characteristics for the noise output of various 
linear and nonlinear devices, as well as the energy spectrum of these noise signals. As was demonstrated in [1], the 
measurement of the amplitude distribution of pulse noise in a channel is conveniently performed using a selector 
unit with an inertial detector and a pulse counter incorporating a threshold relay. For noise signals with a logarith- 
mically -normal law the amplitude distribution is characterized by two parameters which can be determined by 
processing the data at the output of the read-out device [1]. It is of interest to develop methods for computing the 
noise stability of various devices on the basis of data obtained by means of the read-out device. 


The paper studies the effect of pulse noise on two types of linear units: a narrow-band filter at whose output 
the noise becomes fluctuating noise, and a wide-band filter for which the noise at the output becomes pulse noise. 
In the first case the interfering effect of the noise can be adequately evaluated on the basis of the noise power; in 
the second case it is necessary to determine the parameters of the amplitude distribution as well. 


We shall study a nonlinear device consisting of an inertialess limiter; we shall develop a method for com- 
puting a systern consisting of a filter, a limiter, and another filter. The analysis of the passage of pulse noise through 
the enumerated devices requires a more detailed study of the amplitude characteristics for logarithmically -normal 
distribution. We have derived formulas for computing the distribution moments, as well as the expression for the 
characteristic function. 


1. The Parameters of Logarithmically-Normal Distribution 
~The distribution of the amplitude for pulse noise signals in the class indicated above can be written convenient- 
ly on a logarithmic scale (Fig. 1). For u = In x 





(u—ny 


mavre & x0! 1 
@(u) = Vie e ‘ (1) 
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where u = In a and o respectively denote the average and mean-square values of the distribution u. The parameters 
u and o fully characterize the distribution u. As we shall demonstrate below, all of the basic characteristics of the 
distribution x can be expressed in terms of these two parameters. 
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Fig. 1 


It is not difficult'to show that the quantity o is independent of the units in which the noise amplitude is 
measured and is a dimensionless quantity. The quantity U is dimensional, and when the scale varies by a factor of 
k the value of u varies by the amount In k. Note that the difference u—U, just as o,, is dimensionless. 


In the case where we study a logarithmic distribution with the base b instead of natural logarithms, the average 
and mean-square values are increased by a factor of log}, e and the abscissas of the distribution are reduced by the 
same factor. 


Making use of the conventional laws for nonlinear transformation of the distribution law, we find the distribu- 
tion for x 


_ (inx—Iinay 


w(x) = Vixen e Were (2) 





The characteristic points on this distribution (Fig. 1) are the median which is equal to a, and the abscissa of 
the most probable value X;ygx (the mode). From the condition w' (x) = 0 we obtain 


Zmax = ae”. (3) 


The moments for the distribution x can be determined conveniently in termsof the derivatives of the charac- 
teristic function which is equal to 





0, (v) = m, {e**") = is \ ettete oF as (4) 
for a logarithmically-normal law. 
The kth order moment which is determined from the condition 
Be | = itm, {2}. (5) 
is equal to 
m, = we? (6) 


The distribution moments of the first two orders are of greatest significance in describing the distribution x; 
these are equal to 
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m, {x} = Z = ae*?, (7) 
Mz {xz} = a*el2e)"/2, (8) 
The dispersion of the distribution is written as 
M, {x} = Dz = 62 = ae” (e* — 1). (9) 
For certain special applications it is convenient to make use of a characteristic function which is written in 
the form of a MacLauren series in accordance with (5): 
(ko)? 


A my, . mage tier Tae —- 
t2(v) = >) (io = De®. (10) 


k=0 k=0 


2. Parameters of the Clipped Distribution 

In analyzing the passage of pulse noise through nonlinear devices of the inertialess limiter type it is necessary 
to determine the parameters of the clipped distribution. We shall determine the first two distribution moments for 
pulses which have passed through a limiter which limits top and bottom at the level Ug (Fig. 1). 





The average value of the amplitudes for the peaks exceeding Uy can be determined from the formula 


{ exp [- (In z — Ina)? Jaz 














tis . 253 (11) 
‘Di oep 1 exp [-" z-- me) - 
i= yess 
we shall introduce the relative cutoff level 
a In /9 —ina (12) 
Fi 
Substituting (12) and (7) into (11), we obtain 
—- 2 (a—s) (13) 
b V (a) : 
where V(a) is the tabulated function 
co i 
7 eal \ 2 5 
V (a) = Er e dt. (14) 
In analogous fashion we find the average amplitude of those pulses which do not exceed Up: 
= — tIi—Via—>)} (15) 
¥ 1 — V (a) 
The second order moment is found from expression (2). For pulses above the threshold 
_ mV fa — 2s) (16) 
ms, = V (a) 
For pulses below the threshold 
ms [1 —V (a — 23)} 
I sc an 
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The dispersion of the clipped distribution of the noise peak amplitudes is equal to 


on __ mo — mV (a—26) — mV3 (a — 99) 
and >, my My V (a) Vi (a) (18) 





ms [1—V¥ (a —20)) ™(t—V(a—o)) (19) 
1—V (a) [1—V (a, 





Dy = Mey — oe = 


respectively. 


3. The Effect of Pulse Noise on Linear Devices 

We shall study the case for which brief pulses of arbitrary but identical shape with a random instant of appear- 
ance and a random amplitude u,, exist in a channel. We shall assume that the distribution function for the pulse 
amplitudes obeys a logarithmically -norma! distribution law, and that the spectral density for a pulse of unit am- 
plitude is equal to G(w). 


Assume that the noise signals in the channel are measured by a selective read-out device with a frequency 
response f (w)® and with a pass band such that superposition of pulses at the output can be neglected. The amplitude 
of the noise peak at the output of the read-out device is related to the amplitude of the peak in the channel by a 
simple expression 





L= Pa G(@) /, (@) do. (20) 


If within the limits of the pass band of a selective read-out device the spectral density G(w) is constant and 
equal to G(w») (where we is the tuned frequency of the read-out device), then 


x = G(@) A 1, u (21) 


(It is assumed that the read-out device is graduated in such a way that its gain can be assumed equal to unity.) 


It follows from formula (21) that the amplitude of the peak at the output of the read-out device proves to be 
proportional to the equivalent pass band Af, of the read-out device which is expressed as 


12 (22) 
Af, =32\ f, (odo. 


Assume that we know the parameters of the amplitude distribution 1 and o which are determined at the output 
of the read-out device. We shall demonstrate how these parameters vary when a random pulse process passes through 
a linear device. 


Assume that a receiver with an input filter which has a frequency response z(w) that differs from the response 
of the read-out device but is such that pulses at the filter output do not overlap is connected to a communication 
channel. The transfer coefficient for the filter is assumed equal to kg. The amplitude of the noise at the output of 
such a filter when a pulse up is applied to its input is equal to 


u = ky7'z, (23) 
where 
t G (@)z(@)da@ 





J G@)/@do 


* The frequency responses are assumed to be normalized. 
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Since the factor multiplying x does not depend on the amplitude of the peak we can conclude that for linear 


devices the amplitude distribution law for the noise remains the same when the individual pulses do not overlap. 
The distribution parameters can be expressed as follows in accordance with (2): 


Inue = Ina+ Inky’, 


(24) 
or= SG. 


(25) 


If the spectral density of the noise is uniform within the limits of the pass bands of both the receiver input 
filter and the read-out device, then the quantity y is equal to the ratio between the equivalent pass bands of the 
receiver and read-out device and 


# ko A fe 


uc 
f ’ 
4/. 


where Af, is the equivalent pass band of the receiver filter 


(26) 


Afp=x§ z(@)do. 


(27) 
As an example of the second case of linear transformation we shall study a narrow-band filter with a frequency 
response z(w) such that the noise at its output can be treated as fluctuating noise, Under these conditions it is ex- 
pedient to evaluate the noise signals according to their power. 


The energy spectrum for the noise at the filter output is equal to 


F (@)=W (@) z* (a), 
where W(w) is the energy spectrum for the noise in the channel. 


(28) 
amplitude is given by 


We know [4] that the energy spectrum for a chaotic train of pulses that are of identical shape but of different 


OTE 2n x (29) 
- 2 —— 
W (@) = 2n|u ” (@)| af , 

where n is the average number of peaks per second, Substituting (29) into (28) and integrating within the limits of 
the filter pass band, we obtain the following expression for the noise power at the filter output: 





2nA Fi¢ zi 
ef = 
P = > Se . (30) 
pf Af, 
After substituting (8) and (30), we obtain 
2n AF ¢ a2 ~o* 
P of = 





(31) 
Af é 
r 
where AF ef is the energy width of the filter pass band: 


1¢ 
AF =~) 2* (w) dw. 
—o 


(32) 
For the same conditions as those which apply in formula (26) we obtain the expression for the noise power at 
the filter output: 
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AF os 
Por = Por aT.’ (33) 


where Por and Af,, respectively denote the noise power at the output of the read-out device and its energy band 
width. 


Note that in scaling the amplitude characteristics of the pulse noise we use the ratio between the equivalent 
pass bands of the filter and the read-out device; in computing the power we use _ the ratio between the energy (ef- 
fective) bands. Thus, in computing the passage of pulse noise through linear devices it is necessary to know the 
parameters of the amplitude distribution at the output of the selective read-out device, as well as the equivalent and 
effective pass bands of the read-out device. 


4. The Passage of Pulse Noise through a Device with a Nonlinear Element 

We shall study the effect of pulse noise with a logarithmically -normal amplitude distribution law on a receiver 
consisting of a relatively wide-band filter 1 with a pass bandA i and a gain Ko, an inertialess detector, 2, a noise 
limiter 3 which clips the top of the wave and has a limiting threshold Up, and a narrow-band output filter 4 with a 
pass band AF and a gain equal to unity (Fig. 2). 


~THEHEHE+ 


Fig. 2 











In performing our computation we shall assume that the pass band Afof the input unit is such that a) the 
superposition of pulses beyond the detector can be neglected, and b) the output filter has such a narrow band that 
the noise at the output of the receiver becomes fluctuation noise. In that case we can assume that the spectral den- 
sity for the pulses after the detector remains uniform within the limits of the output filter pass band. 


The energy spectrum and the parameters of the amplitude distribution for the pulse noise after the detector 
are determined from formulas (28), (24), and (25); the parameters of the clipped amplitude distribution after the 
limiter which clips the top of the wave are determined from formulas (15) and (19). In accordance with these for- 
mulas we can write the amplitude of the pulse after the detector as 


y = rkoXo, (34) 
where 

Af 
} 35 
To = ay . (35) 

r 

Under these conditions 

(36) 


Iny = Ina +- In ko¥o- 


The energy spectrum for the pulse noise after the limiter which clips the top of the wave can be expressed in 
the form of two terms: one for pulses whose amplitude is less than Up, and the other for clipped pulses in accordance 
with formula (29): 








Fy (o) - 2n{ : ("y1@(0)\* xp[— Co etal dy 








ae 262 
V 2nc 3 ad 
4 i ( |G, (@, Lee |- (In y — In akyy,)? | dy! (37) 
V 2xs y P 23? i 


Here G(w) is the spectral density for a pulse of unit amplitude y = Up; G,(w, y) is the spectral density for a 
clipped pulse y > Up. 
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Finding the energy spectrum in general form from formula (37) is difficult because the spectral density of the 
clipped pulses Gy(w, y) is a nonlinear function of w, y and Up; the shape of this nonlinearity is determined by the 
shape of the input filter frequency response. rs 








Fig. 3 


We shall study the case where the input filter has a frequency response of the single tuned circuit type. Under 
these conditions the pulse noise after the detector (Fig. 3) consists of a series of exponential pulses with a random 
amplitude and an identical shape: 


u (t) = ye—4t, (38) 
where Aw is one-half the pass band of the input filter at the 0.7 level. 
The spectral density of an exponential pulse with the amplitude y at zero frequency is equal to 
(0) — 39 
G(0) = Ao’ (39) 
The spectral density of the portion of the pulse that is clipped from the bottom is equal to 


G, (0) = Ct (40) 





and for a pulse clipped from the top (the shaded portion of Fig. 3) it is 


G (0) = G (0) —G, (0) = Vor + F*. (41) 
Substituting the expression 

{ D 

t= a5 ln7, (42) 
into (41), we obtain 
ey 

Ueln 43 
G,(0) = ——— Ze. (©) 


Ao 


The expression for the energy spectrum of the pulse noise after the limiter which clips the top of the noise at 
zero frequency can be written as 














U, 
_ 2n i __ (Iny— Inak,y,)* 
Fo(0) = gor | aay yexp| sp | dy + 

co (in 

us ¢ (ng) (in y — Inte) (44) 
+ Vine \ pe exp [— SB ates | ay}. 
Introducing the relative clipping level 
oo 2 SS (45) 
6 
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we obtain the following result after integrating expression (44): 


Uzoe—¥* (95 — 2) 


(46 
Viz ’ 





F,(0) = 5 (>. + UW (¥) (s* + (oy—1)"] — 


The noise power at the output of the narrow-band filter 4 (Fig. 2) with the frequency response C(w) is deter- 


mined from the formula 


1. 
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P= Fo(0) 5- | C(w) do = Fo (0) AF « (41) 


—0o 


The author thanks G. A, Shastov for his useful remarks offered in the discussion of this paper. 
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pp. 1210-1219, September, 1961 

Original article submitted April 13, 1961 


The article analyzes the efficiency of the transmission of telemetered information for the 
basic forms of pulse modulation by the criterion of specific rate of information transmission. 
Taking account of the equipment errors of signal transformation makes it possible (in part I 
of the article) to carry out the analysis on the assumption that there is no interference in the 
communication channel. This corresponds in a number of cases to the practical conditions. 

A comparison of the method of transmitting telemetered information is made and is 
illustrated by data on concrete telemetering systems. 


By the efficiency of information transmission we generally mean some complex estimate of the economy or 
optimality of transmission on the basis of all the factors influencing transmission. These latter, of course, in- 
ctude the effects of interference in the communication channel. In the intial stage, however, it is possible to define 
and study the efficiency of transmission of telemetered information in the absence of any interference in the com- 
munication channel, since even in this case the basic parameters which characterize the transmission are in relations 
which can be perfected (optimized) according to a chosen criterion of efficiency. In addition, the transmission of 
telemetered information under conditions where interference plays only a small role is often found in practice. 


We call attention to the fact that in what follows equipment errors are taken into account in the analysis of 
transmission efficiency, which makes it meaningful to carry out such an analysis even in the absence of interference 
in the communication channel.* The unified approach to the analysis of continuous and discrete methods of trans- 


mission is based here on considering the error caused by equipment instability to be identical with the discreteness 
error. 


For the estimation and comparison of the efficiency of information transmission for various methods of tele- 
metering we shall make use of the criterion of “specific transmission rate,” which indicates the ratio of the trans~- 
mission rate R to the bandwidth W occupied in the communication channel. This convenient criterion was proposed 
in [1]. 


In the absence of interference and for a uniform distribution of the information transmitted, the transmission 
tate may be represented by the ratio 


Yad ’ 





R= 


where Na = 1/5, is the precision of the equipment or the number of quantizing levels of the measurement, the 
reciprocal of the equipment error 64, and T is the time necessary for the transmission of one measurement (the 
telemetering response time). Here and in what follows the symbol "log" will mean logarithm to the base 2. 


* By contrast with a number of studies devoted to this question (see, for example, [2]) we have here a distinction be - 
tween the concepts of “interference” (noise), which arises in a communication channel, and “equipment errors," 
which are characteristic of the equipment transforming one form of informations carrying signal into another. The 
necessary presence of equipment error means that even in the absence of transmission interference only a finite 
number of quantizing levels of the initial information are possible. 
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The specific rate criterion in the present case can be expressed by 


R log N 4 2 
(i) =" (2) 


The denominator TW* may be called the frequency inefficiency. The greater TW is for the same response 
time the more poorly the selected bandwidth is utilized in the communication channel. The frequency inefficiency 
TW is of interest in a number of cases as an independent criterion, the reduction of which will lead to the perfecting 
of the telemetering instruments. 


We shall derive explicit expressions for the rate R and the specific rate R/W of transmission for the basic forms 
of modulation used in telemetering. 


1. Pulse-Code Modulation (PCM) 

Let us consider a coded pulse train with n code pulses. The duration of the pulse train r consists of the length 
of the code pulses r,, and the pulse separation rs, as well as the synchronizing pulse. Let us assume that the syn- 
chro pulse and the separations have duration equal to (Tp + T,). Then the duration of the pulse train is defined by 
the formula 





T= (n+ S)(t, + i (3) 


We shall express the duration of the front + ¢ by means of the required bandwidth W, in the communication 
channel: 


and we shall relate the total duration of the code pulse and separation Tp + 7, to the duration of the front r¢ by 
the formula 


° 2k 
tpt % = 2katr— 7, (4) 


where kj =>1. The limiting case k; = 1 can be represented, for example, as transmission by closely adjoining ti- 
angular pulses. 


For PCM the transmission time may be considered identical with the duration of the pulse train, that is T = T. 
On the basis of (1), (3), and (4) we define the transmission rate by 





7 log N 4 log Ny (5) 
Reco = (n + S) (+ tT.) 2k (n + 8) W, 
(for a binary code log Na = n.) The expression for the specific rate (2) will assume the form 
logN, W, (6) 


R\ ew! 
(7) peu hats) W 
We should remark that the bandwidth W actually represented may coincide with the required bandwidth W, or 
may be greater. 
In one channel transmission in most cases S = 2. 


in multi-channel transmission the relative weight of the synchro pulse for an individual channel is decreased, 
and for a large number of channels the synchro pulse may be neglected; then the expression for relative rate for a 
binary code will have the form 


* TW: 2f)/2F where F is the upper frequency of variation of the parameter measured, and f) is the limiting fre- 
quency of the low-frequency filter in transmission without secondary modulation. 
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(WF) par Ry: (6a) 
PCM 2k,W 


2. Pulse-Width and Pulse-Time (Position) Modulation (PWM, PTM) 

We shall first consider transmission by pulses with steep fronts. Such transmission is necessary when no special 
measures are taken in the reception to eliminate the relative instability of the limiting level and the amplitude of 
the video signal (see below), Let us assume that in this case the duration of the front r¢ is equal to the permissible 
absolute error of the measurement. 





The duration of one cycle is defined by the expression 


a as Nate (1) 
t Tat 2t in (Na + 2k), See a 


Here Tq = Nr¢ is the time deviation, T)4,, = kerf is the pulse (or separation) of minimum duration in PWM (equal 


to the working pulse in PTM), 


The time T in which the readings are established in the telemetering receiver in PWM (or PTM) may be equal 
to one or several measuring cycle +, depending on the method of demodulation. 


In the general case 

T = at. (8) 
We shall call a the inertial coefficient.” 
The expression for the transmission rate on the basis of (1), (7) and (8) can be represented in the form 


log N, Ww (9) 
PWM(PTM)~ a(V+2h) °° 
By criterion (2) we obtain the formula 


(=) 8 Na Wo (10) 
W/PWM(PTM) @(N, +2k)W * 


R 


In a number ofcases it is possible to carry out narrow-band transmission by PWM (or PTM) that is, transmission 
by pulses with gradual fronts. Such transmission is more economical by the R/W criterion, but it necessitates some 
additional complexity in order to avoid errors in the measurement of time intervals as a result of a relative instability 
of the amplitude of the video signal and the limiting level. For this purpose the receiver uses either automatic regula- 
tion of the signal level or limiting of the video pulses by the “tracking” level, that is, a level proportional to the 
amplitude. In addition, the change to transmission with very gradual front sometimes requires the application of a 
more stable modulator as a result of the decrease in time deviation (if the same cycle duration is kept), When these 
conditions are satisfied, the duration of a cycle can be conveniently represented in the following manner: 


ai ain ai — Nb + 2hs 
t titer NA + 2kgt, (Nb + 2k) t, “Wo (11) 


where A is the absolute error of the measurement (by contrast with the previous case, here A < rf), b= 4/r¢ is 
the relative steepness factor of the front. 


On the basis of (1), (8) and (11) the expression for the transmission rate can be written in the form 


* For example, if the reception is carried out by the average current of the pulses of PWM (PTM is easily transformed 
with demodulation into PWM) calibrated with respect to amplitude, then a > 1, If the decording of the readings 
takes place during one cycle, then a = 1, 
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R log N , 


PWM(PTM) a(N ,b + 2h) Wo (12) 





and ‘the specific rate criterion will take the form 


(=) 4. log N 4 Wo (13) 
Ww PWM(PTM) a(N ,6+-2k:) W ° 





3. Transmission with Automatic Readjustment 

Transmission with automatic readjustment of the response time and the precision [3] is performed with PCM, 
PWM, and PTM. Let us remember that in those cases when precise telemetering of slow changes of the parameter 
is required but at the same time a cruder transmission of the rapid changes of the parameter is permissible, it is 
possible to use automatic readjustment of the transmission. In automatic readjustment an exchange takes place be- 
tween the response time (duration of the pulse train) and the number of transmitted quantizing levels, depending on 
the rate of change of information. As a result, the rate of information transmission is decreased. If we compare the 
ordinary telemetering devices with the readjusting type, we find that in the latter case there is a gain in transmission 
efficiency, since we require a channel calculated for a lower transmission rate. It is not significant, however, to 
evaluate this gain by means of the specific rate criterion, since transmission with automatic readjustment is charac- 
terized by a lower rate R and a narrower bandwidth W; as a result the rate either does not change (in PCM), or 
changes by an insignificant amount (in PWM and PTM) by comparison to systems without readjustment. 





In connection with this, for comparison of the efficiency of different methods of transmission, including meth- 
ods with automatic readjustment, we introduce as a uniform criterion for readjusting transmission the concept of 
"quasi-rate” of transmission. Let us remember that a readjusting device is characterized by a maximum precision 
Nmax and a minimum transmission time Ty)jin. We shall represent the quasi-rate of readjusting transmission by the 
expression 


= = max 

q min 

The value of Nmax and Tin coincide with the values of N and T for a non-readjusting system. in other 

words, in using the concept of “quasi-rate" we make the assumption that the adjust system transmits the same in- 
formation about the measured parameter as was transmitted earlier by the ordinary non-adjusting system, that is, 
that the rate of information transmission has not changed. In addition to this, changing to a readjusting system led 
to a reduction of the bandwidth necessary for transmission on the communication channel. Therefore, the efficiency 
of readjusting transmission can be evaluated by the ratio of the quasi-rate of transmission to the bandwidth occupied. 
This makes it possible to compare the efficiency of adjusting and non-adjusting transmission by the uniform criteria 
Rg/W and R/W, respectively. The criterion Rg/W may be called by analogy the specific quasi-rate of transmission: 


Rq “a log N ; (14) 
Ww Ee 
min 


4. Pulse-Frequency Modulation (PFM)* 

Let the tracking frequency of pulses in telemetering with PFM change from f; to fy. The lower frequency f,, 
in general, must be so chosen that the condition of telemetering response time is satisfied." * However, we must 
consider the difficulty of producing low-frequency pulses. A very convenient criterion of the degree of realization 
of a given response time T is the inertial coefficient which in the present case is equal to 





* The material of this section is equally applicable to frequency modulation as well, if the latter refers to “low-fre- 
quency telemetering,” that is, a low-frequency change with large relative deviation. However, in practice low- 
frequency systems are usually pulse-frequency systems, since for low-frequency it is easier to form pulses than 
sinusoidal oscillations. 

* * For example, if the reception of pulse~frequency telemetering is effected by the average current of short pulses 
calibrated with respect to amplitude and duration (as in the widely used capacitor-type frequency meter), then the 
frequency f, must be matched with the receiving equipment so as to exclude pulsations in the readings of the instru- 
ment. 
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a=T};. (15) 


Furthermore, the upper frequency f, is defined by considerations of the precision of telemetering. It is obvious that 
the precision may be defined by the formula 


_ fa—h 16 
Na=-=z": (16) 
where A is the absolute equipment error. 
After some simple transformation, we obtain from (16) 


ee . 
N,zo +1 
he 


Using this relation and introducing the notation c = A/f,, we can express T from (15) 
T = (Nac + 4). 


According to (1), the transmission rate is equal to 


log Ny 


Roem = a(N,e+1) /* 


(17) 


We note that Nac + 1 = fz/f,. Obviously the upper frequency f, defines the necessary bandwidth W, in the 
communication channel, For narrow-band telemetering in a communication channel with secondary modulation it 
is sufficient to transmit two side frequencies, that is, Wy = 2f,. Then 


(=) __'8N, We (18) 
Ww PFM 2a(N,c+1) w 


5. Analysis and Comparison of the Efficiency of Telemetering Transmission Methods 
The resulting expressions R/W are, as can be seen, functions of a number of arguments. An analysis of the 
efficiency can, conveniently be made by means of graphs, choosing in the present case, that is, in the absence of 





R/w, binary units 














Fig. 1. Transmission efficiency as a function of telemetering 
precision for PCM and PWM (PTM). The PWM (PTM) trans- 
mission is of the wide-band type, that is, it uses pulses with 
steep fronts for b = 1. 


interference in the communication channel, the equipment precision Nag of the transmission as an argument. For 
simplicity of analysis we shall consider W,/W = 1. 


In Fig. 1, using formula (6), we have constructed for a binary code (log N = n) the curves of (R/W)pcyy 48 & 
function of Na for S = 2, with k, = 1 for curve 1 and k, = 3 for curve 2. As the precision Na increases, the efficiency 
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of coded transmission increases, since the “weight” of the sync. pulse, which does not carry useful information in 
the coded pulse train is decreased. These curves approach the dashed straight lines (S = 0) constructed according to 
formula (6a). From a comparison of curves 1 and 2 it can be seen that an increase in k;, that is, an increase in 
pulse duration as compared to the duration of the fronts, decreases the transmission efficiency, since in this case 
there is an increase either in the time T (for W, = constant), or the necessary bandwidth W, (for T = constant). 


We shall indicate on this graph the values of R/W for concrete telemetering devices. In the Appendix we list 
the data for a number of devices; it follows from these data that the efficiency of the coded devices developed by 
the VNIIE(TK~-2) [4] and the TsNIIKA [5], are defined respectively by points 1 and 3 in Fig. 1. 


We now pass to the case of PWM (PTM). The (R/W)pwM(PTM) Curves calculated by formula (10), charac- 
terizing transmission with fairly steep fronts (b = 1), are also constructed in Fig. 1: curve 4 for ke = 1 and curve 5 
for ke = 3 (in both cases a = 1), These curves have a maximum whose physical meaning is the following: At the 
beginning, for small Na for the duration of the entire pulse train great specific weight is attached to the minimum 
pulse (and separation), which does not in itself carry any useful information, The increase of Nag at first increases 
the efficiency, since there is an increase in the time deviation rg (transmitting useful information) by comparison 
with the minimum pulse, but after this the transmission efficiency decreases, tending to zero, since the duration of 
the entire pulse train for Wy = constant (or the bandwidth W, for T = constant) increases faster than the number of 
binary units transmitted by the pulse train.* We note that the difference in the values of ky is pronounced only in 
the case of small Na: for Ng = 50, curves 4 and 5 become practically identical. It is obvious that an increase in 
the inertial coefficient (curve 6 for ke = 3 and a = 3) sharply reduces the transmission efficiency. 


Examination of the graphs in Fig. 1 shows that the transmission of information by PCM for k; = 1 is always 
more efficient than wide-band transmission (that is, transmission by pulses with steep fronts) using PWM (PTM). 
However, for k; = 3, coded transmission becomes more efficient than wide-band transmission with PWM (PTM) only 
for large values of Na. 


Figure 2 shows graphs constructed by formula (13) for narrow-band transmission by PWM(PTM) for several 
values of the relative steepness factor b of the fronts: from b = 1 to b = 1/20 (the solid lines refer to a = 1, the 
dashed lines to a = 3; ky = 3 in all cases), For comparison we also show here curves 1 and 2 for PCM for the values 
ky = 1 and k; = 3, respectively. 


It is clear from the graphs of Fig. 2 that transmission by PWM (PTM) with b = 1/2 is more efficient than 
transmission,by PCM for k; = 3 (curve 2), while PWM (PTM) with b = 1/20 is more efficient even than PCM for 
ky = 1 (curve 1) for a wide range of values of Na. We note that the decrease of b not only increases the efficiency 
but also displaces the maximum of the efficiency into the region of greater N, values. 


Thus, comparison with coded transmission shows that the transmission of information by PWM (PTM) using 
pulses with gradual fronts is very efficient. Unfortunately, at the present time no PWM(PTM) telemetering devices 
have yet been developed which realize this type of transmission. Existing telemetering devices, for example, those 


* Of course, such a type of functional relation is dictated by the choice of the logarithmic —in the present case, 
binary —units of measurement of information imposed in the definition of transmission rate R. 


Besides this there is another possible estimate of the transmission rate of telemetered data—by the number of 
quantizing levels transmitted per second. The expression for the transmission rate will then have the form Va = Na,/T, 
and the graphs of specific rate Va/W versus precision will look different. For example, the graph of (Va/W)pcy 
will increase monotonically as Na increases, while (VA/W)pwM(PTM) and (VA/W)prm will tend to a finite limit 
different from zero.. 


In a number of cases, especially in the analysis of transmission efficiency where interference is taken into con- 
sideration, it is necessary that the changes in precision N be reflected more clearly in the graphs of the efficiency 
criterion as a function of such transmission parameters as the bandwidth in the channel, interference level, etc. Under 
these conditions an estimate by the criterion V/W = N/TW may be found to be preferable to the R/W = (log N)/TW 
estimate in which the logarithm in the numerator smooths the effect of changes in the precision N (under the in- 
fluence of interference the transmission precision N is determined both by the equipment error and by the error caused 
by interference in the communication channel). 
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of Bristol or the VST-1 (IAT AN SSSR Institute of Automation and Remote Control of the Academy of Sciences of 
the USSR) [6], belong to the class of wide-band devices; their efficiency is defined respectively by points 3 and 4 
in Fig. 2 (the data for these devices are given in the Appendix), 


k/w, binary units 





a6 








~ —— | 
EE SY ee 





teres See ae eS iS _-> = 
eee eT ee “BO 200%, Hantizing 


Fig. 2. Transmission efficiency as a function of telemetering 
precision for narrow-band PWM-PTM. Curves 1 and 2 for PCM 
are shown for comparison, The solid curves refer to a = 1, the 
dashed curves to a = 3. 


k/w, binary units 
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Fig. 3. Transmission efficiency as a function of telemetering 
precision for PFM, For comparison the figure shows curves 1 
and 2 for PCM, Curve 3: a= 15, c= 0.04 (ChTI-1); curve 4: 
a= 5, c= 0.18 (ChIS-1); curve 5: a= 12, c = 0.08 (ChI-1); 
curve 6: a= 81, c= 0.012 (TNCh). 


We shall show for comparison the data on the efficiency of PWM with automatic readjustment (PWM-AR), 
A breadboard model of such a device was made at the(Institute of Automation and Remote Control of the Academy 
of Sciences of the USSR); the data for this breadboard model given in the Appendix. The PWM-AR device transmits 
pulses with fairly steep fronts (r¢ = A, that is, b = 1), but by virtue of its automatic readjustment it has a fairly high 
transmission efficiency. Characterizing the transmission efficiency of PWM-AR by the specific “quasi-rate” of (14), 
we obtain point 5 in Fig. 2. 


We shall now analyze PFM. Formula (18) is analogous in structure to formula (13), and therefore in the case 
of PFM we should expect curves similar to those shown in Fig. 2. Since PFM devices have become widespread in 
practice and a number of concrete developments are known today, it is interesting to construct the efficiency curves 
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for coefficients a and c whose values are close to the real values. In the Appendix we give the data which make it 
possible to determine these coefficients for four telemetering devices: ChTI-1 (IAT AN SSSR) [7], ChIS-1(TsLEM) 
[8], ChI-1 (TsNIIKA) [9], TNCh (TsNIEL—Elektropult) [10]. In these data it is assumed that all the devices have the 
same precision Na = 50 (which is close to the real value) and occupy channels with frequency bands of the minimum 
width necessary for passing the first harmonic for an upper frequency of f,* (that is W9/W = 1). 


Figure 3 shows the curves for the following values of the coefficients: curve 3: a= 15, c = 0.04 (ChTI-1); 
curve 4: a= 5, c = 0.18 (ChIS-1); curve 5: a = 12, c = 0.08 (ChI~-1); curve 6: a= 81, c= 0,012 (TNCh), It is easy 
to observe that the lower the value of c, the farther to the right is the maximum point of the curve. All four curves, 
for the values of Na useful in practice, are much lower than the “standard” curve for coded transmission for ky = 3 
(curve 2), not to mention curve 1 for k; = 1. 


The comparatively low efficiency of pulse-frequency transmission is determined chiefly by the comparatively 
large values of inertial coefficient a. Let us take as an example curve 4, constructed for the coefficients of the 
ChTI-1, If we replace a = 15 by a= 5,°* we obtain curve 7, which up to the value N, = 100 is higher than the 
“standard” curve 2, 


Thus, the lower the value of a, the higher is the transmission efficiency. In the limit a = 1, that is, the most 
efficient reception of readings would be in one period of the lower frequency. Such a reception method is possible 
in principle, but it is obviously analogous to reception by PWM in one cycle, ** * that is, at least with respect to 
the method of reception, PFM would reduce to PWM, 


In addition, transmission efficiency increases as c decreases, or which is the same thing for fixed precision 
NA as the factor Nc + 1= f/f; in the denominator of (17) decreases. Clearly, a decrease in c or in f:/f; (for a given 
value of f;) is equivalent to a decrease in the frequency deviation or the absolute transformation error A, 


A decrease in the inertial coefficient a and the ratio fy /f, is nothing but a decrease in the components of the 
denominator in (17), that is, a decrease in the frequency inefficiency TW, since 


2a(Nac+1)(7) = 2a z (=) = TW. 


In the four telemetering devices taken as examples the transmission precision is the same (Na = 50), and 
therefore the comparative transmission efficiency of each of them (see points 3, 4, 5, and 6 in Fig. 3) is determined 
by the values of TW. 


Summarizing the above, we observe the following facts: 


1. The resulting expressions for the R/W criterion as a function of all the fundamental electrical parameters 
characterizing telemetering transmission make it possible not only to estimate and compare with great clarity the 
transmission efficiency of concrete telemetering systems based on different methods of modulation but also to show 
ways of perfecting these systems and the conditions of their most efficient operation. 


2. A comparison of the telemetering transmission methods PCM, PWM (PTM) and PFM (in a channel without 
interference) by the R/W criterion indicates, in particular, the following: 


a) As telemetering precision increases, so does the efficiency of coded transmission, but for precision values 
of the order of Na = 100 and less, PCM may be less efficient than other forms of modulation; 


* This condition is not always realized in practice, especially in industrial telemetering. However, the given analysis 
is of value only when the question of transmission efficiency is of interest and is isolated from such factors lowering 
the efficiency as excess bandwidth of the filters in the communication channel. 
** This means for the ChTI-1, for example, a change from an inertial output device with T = 3 seconds to a device 
T = 1 secorid, provided, of course, that the receiver of the information—a man or an automatic device —is in a posi- 
tion to provide such a response time. 

For example, if PWM reception in one cycle is affected by a transformation of pulse duration into amplitude by 
means of linear scanning in the receiver, it is possible to make a similar transformation of the duration of the period 
(or a part of the period) of the oscillations with PFM into amplitude by means of hyperbolic scanning, which in a 


number of cases is well approximated by exponential scanning. 
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b) For these precision values (Na €100), narrow-band telemetering (that is, telemetering based on the trans- 
mission of pulses with gradual fronts by PWM (PTM) makes it possible to attain the highest transmission efficiency. 


APPENDIX 


Technical Data for Several Telemetering Devices 











Name of organ- af. 

izationa LN Ww, N "pitts a Ww = Ww 

tials of device |T-sec| ™ | =*-gc} 5 |" 4) _ : sec [Mt = ap tts) _ _ oe Na 
n+8’ 

VNIIE; TK-2 026 6 | 30 | 2] 64) 0.033 1 16 0.375 

TsNIIKA * 0.25 |8} 40 | 2 |100} 0.025 i 20 | 0.33 
































* In the TsNIIKA system a binary decimal code is used; out of n = 8 code pulses, 4 
transmit of the tens in binary code, and the other 4 transmit the units. Thus, here, 
N= 100, while log N ¥ n. 





































































































; -| a| & 
Name of organization " re fl rT. |w riaae. S & be 
and initials of device |gec ‘lee | NA lem] Na’| S] of toe Na 
séc | sec t 4 -—TW- 
sec > Fe (W = W,) 
Bristol 16 14 | 1.8) 16) 50; 1 0.28 1 56 | 0.1 
IAT AN SSSR; VST -1 2 | 1.6)15 4|50; 2 0.032 1 | 120 | 0.047 
~e 
| Name of organ- » we i a 
c jge |; / |sec gc = 7 
IAT AN SSSR 
ChTI-1 5/15; 3) 3 | 50; 15 0.2 0.04 90 0.062 
TsLEM; ChIS-1 1} 10/10) 5 | 50 5 0.18 |0.18 100 0.086 
TsNIIKA*; ChI-1 | 4} 20}; 5| 3 | 50/ 12 0.32 | 0.08 120 0.047 
TsNIEL-Elektropul’t 
TNCh** 27 | 44 |1.63) 3 | 50) 8! 0.34 |0.012| 260 0.02 
* See the remark’ on page 1081. 
Name of organi- * s T { R log 
zation and initials |7max ip E e i*t= ep Wo= 37" ae a max. 
of device sec | sec oor sec gc 
IAT AN SSSR 
PWM -AR 19 2 50 10 0.2 5 0.56 
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THE STARTING-UP TIME AND ITS EFFECT 
ON THE CHARACTERISTICS OF GYRO MOTORS 
WITH HYSTERESIS 


N. Z. Mastyaev and I, N. Orlov 


(Moscow) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 9, 
pp. 1220-1228, September, 1961 

Original article submitted January 21, 1961 


An analytical expression for determining the starting-up time of hysteresis gyro motors, which 
yields results that are in good agreement with experimental data, is derived. The maximum 
calculated power of a hysteresis gyro motor which would secure the assigned starting -up time 
is derived. An analysis of the starting -up time effect on the characteristics of hysteresis gyro 
motors, e.g., power consumption, efficiency, cos y, and heating, is presented. 


Electrically driven gyroscope gained widespread use as position data transmitters and as differentiating and 
integrating elements of automatic systems. 


The gyroscope starting-up time often determines the time when the entire device is ready for operation, and, 
therefore, its magnitude is in many cases critical in estimating gyro motors. In connection with this, the possibility 


of calculating the starting -up time with sufficient accuracy in designing gyro motors is an important practical prob- 
lem. 


In designing gyro motors, it is even more important to determine the maximum calculated motor power which 
would secure the assigned starting-up time. 


The present article presents the solution of the stated problem with respect to synchronous hysteresis gyro 
motors, which are recently being widely used for driving gyroscope rotors in highly accurate systems, At the same 
time, this article presents a general analysis of the starting-up time effect on the characteristics of hysteresis gyro 
motors, €.g., power consumption, efficiency, cos y, and heating. 


1, Mechanical Characteristic and the Starting-Up Time of Hysteresis Gyro Motors 

The mechanical characteristic w = f(M) of the majority of actual hysteresis gyro metors has the shape of the 
curve a in Fig. 1. The maximum electromagnetic moment of a hysteresis motor usually acts at the starting (w = 0, 
s= 1), Even in the case of motors where the effect of eddy currents in the rotor material is negligibly small, the 
maximum synchronous moment My, is smaller than the starting (short-circuit) torque Mgc. 


For well-designed hysteresis motors, the mechanical characteristic coefficient c,,, which is equal to the ratio 
of the starting torque M,,, to the maximum synchronous moment My, usually lies within 1.2 to 1.8, 


The “useful” load of a gyro motor is the moment Mg of mechanical losses—losses in its own bearings and of 
losses in friction between the rotor and the surrounding medium, The gyro motor's operating point will be the point 
of intersection between the motor's mechanical characteristic w =f (M,y,) and the drag (mechanical loss) torque 
curve Mg =f (w) (Fig. 1). 


The kr = Msc/Menp ratio, where Men is the rated load torque, is called the starting torque factor; it deter- 
mines the maximum motor torque for a certain given rated torque. The motor's overload capacity factor kg = 
= Mms/Men, in turn, determines the degree of loading and the motors overload capacity under synchronous con- 
ditions. 





The mechanical characteristic coefficient c,,, the overload capacity factor kg, and the starting torque factor 
km = Cmko are the most important indices of the mechanical characteristic of a hysteresis gyro motor and they 
determine its starting-up time t,, all other conditions being equal. 
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The starting -up process is described by the equation 


d 
M_ {o) — Mg(o) =J—, (1) 
where My,9(w) is the motor's torque in dependence on angular velocity (g-cm), Mg(w) is the load (drag) torque of 
the motor (g-cm), which is a function of the angular velocity, J is the moment of inertia of the motor's rotating parts 
(g-cm sec’), and w is the angular velocity of the motor (rad/sec), 


As is known, the area enclosed by the mechanical characteristic and the drag torque curve of a gyro motor 
(the shaded area in Fig. 1) is used as the measure of the starting -up time t, for certain given J and wg values. 














Fig. 1, Mechanical characteristics w =f (M,,,) and the drag torque 
curves Mg = f (w) of a hysteresis gyro motor. 


The gyro motor’s starting-up time can be determined by solving Eq. (1) with respect to time t: 


ed 
J (2) 
t= da. 
y >: a 





By expressing the mechanical characteristic M,,.(w) and the load curve Mg(w) analytically, we obtain an ex- 
pression for the starting -up time ts of the motor for certain given J, wg, Men, km, and ko values. 


The configuration of the actual mechanical characteristic of a hysteresis motor makes it possible to replace it 
by two straight-line sections (curve b in Fig. 1) and to approximate it by the equations: 
napa teWE] mosey 
0M Ck Min for 9 = @, 


The drag torque curve Mg = f(w) is generally a complex function of the angular velocity, and it has a constant 
component Mqo, which is independent of the angular velocity and which is determined by static friction in the 


bearings as well as a component which is approximately proportional to w'*® and which characterizes friction against 
the surrounding medium. However, for the sake of simplicity, we shall represent the drag torque curve in the shape 
of the straight line d (Fig. 1) by using the equation 


M, (@) = Mz, oye 0<o<q,. (4) 
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Comparison between the Theoretical and Experimental Data Used for Deter - 
mining the Starting-up Time of Hysteresis Gyro Motors 

















‘ Gyro motor specimens 
Quantities | § 2 
ssi! 2 8 f 5 6 1 
K Bec |27 |28 |49.5 |63 230 | — | 4000; 
M2n g-cm 3.0 | 2.6 | 4.3 | 3.0 3.4 | 14 110 
ko - 1.7 | 2.2 | 1.65] 2.4 2.44) 2.08 1.1 
Cm _ 1.6 | 1.5 | 1.2 | 41.53] 4.57) 1.2 1.75 
kin _ 2.7 | 3.3 | 2.0 | 3.67] 3.83) 2.5 1.9 
ts ,(experiment)} sec 5.0 | 4.9 | 9.0 | 7.0 | 26.9 | 60 1080 
t; [calc, b sec 6.1 | 5.17) 9.6 | 8.8 | 27.5 | 76 1390 
theans of (6) 
ts [calc. b sec 6.1 | 5.4 |10 8.7 | 27.2 | 74 1210 
means of (12)] 
t"; /ts oan 1.22) 1.05) 41.41] 1.25) 1.04) 1.23 1.12 



































It is readily seen that, in approximating Myo(w) and Mg(w) by Eqs, (3) and (4), one can expect higher theo- 
retical values for the starting -up time t, in comparison with experimental data for the majority of gyro motors. 


By using Eqs. (2)-(4), we obtain: 
ed 





‘i J dw 
$"M ° (5) 
2 out ae suaes 
Domo * ( k +1) ed 
Integration yields the equation 
K 2,3 hk 
t ae) oe a . > A ee —n_. 6 
P Ms n*¥ mk ot 1 Bio 1 ” 


where K = Jwg(g-cm sec) is the gyroscope's angular momentum, and t, is given in seconds, 


The obtained expression (6) makes it possible to determine the starting-up time t, for a given hysteresis gyro 
motor with known K, Men, ko, and km values, and also to analyze the effect of various gyro motor parameters on 
the starting -up time. 


The table provides experimental data for a number of hysteresis gyro motor specimens and the starting -up 
time ts values, which were calculated by using Eq. (6). The discrepancy between theoretical and experimental data 
does not exceed 25%, Although this cannot be considered as a sufficiently good agreement, it is entirely acceptable 
for approximate estimates, Here, one should mention the difficulties encountered in determining the mechanical 
characteristic of gyro motors, especially with respect to the maximum synchronous moment Myns, which additionally 
reduces the accuracy of calculations. 


For equal K and Men values, the starting-up time is the shorter, the larger the starting torque factor kyy, while, 
for a given kyp, it is the shorter, the larger the motor's overload capacity factor kg (or the smaller the c,, coefficient). 
The dependence of tsMgn/K on kp for different cp, values, which was calculated by using (6) and is shown in Fig. 2, 
confirms the above statement. The minimum starting-up time t, for a given kp, will be obtained for ko = ky, (Cm = 1) 
(curve b in Fig, 1). 


According to (6), for certain given k, and k,,, the starting-up time t, increases with an increase in the angular 
momentum K, This increment is the larger, the smaller the gyro motor's rated torque Mgn, since this signifies a de- 


crease in the excess moment Mp,9—Mg = Men [ky —(&m —Ko + 1) at By expressing K and Mop in terms of the dimen- 
sions and parameters of the gyro motor [1], their effect on the starting-up time can readily be determined. 
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Actually, the expression for the angular momentum for a cylindrical rotor (Fig. 3) is given by 
L D,\4 ; 
K = Jog= 1.047+40-%D* 5° | 1— (5!) | n@-em sec), © 


while the following equation can be used for calculating the drag torque in a certain given range of Do and n: 


P,.,-408 oe = 
Moy =f 5,- = 0.97 (aapo V pp’) Vib Df (4 + 5 59) 10° (g-cm) 8) 


Here, Do, Dy, and Lo are the dimensions given in Fig. 3 in centimeters, y is the specific weight of the rotor 
material in g/cm®, n is the gyro motor's angular velocity in rpm, Bg is the ventilation loss factor, p is the working 
medium density in g sec’/cm*, ut" is the viscosity coefficient of the working medium, ag = 1 + ( APpear/ APyent) 
is a coefficient which takes into account the portion of losses in the bearings in the over-all drag torque, APpear is 


the friction loss in the bearings in watts, and APyen, is the loss resulting from friction in the medium (ventilation 
loss) in watts. 
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Fig. 2. Relative starting -up time t,Mpn/K of a hysteresis gyro 
motor in dependence on k,, for different c,, values, The points 
denote the results of calculation according to Eq. (6), and the 
triangle denotes the results obtained by using Eq, (12). 


By means of (7) and (8), we readily obtain: 


K__ 4.08.107p, —»/ Po_ 1 1 Oy! Pol 


jy (9) 
Men OFF Ty IDoVn abeVon) 





As follows from expressions (6) and (9), for certain given k, and kp, the starting-up time is the shorter, the 
smaller the gyroscope dimension Do, the greater the angular velocity n, the larger the ratio Dj/Do of the rotor 
diameters, the smaller the Lo/Do ratio and the specific weight of the rotor material, and the greater the density of 
the medium. 


The conclusions drawn here apply best in the case of gyroscopes with small bearing losses (ay » 1), when the 
dependence of the drag torque on angular velocity and the gyroscope dimensions can be represented by Eq. (8) with 


1094 





9) 


with 





sufficient justification. If the bearing loss is considerable (a9 > 1), the drag torque is largely determined by the 
rotor's weight, which materially changes the conclusions following from expression (9), The latter applies especial- 
ly to gyroscopes with a small angular momentum and to vacuum gyroscopes, 


2. Theoretical Power of Hysteresis Gyro Motors 
In designing gyro motors, a problem which is opposite to that considered above must be solved: We must find 
the power necessary for securing the assigned starting-up time t,, which is to be used as the basis in designing the 
motor. The practical problem can be reduced, in the first place, to the determina- 
— 6 —> tion of the value of coefficient k,, (or ko) that is necessary for securing the required 
starting -up time t, for certain given K and Men values and for a given mechanical 


IY characteristic coefficient (cm = 1.2-1.8), The value of ko can be determined by 


graphically solving Eq. (6), which can readily be transformed in to the following form: 











tsMon 


‘ k 10 
=f 53K Mt (eg 1) Rol = 18 in FG * tes 








s 
| However, the use of this expression is difficult, especially for general analysis 
purposes, and it is desirable to derive a simpler relationship between k,, and the 
motor's starting-up time t,. It would be more correct to eliminate the logarithmic 
dependence (10) by using expansion in a power series. However, for actual values of 
the starting torque factor (ky, = 1.5-5.0), the series must contain not less than two or 
three terms, which results in a rather complex equation for k,,. It proved to be sim- 
pler to use the approximation of the tsMgn/ K = f (Km, Cr) dependence in the region 
of actual values of Crp (cm = 1.2-1.8) and ky, (km = 1.5-5.0) which is shown in Fig. 2. The approximating equation 
is given by 








\ 





Fig. 3 


0.725¢ 
ky 1+ " (11) 


In this case, the starting-up time t, is determined by the equation 


pe K_ 0.725¢ (12) 
8 MF 1 
Expression (12) was used for calculating the starting-up time tj values given in the table and in Fig. 2, The 
agreement between these values and the results of calculations where (6) is used justifies the use of Eqs. (11) and 


(12) for practical calculations and general analysis, 


Hysteresis motors use the maximum electromagnetic power during the starting (short-circuit) period, If the 
electromagnetic power under the rated conditions is equal to 





er Ns Na ® 


the electromagnetic power at the starting point will be found by using the equation 


1.03nM,k_-10-* 


Ne (w ), (14) 





P 
Peyw= kn" = 


where 2 is the “output” efficiency, which determines the portion of surface losses APs, in the material of the 
rotor that are due to higher harmonics: 
P AP, 
as gt DS ap | oe og 
= eer 


P 


et er (15) 
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The maximum electromagnetic power Pew represents the calculated power of a hysteresis motor. It should be 
emphasized that precisely this point in the starting of hysteresis motors can be most accurately calculated by using 
the energy method. By using (8), we can readily find Mzn with respect to the assigned angular momentum K, the 
starting -up time t,, and the given rotor dimensions, while the starting torque factor k,, can be found by using (11) 
with respect to the assigned mechanical characteristic coefficient cy), after which the power Pew (14) for which the 
gyro motor is designed can readily be determined. 


By using Eqs. (11) and (14), the calculated power Pew of a hysteresis gyro motor can readily be expressed direct- 
ly in terms of the starting-up time t;: 


P,,, = 1,03-10- M,, *(' +a) a (16) 
—s_‘p 
K 


The shorter the starting -up tivne t,, the larger the calculated power P,,, of the gyro motor, all other con- 
ditions being equal. 


By using (7), (8), and (16), the calculated power Pew can be related to the parameters and to the ratio of the 
gyto motor dimensions for a certain given angular momentum K and the starting -up time t,: 








ws.ven) yr _'*52 
> n 
BP, = 0,956 +108 K See PP <5 pear * 
ty 
Do D, 
$c 1814. cae PoDo! [1 (5"] (w) (17) 
(a8, V pp’) ten Ly D, it 
Cee - ae 


The shorter the starting-up time t,, the larger the second term in the large brackets in expression (17) (the 
larger the necessary excess moment of the motor), and the smaller the effect of the surrounding medium parameters 
(Bo V py’), the dimensions (Do,L/Dj, Dj / Do) and the specific weight of the gyroscope rotor on the calculated 
power, since the drag torque portion in the calculated torque of the gyro motor decreases. The remarks made in the 
analysis of expression (9) also apply to the above conclusions. 


3. Effect of the Starting-up Time ts on the Characteristics of Hysteresis Gyro Motors 

The found relationship between the calculated power P,,, of hysteresis gyro motors and the starting -up time ts 
makes it possible to design gyro motors with an assigned ts. Since the calculated power here depends on the starting- 
up time, the latter also determines to a large extent the electrical gyro motor characteristics (the current 1,, the 
power consumption P;, cos y, and the efficiency 7 ) and its heating. 


The effect of the starting-up time t, on the gyro motor characteristics can be conveniently determined by com- 
paring the characteristics of an actual motor with a given kyp (that secures the assigned ts value) with an idealized 
motor which has no torque margin (ko = ky, = 1). The mechanical characteristic of such a motor is given by the 
curve e in Fig. 1. In this case, as can readily be seen from (6), the starting-up time is t, = (the angular velocity 
rise follows an exponential curve). Comparison with such an idealized motor makes it possible qualitatively to 
estimate the effect of any finite value ts #« on the hydraulic motor's characteristics.® 





* We assume that both motors to be compared have equal angular momentums K, equal absolute dimensions and 
their ratios (which are determined by the value of K), equal stator winding specifications, and equal values of the 
hysteresis angle sin 79 in the magnetic materials of the rotor's active part, while it is assumed that these materials 
are used in the optimum manner. The latter meansthat the efficiency n and cos ¢g values are equal at the calculated 


starting point. The motors to be compared differ with respect to the power at the calculated point (lqq =2 and “Ey, 


the air-gap induction B;, the “hardness” of the rotor's magnetic materials, and the operating voltages U, acting in 
the stator windings. 
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By using the fact that, for an optimally designed gyro motor, the required over-all magnetizing force F, in 
the stator winding is determined by the air-gap induction value Bg, which is proportional toVPey , all other con- 
ditions being equal, we readily obtain:® 


fh so Ft B, _ VP 
Tw — Fo By Ve Vat (18) 


By taking into account (11), we obtain the relationship for the stator's phase current: 


Ih = I i+ ce (19) 


a 


The shorter the starting-up time t,, the larger the hysteresis gyro motor current, all other conditions being 
equal, 


The power P49 used by an idealized hysteresis motor can be expressed by the following equation: 
Pro = Py, + ZAPo = Py + kalo + keBi,, (20) 


where LAP = AP, + (APs, + APgyr) is the total loss in the motor, AP, = kylyis the loss in copper (4Pst + APsur) = eB 
is the loss in steel (which takes into account the surface losses in the rotor), and ky and ke are proportionality factors. 


For an actual hysteresis motor under the rated conditions (neglecting small changes in the motor current that 
are due to load variations), the power used is equal to 


P, = P,,+ ZAP = Py, + kyl + kyBS (20a) 


or, by taking into account (18), 


Pr = Pon + hak Tig+ hak Bi, = Pen + k EAP =Pralt + = hea] (21) 


P ZAP 
Here, o> — = ean Pa is the relative efficiency of the idealized gyro motor at the calculated point. (Ac- 


cording to the footnote (* ) on page 1096, the efficiency m9 and cos ¢ for an idealized gyro motor are equal to the 
efficiency n,, and cos gy of an actual gyro motor at the starting point.) 


By using (11), we obtain: 





Pim Pa [14 Ge (t+ Say) (22) 
2n 
or 
Pt = 4 + 0.725 (1 — 19) i. 23) 
Py, t. 2n 


Thus, with an increase in the starting-up time ts, the power P; used by a hysteresis gyro motor under rated 
(operating) conditions increases. 


Correspondingly, by using (23), (18), and (19), we can obtain the expressions for the duty values of the efficiency 
” and cos ¢ for an actual gyro motor: 


* Hereafter, the index "0" will pertain to the idealized gyro motor, where ky = 1, 
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XT Pr | 
No 








> 725 — 4 
1 me are To) ¢ (24) 
s 2n 
and 
4 + 0.725 (14 — wt 
cos@ Pi IU _ Pr (48) aiid an'p : (25) 


cos@ Poli Ui Pr tT) pee ‘aft 
s n 


It can readily be seen that a decrease in the starting -up time t, directly leads to worse energy indices (n and 
cos ¢) of hysteresis gyro motors under rated (operating) conditions. 


Finally, since all of the electric energy used by the gyro motor is transformed into heat energy, which causes 
the heating of the gyro motor, it is readily seen that the overheating temperature 1° of the gyro motor increases 
with an increase in the starting -up time tgs: 


a. ae es 6 
Px 4 + 0.725 (4 No) iM," (26) 


Thus, the presented qualitative analysis of the effect of the starting-up time t, on the characteristics of hy - 
steresis gyro motors makes it possible to draw definite conclusions, which indicate that, with a decrease in the starting- 
up time t,, the duty (operating) values of the current 1,, the power consumption P,, and the overheating temperature 
T° increase, while the efficiency n and cos ¢ decrease, all other conditions being equal. 
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1098 








ng- 





ANALYSIS OF A CIRCUIT FOR CONTROLLING A DC 
MOTOR BY MEANS OF A BRIDGE REVERSIBLE SEMI- 
CONDUCTOR AMPLIFIER 


G. B. Eliasberg 


(Leningrad) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 9, 
pp. 1229-1234, September, 1961 

Original article submitted March 9, 1961 


There is considered the control of a dc motor by means of a reversible amplifier using transistors 
which operate as switches. Mechanical characteristics and a motor efficiency are determined, 
Recommendations are given for the choice of a switching pulse frequency. 


In recent years there has been intensified interest in the possibilities of controlling dc motors by means of 
transistors operating as switches. There are several papers [1-3] devoted to this question in which the basic prin- 
ciples of this control method are analyzed, and concrete amplifier circuits are presented. 


In [2] there was given a detailed analysis of a nonreversible control circuit for a dc motor. For the purposes 
of automatic control, however, there is greater interest in circuits which have the possibility of reversing the motor, 
and also of its operation both as a motor and in a reverse connection as a brake with regeneration. Practically, the 
most advantageous circuit, which insures meeting the indicated requirements, is the reversible bridge circuit [1, 3], 
The present article is devoted to an anslysis of the operations of this circuit. At the same time there is substantial 
interest in the questions of determining the mechanical characteristics and the efficiency of the motor, and also in 
the choice of the working frequency. The principle attention in the article is paid expressly to these. questions, 


1. Basic Operating Principles of the Bridge Reversible Circuit 

In Fig. 1 there is presented a bridge circuit for controlling a dc motor by means of transistors. The circuit is 
controlled with rectangular voltage pulses having a period T and an amplitude Vp (Fig. 2a). This voltage is fed to 
the transistor bases in such a way that for the period of time 
from t= 0 to t= T, transistors Q, and Q, are conducting and 
transistors Q, and Q, are cutoff, while in the period of time t = 
= T, tot = T transistors Q, and Q, are conducting and transistors 
Q, and Q, are cutoff. By altering the ratio of T,/T it is possible 
to regulate the speed of rotation and to reverse the motor. For 
T,/T = 1/2 the average armature current I,, = 0, and the motor 
does not develop a torque. Diodes D1-D4 are provided to en- 
sure uninterrupted commutation (see [3]), The apparatus which 
shapes the control signal can have a number of forms which are 
described in the literature (see, for example, [3]). Therefore 
its operation will not be considered here. 














Fig. 1 


We will adopt the following assumptions: 
a) The current through a transistor which is cutoff, and the switching time are equal to zero; 
b) the change in rotational speed over a period T is insignificant and can be neglected; 


c) the voltage drop of a transistor when in the saturation region, and also the forward voltage drop of a diode 
are equal to zero. 


For these conditions the parameters of the armature circuit of a motor (its ac resistance and inductance) re- 
main constant during the period T, i.e., they do not depend on the direction and path of the current in the circuit or 
on the transistor condition. 
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The equivalent circuit of Fig. 3 can then be substituted for Fig. 1. In this circuit the source of equivalent 
voltage Vg provides voltage pulses of rectangular shape with an amplitude equal to the supply voltage V,,a period T, 
& positive portion of the period of duration T,, and a negative T-T, (Fig. 2b). The counter emf is designated by E 
and the armature current by i. 
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The arrows indicate the positive directions of the currents and voltages in the circuit of Fig. 3. 
For the time interval from t = 0 to t= T, the correct equation is 


Var: ™ di (1) 
pos Rese nl — 
r =i+TH 


where r is the resistance of the armature circuit, r is the time constant of this circuit. 
For the interval from t = T, to t = T the correct equation is 
Vv : 
s+E&_s,; di 
a See =t+ts. (2) 


Solving these equations and taking into consideration that i(0) = i(T) we get: 











for the interval from t = 0 tot= T, Vig rig! whe (3) 
i= (i{—e *)+iQe *, 
for the interval from t = T,; tot = T (4) 
where Feet t (i—e 4 i (71) e 
T—T, T 
ee mJ 5) 
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V.(e —d2e + 1) B (6) 
< 


r(ii—e ) 


It is obvious that i(0) is the minimum value of the armature current during the period T, and i(T;) is its maxi- 
mum value. In the time intervals between these values the armature current changes according to an exponential 
law. A characteristic curve is presented in Fig. 4. 


The values of i(0) and i(T,), each taken separately, can assume either positive or negative values depending 
on the circuit parameters and the ratio T,/T. 


2. Mechanical Characteristics and Motor Efficiency 








Expressions (3) and (4) are cumbersome and inconvenient for further study. Therefore it is expedient to make 
some simplifications which permit the calculations to be simplified while producing little effect on the final results. 


Let us introduce the voltage V, in the form of a Fourier Series: 


fo 2) 
. 2ukt 2ukt 
Ve =Vayt 2 (Axsin a + By cos), (7) 
where 
~Z +7, q 
Vav =| | vy ae \ v, a}, (8) 
T T 
-> —-=+T, 
—-—-+T, t 
Ay = = { Vs; sin cil dt — \ v, sin dt }. (9) 
T T 
~9 “— 
—Tit, , 
2 Qmnke , 2nkt 
By, = a ( Vs cos a aut — \ Vv, COs dt |. (10) 
T T 
_ ry =o +h 


* 


After appropriate transformations we get 


T 
Vay =v, (29-1), (11) 
2V T 
Ax = (— 1)" —8.(1 — cos), (12) 
2v 2nkT (13) 


B, = (— 1)* —#- sin =. 
The average value of the armature current is, in a similar way 


T; 
g2eog 2) Ugo» Be (27—-1)—8 (14) 
AV r oe, r , 
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Expression (14) is analogous to the equation for the armature circuit of a motor being controlled by a dc volt- 
age V if we set V= V,(2T;/T~-1), Therefore the mechanical characteristics of the motor are identical for both con- 
trol methods. 


When T,; = T/2, Vay = 0 and the motor is operating ina regenerative braking condition with increased power 
losses on account of the pulsations. 


From the power point of view it ie of interest to determine the effective value of armature current Ig and the 
motor efficiency K, = Iay/Iggf. 

Designating by 1, the effective value of the kth harmonic of the current i, by V;, the amplitude of the kth 
harmonic of the voltage Vp, by Za, the impedance of the armature circuit for the kth harmonic, we get 
U, 


(15) 
V2Za, 





Ty = 


Ve = Vat Bh a0 


Zak =r V1 + (AS (17) 





Taking into account (12) and (13) we get from (16), after some intermediate transformations 


4v, 
Vy = > sin akg, (18) 
where 


T 
1=7- 


It is shown below that in order to get an acceptable value for efficiency, it is necessary to select a period for 
the control volt which is of the same order as the time constant, or even less. Taking this into account expression 
(17) can be simplified: 


2uktr 
Zak > - (19) 


Let us find the value of the ratio between the current of the kth harmonic I, and the current of the first har- 
monic I,. Taking into account (15), (18), and (19) we get 


Ty VxZA1 _ 14 sinnkg (20) 


SSC 


Study of expression (20) shows that for any values 0 < q < 1 the following relationship is true 


Ty U1 (21) 
St" 
The ratio of heat losses, determined by their corresponding current harmonics, satisfy the inequality 
P,: 4 (22) 
P, => r < r . 


Expression (22) shows that when determining motor efficiency it is possible to neglect all current harmonics 
except the first without much error. 


We obtain the value of current I, from (15), (18), and (19): 
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2) 


I V2VsT (23) 
1 = “ Sin nq. 


Knowing the value of the average armature current I, and having found from (23) the value of the first har- 
monic current, the efficiency can be found from the formula 


K, = P. om (24) 


Vigy tit’ 


3. Selection of the Control Voltage Frequency 








Expression (23) shows that for any given value of q the current I, does not depend on either the load or the 
rotational speed. For q= 0.5, i,e., in the dynamic braking condition, the value of I; is a maximum. Moreover, the 
value of I, is inversely proportional to the frequency of the control voltage f= 1/T, Therefore, for the best motor 
efficiency it is desirable to increase the frequency. However, as is well known, this brings about an increase of the 
switching losses in the transistors and additional heating in them. Therefore the frequency should not be increased 
above the value which is required to give an acceptable motor efficiency. 


The determination of the required value of f must be directed towards the operating condition of the motor for 
which the effective value of the armature current is a maximum. For example, in servo system motors such an 
operating condition is represented when the signal at the input of the system amplifier is somewhat less than the sig- 
nal which starts it in motion. In this instance the armature current is close to the starting current but the speed is 
equal to zero. Such a condition, generally speaking, can exist for a long period. We will assume that the starting 
current I gr isequal to the nominal armature current Ijoyy. Generally speaking, in servo systems Igy < Iyoy: 48 
arule, However, in practice cases are encountered when Ist ~ Iyon. Obviously such a case corresponds to the 
most difficult operating condition of a motor, and therefore it is expedient to consider it expressly as an example. 


Let us denote by Kg7 a multiplier for the starting current 
Vv 


ae ode 5 
Ksr =77*- (25) 
NOM 


From (14) and (25) we have for the above-mentioned condition 





T, 
2 a —!1 
a v, _v, 29+) 
ST 
Whence we find 
inal 
{eo ST 


For a de motor the value of Kgr is generally not less than 5-10, Therefore the value of q as determined from 
expression (27) is close to 0.5. 


Hence for the operating condition being considered the first harmonic of the armature current 1, must have a 
value close to the maximum at the given frequency. Since lay * Iyoyy; it is obvious that this condition is the 
most difficult of all those possible. 


Supposing, for example, that the allowable value of efficiency K, is equal to 0.9, we find from (24) 


a (28) 
i, = AV Vi — Ki 0.48 1,,. 


Ke 
Substituting (27) and (28) in (23) and taking into account (26) we get 
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1 
7 § 4- 
T (29) 
0.48 2 = V2v,1 sin x al 
r 





Kst n*tr 2 
whence 
ee 0.48x? 
& tar: 1 . 
*(1+57) (0) 
KsT V2sin 5 





Expression (30) determines the value of the period T, or the frequency f= 1/T, which is required in order to 
obtain an acceptable motor efficiency. 


To increase the frequency above this value is not expedient because it brings about an increase in the transistor 
losses. It should be noted that the value of T which is found from expression (30) is close to the value of the time 
constant. For example, when Ky = 5, we get T = 0,77 from (30). 


It must be emphasized that expression (30) determines the frequency only for the particular case, given as an 
example, when Igy * INoM: Kst = 5, and Ke = 0.9, In the general case, knowing the armature parameters (r 
and r), and also the current lay for the worst condition, and assuming a permissible value for K,, these values must 
be substituted in expressions (23) and (24) and,by a common solution obtained of the equations, the required value of 
f is then found. 


SUMMARY 
1, The mechanical characteristics of a motor which is controlled by means of a bridge amplifier using tran- 
sistors does not differ in principle from the mechanical characteristics of a motor which is controlled by means of a 
change in the voltage supplied to the armature. 


2. When choosing the frequency of the control pulses the value of motor efficiency found from expressions (23) 
and (24) must come within acceptable limits. The frequency must not be increased above the value determined by 
these expressions, 
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It is shown that there is a possibility for carrying out practically inertialess temperature measure - 
ments of gases and liquids in the following way: the temperaturesare taken while the fluids are 
flowing through a test chamber which is provided at the inlet and outlet with flow impedances 
having characteristics (flow versus pressure differential) that vary with temperature in a different 
way. An analysis of the flow characteristics through such a chamber is carried out. Methods for 
accounting for and compensation of the errors are given which are caused by heat exchange of the 
fluids with the walls of the test chamber. 


Ordinary pneumatic and hydraulic temperature sensing devices which have widespread technical applications, 
specifically in automatic control systems, are suitable only for measurement of stationary or slowly changing tem - 
peratures. The explanation for this is that the basic principle of their operation is functionally connected with heat 
transfer from the gas or liquid to the container walls and through the container walls. 


As a characteristic example, a schematic sketch of one of the pneumatic temperature measuring devices of 
this type [1, 2] is reproduced in Fig. la. The gas whose temperature Ty is to be measured is flown through a measur- 
ing tube in which are mounted the flow impedances 1, 2, and 3 (orifices, baffles, valves, filters, etc.) and the ex- 
changer 4. The temperature of the gas Ty being tested is judged by the difference in pressure 5 p, before and after 
flow impedance 1. This proves to be possible, if the mass flow of gas in the measuring tube is kept constant. This 
is made sure of by maintaining the pressure differential 5p, across impedance 2 constant by adjustment of the open 
passage area of impedance 3, and also by keeping the temperature T, of the gas reaching 2 constant; the latter is 
accomplished by cooling of the gas from its temperature T, to the temperature T, by means of the heat exchanger 4. 
Heat exchangers or heat accumulators of one kind or another are unavoidable accessories of other pneumatic and 
hydraulic temperature sensing devices also [3 to 6 and others]. 


Of special interest is the principle of temperature measurement described in [7], even though the form of 
measuring apparatus outlined in that paper is suitable only for measurement of stationary temperatures. A dia- 
gramatic sketch of this apparatus is shown in Fig. 1b. Change of pressure in chamber 1 as a function of temperature 
of the flowing -through gas is here obtained by the use at its ends of flow impedances of two types having radically 
different temperature characteristics: one of these—2—is a flow impedance producing laminar flow, consisting of a 
plate of porous material; the other—3—is a flow impedance creating turbulent flow. Ahead of impedance 2, a filter 
4 is installed. To maintain the walls of chamber 1 at a temperature equal to that of the gas flowing through the 
chamber, the gas is passed at the same temperature through jacket 5 surrounding chamber 1, where it is continuously 
flowing over the outside of the chambers walls. Inasmuch as the characteristics of this device too are influenced by 
the temperature of the walls, the process taking place is just as much subject to inertia as in other known pneumatic 
temperature sensing devices. 


It is shown below that the system with fwo flow impedances of different type can serve as a basis for develop- 
ment of methods for practically inertialess measurements of momentary values of the temperature of gases and 
liquids in nonstationary processes, when their composition is known and fixed. 
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Fig. 1 


1. Analysis of the Basic Characteristics of Pneumatic and Hydraulic Chambers 
Terminated by Flow Impedances of Different Type when Used as Temperature Sensors 
1. Initial premises. Preliminary considerations regarding the use of chambers of the described type for prac- 
tically inertialess measurement of nonstationary temperatures. Let us examine a chamber which is provided at its 
inlet and its outlet with flow impedances of different type. Let, for instance, the impedance at the inlet be of a 
type causing turbulence, and the one at the outlet of a type providing laminar flow (Fig. 2). The mass flows of a 
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gas or a liquid through these two impedances have different temperature dependences. To illustrate the prevailing 
conditions, in Fig. 3a characteristic curves are given which show the change with temperature of the mass of air 
passed per second through a turbulence type impedance and through a laminar type impedance. The ordinates of 
these curves give the ratios of the mass throughput at the corresponding temperatures on the abscissa to the mass- 
throughput at 15°C under the same conditions. In Fig. 3b are shown experimental data obtained when water was 
passed through the two different flow impedances; however, here, the throughputs per second are plotted with refer - 
ence to the throughput at 8°C, instead of 15°C as in the case of air. In both Fig. 3a and Fig. 3b, the curves marked 
1 were obtained with turbulent flow, and those marked 2 with laminar flow. 
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When a gas or a liquid flows through a chamber provided at its ends with two such different impedances (Fig. 2), 
then the pressure p; in the chamber —with the pressures py and p, remaining constant —varies as a function of the tem- 
perature flowing through the chamber medium. And here heat exchange through the walls is no more a neces- 
sary element of the temperature measurement process. On the basis of the above reasoning, in [7] the conclusion is 
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reached that the use of a chamber with different kinds of impedances at the ends is advantageous for increasing 
precision in determinations of stationary temperatures under conditions of small gas throughputs. 


If it is possible under actual conditions to achieve the chamber characteristics with respect to absence of heat 
exchange with the walls, indicated in this section, or to eliminate the effect of heat exchange altogether by other 
means (see section II), then it is also possible to measure temperatures of gases and liquids in nonstationary, fast 
going processes by means of a chamber terminated by two flow impedances of different type. Indeed, when the effect 
of heat transfer upon the function py= ¢(T) is eliminated, then the time constant of the temperature sensor is deter- 
mined only by the velocity of the liquids flow through the orifices of the flow impedance and through the chamber. 

In the case of gases, the time constant depends also upon the speed with which the pressure in the chamber is built 
up. However, when the dimensions of the chamber are sufficiently small, the time constant —with all these processes 
present —still is of the order of a hundredth's or a thousandth's part of second [8]. 


2. Calculation and analysis of the relationship p, = ¢(T) for gases. Suppose that in the chamber shown in 
Fig. 2 the flow impedances 1 and 2 are of such construction and the pressures pp and p» are so chosen, that the first 
of these impedances is of the turbulence producing kind, the second-—of the laminar flow type. 





We shall use formulas (1, 10), (1, 11), (I, 26) and (1, 83) of paper [8]. Using these formulas, we consider the 
processes in the flow impedances and in the chamber to be quasistationary. We must remark that formula (I, 26) 
is valid for isothermic conditions, while, as the problem is set up here, the flow should be considered of an adiabatic 
nature. However, there will be no serious error incurred by using formula (1, 26), because, when the pressure drops 


involved are small, the throughputs obtained at equal pressure differentials are practically the same in both kinds of 
flow. 


Using the indicated formulas, we write down the equation expressing equality of the quantity of material passing 
through the two impedances per unit time—twice: once for the normal air temperature (288°K = 15°C), and the second 
time for the temperature of the gas under test T. Designating the pressures in the chamber that correspond to these 
two temperatures respectively by py and py, we obtain from the above two equations the following expression: 


0(2)=0(%) em, ay 


where when the flow through impedance 1 is below the critical value we have: 





Pi\ __ (Pr / Po)* — (P2/ Po)* 
(2 = - 





V (pil po)®!* — (pri po) *t* ' 7 
and when the flow in 1 is above the critical, we have: 
(2 = (Pr! Po)® — (pal Po)® 
Po k—1/ 2 \2/-1) (2") 
k+i CH 
From the above we derive: 
» (7) = (421 + 1,627) VT-10-. (3) 


The curve representing the function »(T), plotted according to this formula, is presented in Fig. 4. However, 
while the temperature T in the above formulas is expressed on the absolute scale, in Fig. 4 the temperature is ex- 
pressed in centigrades. 


Curves giving the values of ¢ (p;/pg) for the entire range of feasible flow conditions through the impedance 1 
and for a number of p)/pp» ratios, as calculated by formulas (2) and (2') with k = 1.4, are presented in Fig. 5a. So- 
lution of equation (1) above with the help of these curves is carried out as follows: First we find, using Fig. 5a, the 
value of 9 (py /Po) which corresponds to the pressure py (ki.own from the setup of the problem), and then we mul- 
tiply it by the values of »(T) corresponding to a number of chosen values of T taken from Fig. 4. From the obtained 
in this manner values of the quantity ¢ (ps/p9) [see equation (1) ] we obtain the corresponding values of p;/p» by 
using a suitable curve of Fig. 5a, and thus, since we know pp, we have found the values of p; corresponding to the 
chosen values of T. 
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When interpolating for values of p,/pp between those for which curves are given in Fig. 5a, the following con- 
sideration should be kept in mind which will make plotting of intermediate curves easier: each of the ¢(p,/po) curves 
intersects the abscissa axis at a point whose distance from the origin is equal to the p,/pp ratio of this very curve. 


Let us examine a numerical example. Let us have for the chamber of Fig. 2 the following parameters: pp» = 1.5 
atm (absolute), pp = 1 atm (absolute), and py = 1,1 atm (absolute), Let us determine the. function p; = 9(T) which 
will be valid when the flow impedances at the ends of the measuring chamber are of the above indicated character. 


Carrying out the computation by the described graphic method we obtain the relationship p, = 9(T) shown as curve 1 
in Fig. 5b. 
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The steepness of the curve p; = 9(T) depends substantially upon the value of pp: when, leaving the values of 
P2 and p, the same as in the preceding example, the value of pp is increased to py = 5 atm, we obtain for p; = 9(T) 
the much steeper curve 3 of Fig. 2. It should be noted that, to achieve this with all other characteristics remaining 
constant, it is necessary to decrease the cross section area of flow impedance 1. 


All given above data pertain to the arrangoment of the flow impedances of diverse type 1 and 2 shown in Fig. 2. 
However, the reverse arrangemcnt of these flow impedance is also feasible, that is, having the laminar impedance 
in the inlet, and the turbulent in the outlet. The formulas for computation and the auxiliary curves for this case can 
be obtained in a manner similar to the one used above in connection with the arrangement shown in Fig. 2. 


3. Calculation of the characteristics of p, = ¢(T) for a liquid. Let us examine as before the chamber shown 
in Fig. 2, Let the flow through impedance 1 be a turbulent one, and let us suppose that for this particular im- 
pedance the local resistances at the inlet and the losses at the outlet are of prevalent importance; let the flow through 
impedance 2 be a laminar one, and in it let the friction losses in the channel be of major importance. 


The throughput of liquid through flow impedances 1 and 2 are computed respectively by the formulas: 


G=e/V2e7Vpo— Pi 








nds i 
G = Far (1 — Pe) 5 


Equating the two throughputs expressed by these formulas, we obtain: 


Vro— Pi _ ndtVy 1 
Pai—Ps 12B0efV2gh - (4) 
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(4) 





In distinction from gases, in dropforming liquids viscosity does not increase with temperature, but decreases, 
and in addition their specific gravity changes less with temperature than that of gases. As an example, in Fig. 6a 
are shown curves expressing the change with temperature of the viscosity » of water (curve 1) and of a motor oil 
(curve 2); (viscosity is here measured in practical metric units of kg-wt-m A sec, translator), On the same figure 
another curve (3) is given, which shows the change of specific gravity with temperature for water, This is taken 
from [9]. According to the data of Fig. 6a, the quantity/ y (y = sp. gr.) for water decreases by less than 2, while 
u is decreased by a factor of 6.5, when the temperature increases from 0 to 100°C, 
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Considering thereforeVy in expression (4) as one of the constants, and taking into account the thus found 
dependence of 1: upon temperature we obtain from (4) the following equation: 


V po — Pips — p-) 
a. et a = p(T). (5) 
V Po— Pu /(Pr. — Pe) vO) 





where the asterisk added to the subscript of p, now indicates the value of py at 0°C (and not at 15°C as in the preceding 
part). 


The curve representing the function y(T) for water is shown as curve 4 in Fig. 6a. 


Equation (5) contains only ratios of pressure differences; therefore, this equation is valid, regardless whether 
all three pressures pp, py, and p, are measured in absolute units, or as the excess values above atmospheric pressure, 
or in any other units. 


Let us examine the following example: through a system of flow impedances arranged as shown in Fig. 2, 
water is flowing. The following pressure conditions are in effect: pp= 1.2 atm, py = 1 atm, and p, = 0 atm. With 
these values of po, py , and pp, formula (5) is reduced to the form of: 2,241 2— pi! Pi =(T). Thecurve 
expressing the function @ (p,) = 2,24V'1,2—>, / p, is shown in Fig. 6b. Chosing a number of different values 
of T, we take the values of »(T) corresponding to these values of T from curve 4 of Fig. 6a. Then equating these 
values with the values of y(p,) expressed by the curve of Fig. 6b, we find the values of p; corresponding to the chosen 
values of T, and thus obtain the curve p, = 9(T) shown in Fig. 6c. We see that with increase of temperature the 
value of p; decreases, contrary to what is taking place in gas flow with the same arrangement of flow impedances. 
This is explained by the pointed out above opposite effects of temperature upon viscosity in gases and in liquids, 
When the arrangement of the flow impedances is reversed (the laminar impedance in the inlet, and the turbulent 
in the outlet), then for a liquid the value of p; must increase with temperature. 
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Il, Consideration of Errors Caused by Heat Exchange with the Walls, and Their 
Compensation 

1. Conditions for decrease of errors connected with heat exchange, Let us estimate the temperature change 
of the working substance caused by heat transfer to it from the walls. We shall examine flow of liquids; however, 
the results of our examination can be applied also for an approximate analysis of conditions in the case of gas flow. 


A nearly perfect turbulent flow impedance can be made in the shape of an orifice in a thin plate; and the flow 
through such an orifice is practically adiabatic. More substantial processes of heat exchange take place in the laminar 
type flow impedance and in the chamber between the two flow impedances. We shall limit our examination to the 
heat transfer from the walls to the liquid in a laminar type impedance having a channel of diameter d and length J. 
The conditions of heat exchange in the chamber can be evaluated in a similar manner. 


Let us analyze the case when the walls of the flow impedance are heated to a constant temperature T,,, and 
the liquid entering the channel of the flow impedance is at a temperature T. As a result of heat transfer from the 
walls, the liquid, while flowing through the channel of the impedance with velocity w, is heated, and its temper- 
ature increased by 5T. Our object being merely a rough estimate of 5T, we shall undertake to find only the con- 
stant mean values of the specific gravity y, and of the specific heat c of the liquid passing through the channel; 
we shall consider the liquid in the channel to have a mean temperature of T +1/25T. We will show that under 
the above conditions the quantity 6T/ (T,—T) can be represented as a function of only two parameters: the ratio 
i/d and the Reynolds number Re, which completely determine the flow conditions. 











From the equation of heat balance: 


andi[T —(T+ )|=ew8P yer 


we obtain the following: 





ar 2 
y —7 "-"s 2 (6) 
* ‘+ Ya T° 


In the above formula, a is the coefficient of heat transfer, which in the case of laminar flow can be represented 
in the form of the following function of Re, c, y, and w: 


a = 3 opw 
7 her 


This has been derived from the general relationship between a and the coefficient of friction losses f when 
a liquid flows through a channel with heat exchange, which is: a = ot cyw (see [10]. 


Substituting the above expression for « in terms of the Reynolds number into (6), we obtain the following simple 
formula for determination of §T/(T,—T) from the given values of I / d and Re: 


or 
TT 





(7) 


P=) 
@ 


1+ 


B+] 0 
~~ 
xl 
~ 


It follows from (7) that to reduce 5T to a minimum, the values of Re used in these operations must be as close 
as possible to the upper limit Rejjy, of the Reynolds number beyond which laminar flow is not possible. On the basis 
of (7) it would also appear that 1/d must be chosen as small as possible. This seems to be feasible, since to secure 
laminar flow at small ratios of 1/ d, it is only necessary to have a sufficiently small d*. The desirable characteristics 


* It follows from the conclusions in paper [8] that the limiting value of d at which laminar flow in an impedance 


changes to turbulent flow is d = 270 va yt » Where yu is the dynamic coefficient of viscosity, p —the den- 
Pp 
sity, and &p-—the pressure differential across the impedance; all quantities are measured in units of the kg-m-sec 


systems. 
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can be obtained by using for the laminar type impedance a honeycomb construction with a large number of parallel 
channels. A flow impedance made of a porous material comes very close to a honeycomb flow impedance in its 
characteristics. 


However, when the value of 1 / d becomes very small, in the present application another consideration enters 
the picture: when //d is very small, the process of impeding the flow is carried out not so much by friction losses 
in the channel, as by localized losses of the total pressure at the entrance to the impedance device and by losses 
at the exit from it; and when this condition is reached, then the dependence of material throughput upon temper- 
ature which was the basis of our deductions in section I above is no longer valid. Therefore, actually the relative 
temperature change of the fluid under test while passing through the laminar flow impedance cannot be made very 
small, For instance if we set in formula (7) Re = 2000, and 7 /d = 2, the quantity 5T/(T,,—T), will amount only 
to about 0.03; but, if Z/d has practically acceptable values (of the order of 10 to 20), the value of this quantity goes 
up already to from 0.15 to 0.28, respectively. It is true that these figures are somewhatexaggerated due to the fact 
that in the deduction of formula (7) it was assumed that Ty is constant, while in reality under nonstationary con- 
ditions, when the passing through fluid is warmed up, simultaneously also the temperature of the walls is decreased. 


Taking into account the above relationships, it is possible to decrease the adverse influence of heat transfer 
upon obtaining a desirable shape of the function p; = ¢(T) for a test chamber of the described type. However, as 
follows from the above discussion, to reduce this influence to a completely negligible extent merely by proper dimen- 
sioning is not possible. Therefore, the method of measurement of nonstationary temperatures examined in this paper 
will be advantageous only, provided a measurement procedure can be developed by which the multiple~valuedness 
of the function p; = 9(T) due to heat exchange can be eliminated, and the effect of heat transfer as a factor slowing 
down the process can be excluded, In the following, two measurement procedures are presented which fulfill these 
conditions. 


2. A method of practically inertialess discrete measurements of instantaneous temperature values in gases and 
liquids, The first necessary condition that must prevail, if temperature measurements using the principle exposed in 
section I of this paper are to be carried out, is the removal of the multivaluedness of the function p, = ¢(T), which 
is caused by the fact that in a nonstationary process at every given instant the temperatures of the conduit walls 
may be different, even though the fluid passing through the inside is at a certain given temperature. In order for the 
function p; = 9(T) to be strictly singlevalued in the presence of heat exchange between the fluid and the walls, it 
is sufficient that the walls be at a certain fixed temperature at the start of every measurement, and that the entire 
measurement be completed in a certain predetermined time. The temperature of the walls can be maintained prac- 
tically constant, if the measurements are carried out not continuously, but in such a way that the flow through the 
chamber of the fluid under test is switched on only for short intervals of time, and sufficiently long intermissions 
are left between these intervals for the wall temperature to be restored to the fixed initial value. This restoration 
may be carried out by an external heat source or sink, For instance, in the setup shown in Fig. 1b, the external 
jacket surrounding the chamber may be separated from the test fluid flow, and a cooling liquid passed through it. 
However, it is more advantageous to use internal flushing of the test chamber by a cooling substance, and if the 
flushing is done in reverse to the direction of the test flow, then in addition to the cooling, the function of cleaning 
the chamber and flow impedances of contaminating particles is carried out. Under these conditions, the character- 
istic function p; = y(T) can be computed with precision taking into account the thermal properties of all used ma- 
terials, This function can also be made more precise by laboratory calibration of the sensor with the aid of optical 
temperature measuring devices or high sensitivity thermocouples, which are unsuitable for continuous operation under 
field conditions. Such a calibrated pneumatic or hydraulic sensor can be used for measurement of instantaneous tem - 
perature values under field conditions, By using several such sensor chambers, switched on one after another, the 
intermissions between successive measurements can be made as short as desired. 





The offered here method of temperature measurement is of interest in connection with the general develop- 
ment of discrete, step by step techniques in control and management of operations. 


3. Compensation of errors connected with heat transfer in continuous measurements. As was already pointed 
out, in the case of continuous measurements, the transient wall temperatures can be different for the same temper- 
atures of the fluids being tested. However, in this case also, the measurement errors connected with heat transfer 
can be excluded. For this, it is only necessary to make a proper selection of the material of the sensor chamber or 
of one of the flow impedances: with this done, when the walls are subjected to a warming up by a heat transfer 
which would change the value of p;, a change of crossection of the passage area of one or both of the flow im- 
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pedances by heat expansion or contraction will occur which will compensate the effect of the heat transfer. For 
example, since the cross section of the passage through flow impedance 2 (Fig. 2) increases with increase of the wall 
temperature, the increase of py which would take place because of heat transfer from the walls to the test fluid if 
the passage area would be unaltered, is compensated for by a decrease of p; because of the increase of passage area 
of the given impedance due to heat expansion, In practice it may prove more convenient to carry out this com- 
pensation in the turbulent flow impedance. 


Both opposing processes, the temperature change of the walls affecting the heat exchange with the fluid under 
test and the change of passage cross section of the flow impedance, proceed slowly. However, in view of the fact that 
these two factors are interconnected and act simultaneously but in opposite directions, for the fast process represented 
by the function p; = y(T) the compensation can be complete if the parameters are chosen correctly. 


In conclusion, we will point out that in this paper only the principle of this method of temperature measure- 
ment of gases and liquids was described, and the feasible design modifications of suitable pneumatic or hydraulic 
transducers were not analyzed. For instance, depending on the range of temperatures to be measured, there will 
have to be a choice made of suitable materials for the components. To obtain sufficiently high material-throughput, 
it may be advantageous to install a honeycomb cell for the laminar flow impedance, or the part of the laminar im- 
pedance may be carried out by the walls of the chamber if they are made of a porous material. It may prove ad- 
vantageous to install a special compensating element (for instance, a needle, inside of the flow impedance channel) 
for the above discussed compensation of errors. Special consideration must be given to the problems of filtrations 
etc. All these problems must be solved in connection with the actual technical conditions. 
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THE STABILITY OF A CLASS OF SYSTEMS WITH VARIABLE 
PARAMETERS WHICH VARY PERIODICALLY AND STEPWISE 


I, V. Pyshkin 


(Moscow) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 9, 
pp. 1244-1247, September, 1961 

Original article submitted December 30, 1960 


A characteristic equation is derived from which the stability of a class of variable parameter 
systems is determined, the parameters varying periodically and stepwise. The equation can 
also be used to investigate the stability of periodic oscillations in nonlinear systems with sim- 
ple nonlinearities, 


The systems with variable parameters which vary periodically and stepwise were investigated in a number of 
papers [1, 2] and the conclusions were drawn that the investigation of such systems in their general form leads to 
extremely involved relations. In [8, 4] special cases were considered of systems with variable parameters, in par- 
ticular, systems with a switchkey which can be regarded as ones whose feedback coefficient varies periodically and 
stepwise. 


In the present paper a method is given which is based on the results obtained in [3] when investigating the 
stability of a class of systems with variable parameters, the latter varying periodically and stepwise. The charac- 
teristic equation which determines the stability of this system with variable parameters, that being the most en- 
countered case, is obtained in a compact form. 


The obtained results can also be used in the investigation of the stability of auto-oscillations in systems which 
contain simple broken -line nonlinearities, in particular, those with saturation. 


The system of simultaneous equations describing the behavior of a linear control system in the absence of a 
disturbing force, can be given in the form 


dq hn 
sr = >) Ove () Pe (k= 4,2, ...,m). (1) 


a=] 


In the case of constant coefficients, in order that the corresponding system be stable, it is necessary that the 
roots of the characteristic equations 


P—bn = Os i da — bin 
sas a al 
Ter ik aevtget nation .. |= @) 


have negative real parts. 


We shall consider a class of systems with variable parameters, the latter varying both periodically and step- 
wise; in these systems the increment Aby,q of an arbitrary coefficient by, can be represented in the form of a prod- 
uct of two sequences of constant quantities hy, and y, (k, a = 1, 2,...,m) such that 


Ab, nos hye ’ @) 
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i.e., 


bee for aT <Lt< aT +7i, (4) 


Ora = | 
bea +hyte for nT +71 St<CnT +71 + Ts = (n + 1) 7. 


The system of differential equations (1) with variable parameters can be written in this case as 


a D PnaPa th (hy Dd) tee (R= 1, 2,-..., m), ' 
a=} a=} 
where 
for aT <t<nT + 71, (6) 


0 
r0-| 
i for aT +7Ts<t<nT+T. 


It was shown in [3] that the characteristic equation which determines the stability of the system (5), takes at 
the discrete time moments 0, T, 2T, . .. the form 

















ePT LL hhTr tT ePT  phaTr tits ePT gm Tet iTs 
Hai — An Ha — In iid Hm — Ma 
ePT ebaTretAeT . eT bad ebsT i tAaTs ePT a ebmTrtATs 
T > os =0, 7 
shales wun” in— Ta Him a 9 


oe 2 a ees 6S Oe err eee ere. “65160 82 658 2.0e¢ © © © 6 


PT Trl: PT ptsTit Ts ePT _ ghmTrt}AmTs 
Hi — Am Ha — Am mE tray 














where 4, Ha, +++ Hyp are the roots of the characteristic equation corresponding to the system (1) for the first values 
of the parameters, or to the system (5) when f(t) = 0, that is they are the roots of the equation(2); Ay, Az, -- + + Am 
are the roots of the characteristic equation corresponding to the system (1) for the other values of the parameters, or 
to the system (5) when f(t) = 1, that is they are the roots of the equation 


4 —fbu — ant — bin — Aya ui — dim — Pim 
— ba — A — beg — ait —b. — 
SM crinccatagt Sr ials, peel Pt {mo ) 
— by — Amt — bing — hts A — Din — PenTen 


The characteristic equation (7) is essentially derived as follows [3]. The initial form (5) of the equation is 
reduced to A. I, Lure canonical form for the time intérval nT < t < nT + T; and independently also for nT + T; < 
<t<nT+T, A system of difference equations is found between the canonical variables at the discrete instances 
nT and (n + 1)T; these difference equations are obtained by adjusting the solutions of the canonical equations in the 
above time intervals also taking into account the relations between the canonical variables. The characteristic 
equation (7) is now obtained in the usual manner from the system of difference equations. 
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In order that the system with variable parameters be stable, it is necessary that the real parts of the roots p;, 
Pa» + - » » Pm Of the equation (7) be negative. 


Thus, if the original system (1) with variable parameters can be brought to the form (5), the following steps 
should be taken when investigating its stability: 


1) to find the roots of the system characteristic equation for one set of the values of the parameters: 4 ;, M2, 
eee Lm: 


2) to find the roots of the characteristic equation for the other set of the values of the parameters : A,, Ag, 
“ees Ami 


3) to substitute these values into the characteristic equation (7) given in its determinational form; 


4) to represent the characteristic equation (7) in the form of a polynomial of the variable ep? by expanding 
the determinant, and subsequently to apply any of the stability criteria for pulse automatic systems as given in [5}, 


The restricting condition (4) is satisfied in a wide range of problems of practical interest. This is particularly 
true for closed-loop systems consisting of a constant parameter part of arbitrary order in cascade with an aperiodic 
member whose parameters, the gain and the time constant behave in the periodic stepwise pattern, 


The results obtained can also be applied in order to study the stability of auto-oscillations in a nonlinear system 
with a broken line nonlinearity N, as shown in Fig. 1. It is known [6, 7] that the investigation of stability of auto- 
oscillations in systems of this kind can be reduced to the study of stability of the equilibrium position of a system 
with periodically varying parameters. In such a case, the new system is obtained by replacing the nonlinearity N 
by an amplifier whose gain assumes periodically the values tan a and tan 6 (with the period equal to the period 
of auto-oscillations), Such an equivalent system can bedescribed by simultaneous equations which can easily be rep- 
resented in the form (5). 


EXAMPLE 
We shall consider a typical automatic control system consisting of one integrator and one aperiodic member 
(Fig. 2). The equations of the system (in the absence of a disturbance) are 


To RQ +H=EOO TH, =—M 


%, = ri (r+ (OG, Gy = ker, (9) 


- i k (t) i 
> pr l(t =— FH a= hip ae ees 


Let the time constant T(t) and gain k(t) vary periodically and stepwise. We denote the increments of the 
parameters p, and k, due to the increments of the parameters Ty and K by Ap, and Ak;. By putting 


where 


m=An, 12 = Ak, hy =1, hy=0, 


we obtain the system (9) in the form (5). Consequently in view of (7), the characteristic equation can, in this case, 
be written in the form 














ePT eit rthiTs ePT eT Ts 
ePT _ hTrtuTs PT LL eT hats =0, 
ai — Aa Ha — Ma 


where 1; and {lz are the roots of the system characteristic equations for one set of values of the parameters and for 
the other set. 


We put T = 1 sec, T, = Tz = 0.5 sec, T, = 1 sec. 


1115 








k() = va for n<tcn+ 0.5, 

1.2 for at05<t<cat i, 
r= {5 for n<t<at0.5, 

0.2 for n+ 05<t<ant i. 


Then for n < t < n + 0.5 the characteristic equation becomes 


Oip +1 mont 
| ; 


‘ = O.ipt+p + 0.9 = 0, 
that is,=-1, Mg = —9. 


For n+0.5<tene+1 
pal —1i,2 


= 4 = 
{ h 0.28 A+ 1.2 = 0, 


that is Ay = —2, A, = -3. 
By substituting these numerical values into (10) we obtain the following characteristic equation 














eP — ¢g0-5-1 eP — e451 
oP — ¢70-6-1.5 oP — ga h-b-15 
—TF3 —o+3 


4e?? — 1.158 e? + 0.00222 = 0. 


In order that the roots of the equation 
age*?T +4 aye?T 4.4, = 0 


be negative the following inequalities must be valid [5]. 


@+a;+a>0, ag—ai+a>0, a—a>0. 


In our case a = 4, ay = —1.158, ay = 0.00222. This implies that the varying-parameter system of our example 
is stable, 
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DETERMINING THE TRANSFER FUNCTIONS OF MULTI- 
DIMENSIONAL LINEAR SYSTEMS FROM THE STATISTICAL 
CHARACTERISTICS OF THE INPUT AND OUTPUT 
QUANTITIES OF THE SYSTEMS 


J. Maty4s and J. Silha4nek 


(Pardubice, Czechoslovakia) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 9, 
pp. 1248-1252, September, 1961 

Original article submitted January 10, 1961 


The method described in this paper makes it possible to find a mathematical description of linear 
systems from the spectral densities of the random input and output processes. If internal noise 
sources are present in the system, then the proposed method can be used to determine the spectral 
density matrix for the noise at the system output. Thus itis possible to determine the effect of in- 
ternal noise sources on the behavior of the system. 


In investigating automatic control systems using simulators it is necessary to know the transfer functions or the 
differential equations for the system. 


In practice it is often required to investigate systems whose mathematical description is unknown. Finding 
such a mathematical description is in many cases a rather complex problem. 


The input and output quantities for such systems can be treated as stationary random functions. The method 
described in this paper can be used to determine the transfer functions of the investigated dynamic systems from the 
statistical properties of the input and output quantities. 


If the investigated systems are not rigorously (exactly) linear and if their exact behavior cannot be described 
simply, then the proposed method makes it possible to find the linear approximation for these systems. 


Statement of the Problem 








Systems with several inputs and outputs shall be called multidimensional systems. Figure 1 shows a multi- 
dimensional system with n inputs and m outputs. Assume that stationary random processes (a random vector) act at 


the inputs of this system: 


v = [v, Vag «9 v,]- 


Assume that the system contains internal noise sources which are uncorrelated (for example, statistically in- 
dependent) with respect to the input processes. Then the output random vector 


u=([u, Ua, e+e Um!) 
can be described by a sum 


user+y= (2+, a+ Yay «+s lm + Ym) (1) 


where x are the output components without the noise and y are the components produced by the internal noise sour - 
ces in the system. Thus the investigated system can be represented in the form of the equivalent block diagram 
shown in Fig. 2. 
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We shall use 
S = [S,,) (i, kori, 2, ..., a) (2) 
to denote the spectral density matrix for the input quantities (the input random vector) v. The expression 
G =[G,,] (i, k= 1, 2, ..., m) (3) 
shall denote the spectral density matrix for the output quantities (the output random vector) u, and 
H =[Hy) (i=1, 2, ..., a) h=4, 2, ..., m) (4) 


denotes the matrix for the mutual spectral densities of the random vectors v and u (cf. [1)). 








aw 
“~ 

| a 
fay “Se 
3 




















v, o—: }——ou,, ye— fo TF] 
Fig. 1. Multidimensional system of the Fig. 2. Equivalent block diagram for a mul- 
general type. tidimensional system containing internal 
noise sources. 


Assume that the investigated system is a multidimen- 
sional linear system of the general type. A portion of this 
system (that corresponding to the kth output) is shown in 





Fig. 3. 
The transfer functions Y;), of the investigated system 
form the matrix 
Fig. 3 
¥=[Y,] G@=1, 2, ..., kei, 2,...,0). © 


Our problem now consists of the following. 


1, We must determine the matrix Y for the transfer functions of the investigated system C from the specified 
spectral density matrices S, H, G. 


2. We must establish whether or not internal noise sources exist in the system C. 


3. In the case y » 0 we must determine the spectral density matrix for the vector y and construct shaping 
filters for generating the random vector y (the random quantities y;, yz, . . . » Yp, which are the components of the 
vector y). 


The Fundamental Matrix Relationship 

Since the random vectors v and y are uncorrelated, it is not difficult to prove that the matrix H (4) is also the 
matrix for the mutual spectral densities of the vectors v and x. Therefore (cf. [1], formula (48)] the matrices S, H, 
Y are related by the expression 





SY = H (6) 
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If the random quantities vj (the components of the input vector v) are not mutually correlated, then the 
matrix S is a diagonal matrix. 


Therefore in this case the transfer functions Y;), are determined very simply from (6): 


Hy 
Yunge. (7) 


In the general case the matrix relationship (6) represents m systems of linear algebraic equations with a general 
matrix S for the system and n * m unknowns Yij- The system (6) has a unique solution only in the case where the 
determinant for the system is non-zero (i,e., when the matrix S for the system is non-singular), Since S is a Her- 
mitian matrix, it follows that the uniqueness condition for the solution of (6) reduces to the form 


|}S|>0. (8) 


In this paper we shall limit ourselves to the case of a non-singular matrix S [i.e., we shall assume that con- 
dition (8) is satisfied]; In that case there exists an inverse matrix s~! and the solution of the matrix equation (6) is 
determined from the formula 


Y=s "nH (9) 


For the case of a singular matrix S(| S| = 0) either Eq. (6) is contradictory, or an infinite set of solutions 
exists. The transfer function Y;; of the investigated system cannot be determined in single-valued fashion in this 
case. Therefore this case is not of interest as far as practical applications are concerned. 


The Solution of the Matric Equation (6) 








The matrix equation (6) represents m system of linear algebraic equations with an over-all system matrix S 
with n - m unknowns. In order to solve it we may use different methods cited in the literature (cf., for example, 
[3]). Thus by applying Cramer's rule we obtain each transfer function Yj}, in the form of a ratio of determinants. 
But the computation of higher order determinants whose elements are the functions of a complex variable requires 
difficult computations. The over-all matrix S for the systern (6) is by definition a non-singular Hermitian matrix, 
This fact permits considerable simplification of the process involved in the solution of the problem. It is not dif- 
ficult to prove that the matrix S can be written in the ~ of the product 


S = BB, (10) 
where B is a triangular matrix. The matrix S can be reduced to the form (10) using, for example, the method de- 
scribed in [6, 7]. The elements Bj, of the matrix B are, with the exception of phase filters, the transfer functions 


of the system of shaping filters which generates the random processes with the spectral density matrix S. The matrix 
B" is the transposed matrix with conjugate elements relative to B, 


Therefore the matrix equation (6) can be replaced by the equivalent system of equations: 


BX=H, BY=X, (11) 
Since the matrix B, and therefore B', are triangular matrices, it follows that the solution of matrix equations 
(11) is not difficult. First of all we use the first equation to compute the auxiliary matrix S, and then use the second 


equation to compute the desired matrix Y for the transfer functions of the investigated system C. In view of (10) it 
is possible to write the inverse matrix S~* in the form 


S-1 = B- (B’)". (12) 
Thus by substitution into (9) we obtain the solution of the fundamental matrix relationship (6) in the form 


Y = B- (BH. (13) 
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Formula (13) can be used expediently for solving the matrix equation (6), since Eq, (6) represents m systems 
of linear algebraic equations with an over-all matrix S, and the inversion of the triangular matrices B and B" is not 
difficult. 


At the beginning of this paper the problem was limited to an investigation of linear stable systems whose in- 
put and output quantities are stationary random processes. We also adopted condition (8) which requires that the 
matrix S must be non-singular. 


Assuming that these conditions are satisfied, it is possible to use the method described above to determine 
the matrix Y for the transfer functions of the investigated system. This matrix describes the behavior of the in- 
vestigated system. This matrix describes the behavior of the investigated system exactly. Such a statement is valid 
only on the assumption that all of the measurements and computations which are used to determine the spectral den- 
sities Sj, and Hy, are performed exactly, and that the computations used to solve the system (6) are also formed ex- 
actly. In practice it is necessary first of all to use a special computer (correlator) to find all the required correlation 
functions; then some approximate Fourier transform method is used to compute the corresponding spectral densities. 
These measurements and computations cannot be performed exactly in practice, Thus, the limitations cited above 
will not be rigorously satisfied. Therefore the solution Y in this case determines a certain linear approximation of 
the mathematical description of the investigated system. 


Internal Noise Sources 
If the system C contains internal noise sources, then it is possible to determine the matrix for the spectral den- 
sities of the random vector y whose components y; correspond to the noise signals at the system output. 





By definition the vectors v and y are uncorrelated; i.e., the matrix for the mutual correlation functions will be 
a zero matrix (cf. [1]): 


M {v(t)’ y(t + t)} =0 (14) 
(here M denotes mathematical expectation). It is not difficult to demonstrate the fact that the matrix 
M {x (t)' y(t +))} =0, (15) 


will also be a zero matrix; i.e., the random vectors x and y are also uncorrelated. 


In this case the matrix G can be written in the form 
G=6,+6,, (16) 


where G, and Gy are the matrices for the spectral densities of the random vectors x and y. 


The matrix Gy, can be written as (cf. [1]) 
G, = Y’SY. (17) 
The matrices S, G, Y are assumed known, and therefore we determine the matrix Gy from the formula 


int 18 
G,=G— G,=G—Y'syY. (26) 

Formula (18) determines the spectral density matrix for the random vector y corresponding to the noise sig- 
nals at the system output. The shaping filters which generate the quantities y; with the specified spectral density 
matrix Gy can be determined using the methods described in [6, 7}. 


Example 


It is required to determine the matrix Y for the transfer functions of a linear system with two inputs and three 
outputs from the specified spectral density matrices: 
1 


i 
4—s8 (i—s) (2—s) 
i 2 y 
(+s) (2+s)i—s 





S= 
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1 1 2 
na| @+9@—9) 4—#BFHR—5 
a ‘ ‘ ‘ rf 
(1 + 8) (2 + s) (1+-s)*(3 + 8) (1 + 8) 
First we determine the matrix B by the method described in [6]: 








1 | 
2+s i-s 
i 
i+s 


B= 





0 


The transfer functions Y;, could be computed by solving the system of matrix equations (11). Here it is ex- 
pedient to use formula (13) to solve the problem. 


We shall compute the reciprocal matrix 





i—s 


i-+s 
B-les 2+ 8s — (2 +8) 
0 i+s 
and from formula (12) we obtain 


atin 2(4—s) — (1+ 5) (2 + 8) | 
7} csesoyeraree oe ; 


Substituting into (9), we find 














3+s s 3 (2 +s) 

i+s i+s 3+8 
i—s i—s i—s |’ 
—W+9)2+5 G+)e+8 —3Fe 


and thus we have determined the mathematical description of the investigated linear system. 


In the case where the spectral density matrix G for the output random vector u is also specified, we first of all 
establish whether or not internal noise sources exist in the system. If they exist, then from formula (18) we deter- 
mine the spectral density matrix Gy for the random vectory y (the output noise signals of the system). The shaping 
filters which generate the corresponding random processes can be defined and constructed by the methods described 
in (6, 7]. 


SUMMARY 
The paper describes the method for determining the matrix for the transfer functions of multidimensional 
linear systems from the spectral densities of the random input and output quantities. These transfer functions are 
the solution of the matrix equation 


SY=H 


i.e., they are a solution of m systems of linear algebraic equations with an overall system matrix S. 


If the system contains internal noise sources, then the method described above can be used to compute the 
spectral density matrix for the random vector y (the output components of this vector correspond to the internal 
noise sources in the investigated system). 


In practice this method makes it possible to compute the linear approximation of multidimensional systems 
whose mathematical description is unknown. 
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Original article submitted August 15, 1960 


The present article is concerned with the effect of Coulomb friction and of the feed pump 
capacity limitations on the stability of duplicating -machine hydraulic servosystems. The 
solution which was obtained by means of a mock-up of the device showed that the previ- 
ously used solution, which was given by A. A. Andronoy and A, G. Maier, is incorrect in 
this case because of the difference between the initial conditions. The solution obtained 
here makes it possible to render the known methods for calculating hydraulic servosystems 
more accurate. 


Coulomb friction in the guides of operating parts exerts a considerable influence of the stability of hydraulic 
duplicating machines. An attempt was made in [1] to take into account the Coulomb friction forces in deriving the 
stability conditions for duplicating machines. In this, the results obtained by A. A. Andronov and A, G. Maier [2] 
in applying the direct control theory were used. 


E. M. Khaimovich [3] and T. M. Bashta [4] later used the same results. However, further investigations showed 
that, strictly speaking, the direct application of the results obtained by A. A. Andronov and A, G, Maier to the prob- 
lem of the stability of a hydraulic servosystem was not justified due to the difference between the initial conditions, 
although the obtained stability condition was in many cases in good agreement with the results obtained in practice. 


In connection with this, it became necessary to analyze this problem more rigorously and to formulate the 
stability conditions in a more precise manner. For a more accurate description of the system's behavior under os- 
cillating conditions, the characteristic of hydraulic drive saturation, which is determined by the finite capacity of 
the pump feeding the system, was introduced in the consideration. An MN-M electronic simulating device was used 
in solving this problem. 


For checking the accuracy of the obtained solution, we used the same mock-up for producing the solution of 
Vyshnegradskii's problem, which coincided with the exact solution given by A. A. Andronov and A, G. Maier. 


The derivation of the basic equation of motion is given in [1], Here, we shall consider only the effect of 
Coulomb friction and of the saturation characteristic, which is determined by the pump capacity limitations. 


The linearized equation of the system for a linear input function, i.e., for a constant velocity of the pickup 
end-piece, where the friction force is not taken into account [1], has the following form: 


a dB db -~ pay (1) 
LEV +MI4+D5 +b =Dz+R. 

Here,L is a coefficient which characterizes the effect of the compressibility of the operating liquid and of the 
elasticity of piping, M is the mass of the machine's operating element, D is the internal damping coefficient, E is 
the rigidity coefficient, R is the constant component of the cutting force, 6 is the error in duplicating, and y is the 
input coordinate (the coordinate of the pickup end -piece); 

6 = y-x, 


where x is the coordinate of the machine's operating element. 
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The problem of the stability of this system can be reduced to Vyshnegradskii's problem in the theory of direct 
control. The Vyshnegradskii parameters for the system under investigation are 


M (2) 


Gay". 





VEL VEE 


The stability condition is given by the inequality UZ = 1, which corresponds to the points that lie in the plane 
of parameters U and Z above the Vyshnegradskii hyperbola and on it. 


Let us consider the system's behavior in the presence of Coulomb friction in the case where it moves about the 
equilibrium position. 
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Fig. 1. Simulation scheme. 


Assuming that R = 0 and y = 0, the system's behavior can be described by the following system of equations: 


ax , Lap dz 
Moa =P+ile’). op Gr tP=—ke D>, (3) 


where P is the motive force; 


f Pr 
+> for *<° oF 
P, 
for * =Vand PI< >": 
P 4 
te") = Ew ose... (4) 


P, 
for 2’ =Oand P|>—-> 





, P, 
iP for Onl PIS 
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For the sake of brevity, we shall denote this function by the symbol 
P, 
{(z’)=— > sign z’. 


By dividing the second equation in system (3) by E, we obtain: 


L dP 4 D dz 


WEatE?"-*-ta° 


We shall now reduce the system under investigation to the dimensionless form by substituting the variables 
according to the method used by Vyshnegradskii and by A. A. Andronov and A. G, Maier: 


3/ kh. yl (5) 
z=, P=1M V P,=m yf #2, t= Et 


Here € and » are the dimensionless coordinates of the system, and r is the dimensionless time. 
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Fig. 2. a) Boundaries of the stability region for different initial conditions for the case con- 

sidered by A. A. Andronov and A, G, Maier; b) boundaries of the stability region of the dup- 

licating system for different initial deflections of the pickup valve. 1) &»= 0.5; 2) E, = 1.0; 
3) £9 = 2.0; 4) Eg = 2.5; 5) & 9 = 5.0; A) Vyshnegradskii hyperbola, 


After performing the substitution and the necessary simplifications, system (3) assumes the following form: 





i d ° 
= =O— = 7 Sign = 2, + Uy = —t—s&, (3°) 
where U and Z are Vyshnegradskii's parameters. 
d 
Assuming that - » we obtain: 

S on = = 
a7 =% > sign n, (3") 
ep ne 


The initial conditions are: E59 #0, 19 = 0. 
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Due to the fact that the capacity of the pump that feeds the system is a finite quantity, the velocity of the 
copying support cannot exceed the value 


, Q 
* max E F : (6) 


where Q> is the pump capacity and F is the cross-sectional area of the piston, if the condition of the liquid flow 
continuity is to be preserved. 


In dimensionless form, 


Q " 
(FE) naa "max — FE (7) 


tlh 


Thus, the following additional condition is imposed on the solution of system (3"): 
in|} S9 max (8) 


In order to solve the system of equations (3") by taking into account condition (8), we composed the simula - 
tion scheme shown in Fig. 1. 


The variation of the U and Z parameters in this system is realized by means of potentiometers 1P and 2P. 


The simulation method consisted in letting the U parameter change to a value for which the phase trajectory 
projection on the plane n = 0 became a closed curve for a given value of the Z parameter and for the assigned 
initial conditions. 


For larger values of the U parameter, the phase trajectories “converge” toward the equilibrium zone, which 
corresponds to a damped process; for smaller U values, these trajectories "flare out,” which corresponds to a divergent 


process. 
The thus obtained parameter values are given in Fig. 2b in the shape of a family of curves. 


Curves 1-5 represent the stability boundaries of the system for different initial deflections of the pickup valve 
from the neutral position. For the sake of comparison, this figure also shows the Vyshnegradskii hyperbola (Fig. 2b) 
and the stability boundaries for different initial conditions for the case which was considered by A, A. Andronov and 
A. G, Maier (Fig. 2a). 


According to condition (8), only a limited amount of energy from the outside (from the pump) can be supplied 
to the system. The energy supplied to the system is regulated by the magnitudes of the cross-sectional areas of 
the slide valve passages. Beginning with a certain magnitude of the passage cross section, which is determined by the 
maximum displacement |€ 4x !, the supply of energy does not change, and, consequently, if the system is stable 
for |€ max!» it will also be stable for|&|>|£,.x/. | Such a system can be considered as absolutely stable. 
The value of | max! is given by 





-_ 9 
al ati Uy max “7 
By expressing the relative parameters in terms of the absolute parameters, we obtain: 
M EQ \ 
l€ nay! = ( Pip}, (10) 
max P, VDE FCo max 


where, according to the linearization conditions for the discharge equations [1], 
P nal 0,5 pF, 


where p is the pump pressure. 
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As was shown by mock-up tests, the amplitude of self-oscillations which can arise in the system if it is outside 
the stability boundaries is directly dependent on the pump capacity. 


SUMMARY 
The results obtained by A. A. Andronov and A, G. Maier cannot be directly extended to the problem of the 
stability of hydraulic servosystems in duplicating machines by taking into account Coulomb friction in the guides. 


As a result of investigations of a hydraulic servosystem where the limitations with respect to the pump capacity 
were taken into account, we have derived the absolute stability conditions for the system in dependence on its para- 
meters and the pump capacity, which made it possible to formulate the previously derived stability conditions with 
a greater degree of accuracy. 
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1127 








EFFECT OF LEAKAGES ON THE CHARACTERISTICS 
OF PNEUMATIC FORCE-COMPENSATING UNITS 


A. 1, Gitel’man and V, M,. Syrodoev 


(Leningrad) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 9, 
pp. 1257-1261, September, 1961 

Original article submitted February 23, 1961 


This article presents theoretical and experimental data on the operation of a membrane -nozzle- 
valve complex under air pressures of up to 10 atm. Recommendations concerning the stability of 
the characteristics are given. 


Along with pneumatic automation devices which operate under air pressures of up to 1 atm abs (AUS devices, 
etc.), devices using pressures of up to 10 atm are also used to a large extent. Asa rule, the use of high air pressure 
is determined by the characteristics of the systems to be automated. 


The force-~compensating unit is one of the basic units in modern pneumatic automation devices. The com- 
plex consisting of a membrane, a nozzle, and a valve represents a widely used basic design of this unit. The opera- 
tion of this complex at high air pressures is connected with a number of peculiarities which make it different from 
operation at low pressures. 


The present article provides the generalization of some experience gained in factory adjustment of automatic 
control systems with force-compensating units that operate under a pressure of 10 atm. 


Figure 1 shows the schematic diagram of the force-compensating unit. Its operating principle is based on the 
compensation of the command force Q, by the force Q,,, which the air pressure p in chamber A exerts on the mem- 
brane.” 


Let us consider the equilibrium of the forces acting in the force-compensating unit. 


The possible force distributions can be reduced to three cases. In the first case, the nozzle is pressed against 
the valve, and the valve is not in forced contact with the valve seat. In this case, the equilibrium of forces is deter- 
mined by the relationship 


—Q —Q = 
Q ge, + @ _ 0, (1) 
where Q, is the command force, Q, is the force which acts on the valve from the side of compressed air and which 


is composed of the spring-action force and the force resulting from the difference between the air pressures acting 
on the valve, and Q,,, is the force which the air pressure p exerts on the membrane. 


In the second possible case of force distribution, the valve is pressed against the seat, and the nozzle is not in 
forced contact with the valve. In this case, the force Q, does not figure in the balance of forces, and the equilibrium 
condition is determined by the relationship 

eS 6. (2) 

In the third possible case of force distribution, the action of force Q, is divided between pressing the valve 

against the seat and pressing it against the nozzle. 


* The case where the air is not used by a device is considered. 
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The equilibrium condition is determined by the relationship 


Here, Q, is the force with which the valve is pressed against the sealing edge of the seat. 
Let us consider the conditions under which the above cases of force distribution can arise. 


For an ideal airtightness of the contact between the valve and the seat and between the valve and the nozzle,* 
the force distribution that is connected with an increase in the force Q, within the range 0 < Q. < Qyp (section OC 


in Fig. 2) corresponds to the third case of force distribution: Q,, is spent on overcoming the force Qy» with which 
the valve is pressed against the seat. 
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Fig. 1. Schematic diagram of the force-compensating unit. 


Due to the ideal airtightness of the valve—seat contact, the air under pressure is prevented from entering cham- 
ber A along the section OC, and the air pressure p is equal to zero. 


If the force Q,, further increases (section CB), we obtain the first case of force distribution: Under equilibrium 
conditions the nozzle is pressed against the valve with the force Q,,** while the force which presses the valve 
against the seat is Q, = 0. In this case, the increase in pressure p 
occurs as a consequence of feeding compressed air into chamber A 

AyD po through the clearance between the valve and the seat, which is 
max e formed at the instant of time when the force Q,, increases since the 
relationship Q¢ > Qrm + Qy holds until equilibrium is established. 


At the point B, pressure p attains the value pyyax = Ps, and 
the subsequent increase in Q,, (the section to the right from the point 
4 I! B) is spent on compressing the spring. Over this section, the distribu- 
tion of forces which corresponds to the first case is preserved not only 
when the force Q,, increases, but also when it decreases. 


Along the section BA, we again have the third case of force 
ic distribution: A decrease in Q, is compensated by an increase in the 
——f el a. force Qs; with which the valve is pressed against the seat; this force 
vo increases to its maximum value Qgnax = Qy1- 











Fig. 2. Static characteristics of the force- fvdiin 

compensating unit in the case of ideal air- * It is considered that the contact remains airtight even if the force 

tightness of the sealing edges. that keeps the sealing edges in contact is equal to zero, 
** Since the force Qy depends on the difference between the air pres- 
sures acting on the valve, the pressure difference decreases as the pres- 
sure in chamber A increases, while Q, also decreases Qy(Qyy > Qy;). 
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In this due to the ideal airthightness of the contact between the valve and the seat and between the valve and 
the nozzle, the air cannot either enter or leave chamber A, and the pressure is p = pmax. 


If the force Q,, further decreases (section AO), we have the second case of force distribution: Under equilibrium 
conditions, the valve is pressed against the seat with the force Qy, while the force with which the nozzle is pressed 
against the valve is equal to zero. In this case, the pressure is lowered due to the leakage of air from chamber A 
into the atmosphere through the clearance between the valve and the nozzle, which appears when the force Q,, is 
reduced, since the relationship Q,,, > Q, applies until equilibrium is established. 


Thus, in the case of perfect airtightness, the force-compensating unit is insensitive when force Q,, changes its 
direction of action. The magnitude of this insensitivity is characterized by the horizontal plateau between the lines 
OA and CB (for instance, AC’), i.e., it is equal to the force Qy. 


If the contact between the valve and the seat is not airtight, the valve is pressed against the seat, however, 
compressed air continuously enters chamber A, and, when the force Q,, changes, the clearance between the nozzle 
and the valve changes in such a manner that a pressure p which secures the required equilibrium of forces is estab- 
lished in the chamber A between the throttles. 


Since, in this case, there is no dynamic interaction between the nozzle and the valve, the force Q, does not 
affect the balance of forces, i.e., the force distribution corresponds to the case (2) regardless of the direction in which 
the force Q. changes. The static characteristic for this case is shown in Fig. 3a. 
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Fig. 3, Static characteristics of the force-compensating unit in the 

absence of airtightness. a) Leakage through the clearance between 

the valve and the seat; b) leakage through the clearance between the 
valve and the nozzle; c) the leakages are commensurable. 


In the absence of an airtight contact between the valve and the nozzle, the nozzle is pressed against the valve, 
however, the air from chamber A continuously escapes into the atmosphere, and, when the force Q,, changes, the 
necessary equilibrium or forces is secured by varying the clearance between the valve and the seat. 


In this case, after the valve moves away from the seat, the nozzle and the valve are continuously in contact 
due to the pressing force Q,,, i.e., the distribution of forces for Q, > Qy (section CB in Fig. 3b) corresponds to the 
case (1) regardless of the direction in which Q, changes. 


Besides the cases considered above, we can also have the case where both the seat and the nozzle contacts are 
not airtight. The shape of the force-compensating unit characteristic will be then determined by the ratio of clear- 
ance areas which are formed in forced contact between the valve, the seat, and the nozzle. If the clearance area 
between the valve and the seat is: much larger than the clearance area between the valve and the nozzle, the charac- 
teristic will correspond to that shown in Fig. 3a. If the clearance area between the valve and the nozzle is much 
larger, the characteristic will correspond to that given in Fig. 3b. If these clearances are commensurable, the charac- 
teristic will consist of the two sections OA" and C'B (Fig. 3c). On the section OA", an increase in the force Q, leads 
to a reduction of the clearance between the nozzle and the valve, due to which the pressure p increases. The second 
case of force distribution applies here. At the point A’, the nozzle is brought into forced contact with the valve, 
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and a further increase in the force Q,, does not affect the ratio of the cross-sectional areas of passages at the inlet 
and outlet of chamber A, The force Qc is spent on overcoming the force Qy. ‘The relationship between forces 
corresponds to the third case (3). 


Along the section C’B, an increase in Q, leads to an increase in the clearance between the valve and the seat, 
due to which pressure p further increases. The relationship between forces is determined by expression (1). 


Thus, if the clearances between the elements in a force-compensating unit are mutually commensurable, the 
characteristic can have an insensitivity region, as in the case of ideal airtightmess. The magnitude of this insensit- 
ivity, which is equal to Q,, is present in both cases, however, its character is different: If the airtightness is ideal, 
the insensitivity arises when the direction of force Q, is reversed, while, when the clearances are commensurable, 
the insensitivity is determined by their ratio and is independent of the direction in which Q,, changes. 


Numerous experiments that were performed on a pneumatic command unit with metallic nozzles and seats, 
which had different hardness values, shapes of the sealing edges, and diameter ratios (for p, and p pressures of up to 
10 atm), showed that stable airtightness at both sealing contacts of a force-compensating unit cannot be guaranteed, 
Therefore, the operation of a force-compensating unit in the case of ideal airtightmess is mainly of theoretical interest. 
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Fig. 4. Static characteristics of an experimental pneumatic com - 
mand device for different leakage ratios (p, = 10 atm), The cir- 
cles correspond to air discharge G,., = 0 through the contact be- 
tween the valve and the seat and to discharge Gnoz = 60 liters/min 
through the valve—nozzle contact; squares: Gseg = 16 liters/min 
and Gpoz * 20 liters/min; triangles: G,., = 50 liters/min and 
Gnoz * 6 liters/min. Gs, values were measured for p= 0, and 
Gnoz Values were for p = 10 atm; the white points correspond to an 
increase in the command spring tension, and the black points corre- 
spond to a decrease in the spring tension. 


As a result of numerous experiments, it was established that not only the shape, but also the stability of the 
characteristics, depends on the ratio of leakages in force-compensating units. This applies especially in those cases 
where the characteristic has an insensitivity section that is located in the operating range, i.e., in the 0 < p < Ps 
range. The values of p and Q, that correspond to the appearance of insensitivity and also the slope of the insen- 
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sitivity section vary in a random fashion for different measurements without any external causes. This can be ex- 
plained by changes in the conditions which determine the ratio of leakages (clogging and division of the clearance, 
nonalignment, etc.). 


Figure 4 shows the static characteristics p = f(S) of an experimental pneumatic device for different ratios of 
leakages through the sealing contacts in the force-compensating unit. (Here, S is the deflection of the compressed 
command spring. The command force is Q, = cS, where c is the spring stiffness.) In our experiments, c was equal 
to 1.76 kg/mm, and the compressed~-air pressure was p, = 10 atm. The leakages into the seat and nozzle bypass 
were secured by suitably mounting special adjusting screws. 


If air leakages at the seat or at the nozzle are positively provided, it is possible to secure a stable linear 
characteristic of the force-compensating unit throughout the entire 0 < p < 10 atm range without any appreciable 
insensitivity. 


From what has been said above, it follows that, for a reliable operation of the force-compensating unit, it is 
desirable to provide a throttle in the seat or the nozzle bypass, whereby the existence of leakages would be guaran- 
teed. 


In choosing the throttle location (in the nozzle or the seat bypass), the requirements for the force-compen- 
sating unit characteristic must be taken into account. Moreover, the ratio of the pressure p, and p for the longest 
operating periods must be considered, since this will determine the throttle dimensions and design in dependence on 
the maximum allowable air discharge. In practice, the following considerations are of the greatest importance in 
selecting the throttle location. The results of many experiments indicate that reliable airtightness of the valve — 
seat contact can be secured by simple technological expedients (grinding, smoothing -out of sharp edges, etc.), In 
this, the airtightness of contacts is not disturbed even in operation (if the air is well filtered and the valve does not 
vibrate). As a rule, it is not possible to secure an airtight contact between the valve and the nozzle. This can ob- 
viously be explained by unavoidable misalignments of the membrane —spring group, by the nonuniform pressure 
with which the nozzle's edge is pressed against the valve surface, etc. Moreover, the escaping of the air into the 
atmosphere is determined, besides the valve —nozzle contact, also by leakages through many seals of the air duct 
between the membrane cavity and the device which uses compressed air (at the location where the membrane is 
built in, at the connecting pipe joints, etc.). These leakages can increase to a considerable extent during operation. 


Thus, if the throttle is to be installed in the seat bypass, one should strive for a greater air discharge in order 
reliably to secure the prevalence of leakages through the seat over possible leakages into the atmosphere. 


If the throttle is to be installed in the nozzle bypass, one can work with a minimum air discharge, since the 
probability that the leakages will be equalized during operation is in this case very low. Moreover, the servicing 
of the throttle then becomes much simpler, since the force-compensating unit does not have to be disconnected 
and dismantled for the checking, cleaning, etc. of the throttle. 
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INFORMATION 


Translated from Avtomatika i Telemekhanika, Vol. 22, No. 9, 
p. 1262, September, 1961 


The regular session of the IFAC Executive Council under the Chairmanship of the IFAC President, Prof. A, M, 
Letov, was held in March 1961 in Bergen. The Executive Council adopted the following resolution concerning the 
organization of the Second IFAC Congress. 


The congress will be held in Basel (Switzerland) in September, 1963. 
A. The program of the Second Congress is concerned with the following scientific fields, 
1, Automatic control theory: 
a) discrete systems 
b) stochastic processes 
c) optimal systems 
d) self-adjusting systems 
e) reliability theory 
2. Automatic control application: 
a) iivestigation of the dynamics of processes 


b) investigation of industrial automation problems by means of simulating devices and digital com- 
puters which are included in the process as well as those which do not form a part of the process 


c) application of optimizing and self-adjusting control systems 
3. Elements: 

a) new and efficient devices 

b) estimate of the reliability of elements 
4, Other subjects: 

a) education 

b) terminology 

c) bibliography. 


This program does not exclude other possible subjects; however, not less than 80% of all reports must pertain 
to the theory and application, and not more than 20% must be devoted to elements and other subjects. 


B, The total number of reports submitted to the Congress must not exceed 100. The volume of each report 
must not exceed 30,000 type-written letters, including a summary in two or three languages and illustrations. 


C, The selection of reports will be performed by the IFAC National Committees, which will be entrusted with: 
a) contacting various specialists and requesting them to send in reports on the above subjects, 


b) finding specialists who will be able to make contributions and to perform the preliminary selection of 
reports, 


¢) to inform the authors that the reports must be submitted to the National Committee not later than the 
end of 1961, and 


d) to inform the authors that the reports may be written in Russian, English, French, or German; every report 
must include a summary of not more than 200 words in the language used in the report as well as in Russian or English. 
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The final selection of the reports which will be presented to the Congress will be made by the IFAC Com- 
mittee on Reports, which will be formed by the presidents of the IFAC technical committees under the chairmanship 
of Prof. Goerecke. In this, the IFAC Committee on Reports will be guided by the following principles in making the 
decisions: 


1) the timeliness of the topic, 
2) the value of the obtained results and their novelty, and 


3) clarity of presentation. 


In connection with the above resolution of the IFAC Executive Council, the National Committee of the USSR 
on Automatic Control hereby asks the Soviet specialists to send their suggestions to the National Committee (15a 
Kalanchevskaya Street, I-53, Moscow). 


National Committee of the USSR for Automatic Control. 
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